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INTRODUCTION 


It  is  now  understood  that  the  sun’s  outer  atmosphere,  the  solar  corona,  expands  supersonically 
to  become  the  solar  wind  and  that  this  solar  wind  carries  the  solar  magnetic  field  with  it  into 
interplanetary  space.  Upon  reaching  the  Earth,  the  solar  wind  plasma  and  magnetic  field  interact 
with  the  Earth’s  magnetic  field  and  through  this  interaction,  solar  wind  energy  is  coupled  to  the 
Earth’s  magnetosphere  and  ionosphere.  Roughly  a  few  percent  of  the  kinetic  power  carried  by 
the  solar  wind  plasma  impinging  on  the  dayside  magnetopause  (between  10n  and  1012  W)  enters 
the  magnetosphere  to  drive  currents,  energize  plasma,  and  produce  a  complex  pattern  of  plasma 
convection  (e.g.,  Hill,  1983;  Baker  et  al.,  1985).  Much  of  the  extracted  solar  wind  energy  is  thought 
to  be  ultimately  dissipated  as  heat  in  the  Earth’s  upper  atmosphere,  while  another  portion  of  the 
extracted  energy  is  thought  to  be  returned  to  the  solar  wind  by  plasmoids  which  are  ejected  down 
the  magnetotail  during  intervals  of  substorm  activity. 

The  coupling  processes  which  occur  at  the  magnetopause  can  be  divided  roughly  into  two  cat¬ 
egories.  The  first  is  electromagnetic  coupling  which  involves  interactions  between  the  solar  wind 
magnetic  field  and  the  geomagnetic  field  (e.g.  magnetic  reconnection).  The  second  coupling  is  non 
magnetic  and  is  generally  referred  to  as  a  “viscous  interaction”  (e.g.  Kelvin-Helmholtz  and  impul¬ 
sive  plasma  penetration).  It  appears  that  electromagnetic  coupling  provides  the  major  momentum 
exchange  between  the  solar  wind  and  the  magnetosphere  (Lee  and  Roederer,  1982;  Podgorny  et  al., 
1978),  though  there  remains  considerable  debate  regarding  the  possibility  that  the  viscous  inter¬ 
actions  could  be  an  important  or  dominant  coupling  mechanism  under  some  conditions  [Reiff  and 
Burch,  1978;  Heikkila,  1986  and  references  therein], 

A  consequence  of  the  energy  coupling  at  the  magnetopause  is  a  system  of  field-aligned  currents 
whose  strength,  polarity  and  configuration  vary  with  the  strength  and  orientation  of  the  interplan¬ 
etary  magnetic  field  (IMF).  These  currents  provide  the  electrical  coupling  between  the  solar  wind 
MHD  dynamo,  the  magnetosphere,  and  the  conducting  ionosphere.  The  high  latitude  region  1 
currents  observed  at  the  polar  edge  of  the  auroral  zone  connect  to  the  boundary  layer  near  the 
magnetopause.  These  currents  are  thought  to  be  driven  by  the  solar  wind  MHD  dynamo  electric 
field  and  also  by  polarization  electric  fields  produced  by  the  boundary  layer  flow.  The  region  2  cur¬ 
rents  are  thought  to  connect  the  equatorward  side  of  the  auroral  electrojets  to  the  charged  Alfven 
layers  of  the  inner  magnetosphere  ring  current  region.  This  view  of  the  topology  of  the  region 
1  and  2  currents  is  common  and  is  shown  schematically  in  Figure  1.  Recently,  however,  studies 
suggest  that  the  dayside  and  nightside  sections  of  the  region  1  and  2  currents  may  be  separated 
and  may  be  produced  by  different  sources  (Friis-Christensen  et  al.,  1985;  Spence  et  al.,  1987).  On 
the  dayside,  a  variable  system  of  currents  flows  in  the  vicinity  of  the  polar  cleft.  The  cleft  system 
has  been  modeled  as  two  independent  currents  controlled  respectively  by  the  Y  and  Z  components 
of  the  interplanetary  magnetic  field.  These  two  systems  are  responsible  for  producing  the  well 


described  DPY  and  DPZ  ground  magnetic  perturbations  which  result  from  the  ionospheric  Hall 
currents  associated  with  these  systems  (D’Angelo,  1980;  Banks  et  al.,  1984;  Clauer  and  Banks, 
1986;  Clauer  and  Friis-Christensen,  1988).  These  cleft  currents  are  thought  to  be  driven  through  a 
direct  electrical  connection  to  the  solar  wind  MHD  dynamo  at  the  magnetopause. 

Apart  from  this  system  of  large  scale  currents,  systems  of  small  scale  currents  are  suspected 
to  be  associated  with  other  coupling  processes  which  are  sporadic,  localized,  and  in  some  cases, 
impulsive.  These  processes  include  the  so  called  Flux  Transfer  Events  associated  with  brief,  local¬ 
ized  reconnection  at  the  magnetopause  (Russell  and  Elphic,  1978;  1979;  Paschmann  et  al,  1982; 
Saunders  et  al.,  1984)  or  impulsive  plasma  penetration  in  which  bursts  of  plasma  penetrate  into  the 
magnetopause  boundary  layer  (Heikkila,  1979;  Lemaire,  1978).  The  relative  importance  of  such 
processes,  in  addition  to  Kelvin  Helmholtz  interactions,  is  currently  an  area  of  controversy  and 
vigorous  investigation.  For  example,  recent  theoretical  work  (Southwood,  1987;  Lee,  1988)  as  well 
as  some  experimental  evidence  (Cowley,  1982)  suggests  that  localized  reconnection  or  FTEs  at  the 
dayside  magnetopause  may  play  a  key  role  in  the  electrodynamic  coupling  between  the  solar  wind 
and  the  magnetosphere  and  the  ionosphere. 

The  magnetc«pheric  polar  cleft  is  a  dayside  topological  feature  of  the  magnetosphere  through 
which  magnetosheath  plasma  may  penetrate  to  low  altitudes  along  magnetic  field  lines  which  thread 
the  entry  and  boundary  layers  of  the  magnetopause.  This  is  illustrated  schematically  in  Figure  2. 
The  cleft  is  a  latitudinally  narrow  and  longitudinally  extended  region  of  solar  plasma  entry  to 
ionospheric  altitudes. 

Scientific  appreciation  of  the  unique  role  of  the  cleft  within  the  magnetosphere  -  ionosphere 
system  is  reflected  in  the  wide  variety  and  increasing  numbers  of  ground-based,  rocket,  balloon, 
and  satellite  investigations  which  are  each  contributing  new  information  about  the  various  different 
physical  processes  which  occur  within  this  relatively  small  portion  of  the  dayside  auroral  region 
and  polar  cap.  Within  the  overall  context  of  magnetospheric  research,  studies  of  the  cleft  and 
the  phenomena  associated  with  it  are  beginning  to  elucidate  the  full  complexity  of  the  interaction 
of  the  solar  wind  with  the  magnetosphere.  There  is  substantial  evidence  that  the  observations 
within  the  cleft  provide  the  most  direct  information  about  the  electrodynamics  of  the  solar  wind  - 
magnetosphere  -  ionosphere  interaction  since  current  systems  and  electric  fields  within  this  region 
appear  to  respond  directly,  with  little  time  delay,  to  variations  in  the  solar  wind  (Clauer  et  al., 
1984;  Clauer  and  Banks,  1986;  Clauer  and  Friis-Christensen,  1988). 

For  these  reasons,  our  research  efforts  have  focused  upon  the  coordinated  analysis  of  a  variety  of 
electrodynamic  observations  utilizing  correlative  satellite  and  ground  data  acquired  in  the  vicinity 
of  the  polar  cleft.  In  this  report,  we  summarize  the  results  of  a  three-year  analysis  and  modeling 
effort  supported  by  AFGL  Contract  F19628-85-K-0001.  The  acquisition  of  the  data  utilized  in 
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SOLAR  WIND-MAGNETOSPHERE  DYNAMO 
AND  DP2  EQUIVALENT  CURRENT 

Figure  1-  Schematic  representation  of  region  1  and  region  2  currents  driven  by  the  solar 
wind  imposed  emf  Em.  Alfven  layer  electric  field  shielding  gives  rise  to  the  potential  drop 
between  the  magnetopause  boundary  layer  and  the  ring  current  inner  edge.  The  ionospheric 
Hall  currents  which  result  are  similar  to  the  DP2  equivalent  current  pattern. 
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these  investigations  has  been  supported  by  other  sources.  The  AFGL  Contract  F19628-85-K-0001 
supported  the  coordinated  analysis  of  these  data.  The  objective  of  the  research  effort  has  been  to 
develop  an  electrodynamic  model  of  the  connectivity  of  the  dayside  high  latitude  ionosphere  with 
the  solar  wind  and  magnetosphere.  To  reach  these  objectives,  our  research  has  been  directed  along 
three  major  divisions:  1)  Observations  and  model  simulations  of  the  dayside  large  scale  high  latitude 
current  systems,  2)  Detailed  coordinated  analysis  of  combined  satellite  and  ground  based  data  to 
elucidate  current  strengths  and  configurations  under  particaular  solar  wind  conditions,  and  3)  The 
investigation  of  small  scale  coupling  phenomena  through  a  detailed  examination  of  phenomena 
occuring  in  the  vicinity  of  the  dayside  high  latitude  convection  reversal  boundary  which  we  feel  are 
the  ionospheric  manifestations  of  processes  occuring  at  the  magnetopause  and  within  the  entry  and 
boundary  layers.  Progress  within  these  three  areas  of  investigation  are  discussed  in  detail  in  the 
following  sections. 

LARGE  SCALE  HIGH  LATITUDE  CURRENT  SYSTEMS 

THEORETICAL  INTRODUCTION 

The  origin  of  the  polar  cap  electric  field  is  generally  described  in  terms  of  two  dominant  mecha¬ 
nisms  acting  at  the  outer  edges  of  the  magnetosphere.  The  first  of  these  is  magnetic  merging  at  the 
dayside  magnetopause,  a  process  which  depends  strongly  upon  the  interplanetary  magnetic  field 
(IMF)  orientation.  The  second  mechanism  is  related  to  a  viscous-like  interaction  at  the  magneto- 
spheric  boundary  layer  and  is  responsible  for  the  two-cell  convection  pattern  encompassing  the  night 
side  polar  cap  and  auroral  zone.  This  process  is  much  less  dependent  upon  the  IMF  orientation.  A 
number  of  studies  [Friis-Christensen  and  Wilhjelm,  1975;  Maezawa,  1976;  Friis-Christensen,  1979; 
Wilhjelm  et  al.,  1978;  Sergeev  and  Kuznetsov,  1981;  and  others]  demonstrate  the  separability  of 
the  two  mechanisms  and  provide  information  about  their  dependence  upon  the  IMF  and  solar  wind 
velocity  as  well  as  their  relationship  to  the  polar  cap  and  auroral  zone.  The  results  of  these  studies 
indicate  that  the  magnetic  field-aligned  currents  associated  with  merging  are  concentrated  near  the 
dayside  cleft  region.  Those  arising  from  the  boundary  mechanism  extend  towards  the  midnight 
sector  along  the  sides  of  the  auroral  oval  and  have  been  identified  as  the  region  1  currents. 

Banks  et  al.,  [1984]  suggest  a  model  in  which  the  pattern  of  dayside  high  latitude  field-aligned 
currents  and  ionospheric  plasma  convection  is  directly  determined  by  the  interplanetary  electric 
field  (IEF).  This  model  is  the  result  of  a  synthesis  of  a  variety  of  observations  and  the  recognition 
that,  due  to  the  geometrical  properties  of  the  magnetosphere,  the  east-west  component  of  the  IEF 
could  drive  field-aligned  currents  which  connect  to  the  ionosphere  at  points  lying  on  either  side  of 
local  noon,  while  currents  associated  with  the  north-south  component  of  the  IEF  will  connect  the 
two  polar  caps  as  sheet  currents  centered  at  local  noon.  Figure  3  from  Banks  et  al.,  [1984]  shows  a 
schematic  description  of  these  currents  and  their  relationship  to  the  solar  wind  velocity,  magnetic 
field,  and  motional  electric  field. 
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PANEL  A  PANEL  B 


Figure  3:  A  general  view  of  the  way  deft  currents  driven  by  the  interplanetary  electric  field 
connect  across  the  magnetosphere.  In  Figure  3 A,  line  currents  arising  from  Ey  connect 
from  the  flanks  of  the  magnetosphere  to  two  discrete  zones  on  either  side  of  local  noon. 
The  direction  of  the  current  is  determined  by  the  sign  of  B,  and  the  case  for  B,  >  0  is 
shown.  In  Figure  3B,  currents  arising  from  E,  connect  to  each  cleft  through  sheet  currents 
centered  at  local  noon.  Completion  of  the  current  circuit  through  the  opposite  hemisphere 
is  accomplished  by  means  of  a  sheet  current  along  closed  field  lines  lying  at  or  near  the 
magnetopause.  The  direction  of  currents  shown  here  corresponds  to  By  <  0.  Oppositely 
directed  current  would  result  for  By  >  0.  After  Banks  et  a/., [1984], 
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Thu  model  has  guided  our  research  effort  to  determine  the  configuration  of  large  scale  current 
systems  at  ^'gh  latitudes  on  the  dayside  of  the  magnetosphere  and  their  control  by  solar  wind 
parameters.  These  currents,  according  to  the  model,  are  driven  by  the  solar  wind  emf  imposed 
across  the  magnetopause.  The  model  follows  a  simple  decomposition  of  the  dayside  high  latitude 
field-aligned  currents  into  systems  which  are  separately  controlled  by  the  IMF  By  and  Bt  compo¬ 
nents  (recalling  that  apart  from  the  solar  wind  velocity,  these  are  equivalent  to  the  IEF  Ez  and  Ev 
components). 

These  field-aligned  current  systems  (shown  in  Figure  3)  are  referred  to  as  the  DPY  and  DPZ 
current  systems  as  they  are  associated  with  the  ionospheric  systems  which  have  originally  been 
labeled  with  these  names  [Friis-Christensen  and  Wilhjeim,  1975].  Here,  however,  we  modify  slightly 
the  definition  of  the  DPZ  disturbance.  In  the  original  definition,  the  DPZ  disturbance  included 
the  entire  two-cell  current  system  also  known  as  DP2.  We  now  limit  the  definition  of  the  DPZ 
system  to  the  more  local  two-cell  pattern  observed  on  the  dayside  at  times  when  |£*|  >>  |By|. 
This  system  corresponds  to  the  cusp  currents  statistically  described  by  Iijima  and  Potemra,  [1976] 
for  IMF  Bx  southward  and  to  the  system  observed  by  Araki  et  al.,  [1984]  and  termed  NPZ  currents 
by  I(jima  et  al.,  [1984'  Zanetti  et  al.,  [1984],  and  I(jima  and  Shibaji,  [1987].  The  DPY  system  is  a 
single  ionospheric  system  produced  by  Hall  currents  flowing  between  a  pair  of  oppositely  directed 
field-aligned  current  sheets  oriented  longitudinally  across  and  centered  near  the  noon  meridian,  as 
observed  by  Wilhjeim  et  al.,  [1978]. 

OBSERVATIONS,  ANALYSIS,  AND  MODEL  SIMULATIONS 

Our  investigation  of  the  large  scale  electrodynamic  current  coupling  between  the  dayside  high 
latitude  ionosphere  and  the  solar  wind  has  taken  two  approaches.  Both  rely  upon  high  latitude 
electrodynamic  measurements  from  the  Sondre  Stromfjord,  Greenland  incoherent  scatter  radar  as 
the  primary  data  set. 

The  first  approach  uses  a  numerical  computer  model  to  examine  the  ionospheric  electric  fields 
produced  by  different  field-aligned  current  distributions  in  the  cleft  [Rich  and  Kamide,  1983;  Banks 
et  al.,  1984].  The  model  can  then,  in  turn,  be  compared  with  measures  of  ionospheric  electric  fields 
obtained  by  high  latitude  radars.  In  our  research,  we  have  adapted  the  computer  model  used  by 
Banks  et  al.,  [1984]  to  adjust  the  strength,  direction,  and  location  of  the  DPY  and  DPZ  field- 
aligned  currents  in  the  polar  cleft  based  upon  IMF  measurements  obtained  by  the  IMP-8  satellite 
upstream  in  the  solar  wind.  Then,  through  a  series  of  steady  state  calculations,  we  simulate 
the  temporally  varying  ionospheric  convection  pattern  produced  by  the  model  during  a  particular 
interval  of  time.  The  simulated  convection  is  then  directly  compared  with  the  actual  convection 
measured  by  the  Sondre  Stromfjord,  Greenland  radar.  Results  obtained  from  the  first  approach 
have  been  presented  in  two  papers,  Clauer  and  Banks  [1986]  and  Clauer  and  Friis-Christensen 


[1988]  which  are  attached  in  Appendix  1.  A  brief  summary  of  the  results  presented  in  these  papers 
follows. 

The  initial  investigations  presented  by  Clauer  and  Banks  [1986]  utilized  a  distribution  of  currents 
which  are  static  in  location  but  changed  in  strength  and  direction  according  to  measurements  of 
the  IMF  By  and  Bz  components.  Figure  4  shows  a  polar  projection  of  these  field-aligned  current 
distributions  in  the  northern  hemisphere.  The  current  distribution  includes  the  region  1  and  region 
2  currents  using  the  approximate  statistical  distribution  and  spatial  arrangement  from  Iijima  and 
Potemra  [1976]  and  the  suggested  DPY  and  DPZ  field-aligned  currents  in  the  polar  cleft  from 
Banks  et  al.,  [1984].  This  figure  shows  the  current  flow  directions  for  the  case  when  Bz  >  0  and 
By  >  0.  When  Bz  =  By  =  0  no  cleft  currents  would  flow  and  for  Bz  <  0  and  By  <  0  the  polarities  of 
the  cleft  currents  would  reverse.  The  region  1  and  2  currents  are  constant  in  strength  and  location. 

For  each  measurement  time  at  the  Sondre  Stromfjord  radar,  we  take  the  measured  IMF  at  IMP- 
8  using  the  appropriate  delay  time  determined  from  the  actual  solar  wind  conditions  during  the 
event.  Based  upon  these  IMF  values,  we  adjust  the  strength  and  direction  of  the  DPY  and  DPZ 
currents  and  compute  the  hemispheric  electric  field  distribution.  Using  the  computed  electric  field, 
the  E  x  B  drift  at  the  radar  measurement  points  are  determined  and  the  simulated  convection 
velocity  pattern  is  plotted. 

Plates  3,  4,  and  5  in  Clauer  and  Banks  [1986]  attached  in  Appendix  1  show  measured  and 
simulated  high  latitude  ionospheric  convection  on  three  days  which  were  investigated.  The  corre¬ 
sponding  IMF  measurements  obtained  by  IMP-8  are  shown  in  Figures  5,  6,  and  7  of  that  same 
paper.  The  first  of  these  days,  April  25,  1983,  shows  a  changes  in  sign  for  the  IMF  By  component 
from  negative  to  positive  to  negative  which  results  in  changes  in  the  direction  of  the  high  latitude 
ion  drift  from  eastward  to  westward  and  back  to  eastward.  The  model  reproduces  this  change  in 
drift  associated  with  the  change  in  sign  of  By. 

The  second  interval,  March  16,  1984,  is  particularly  interesting  because  of  an  extended  region  of 
poleward  flow  from  1300  UT  through  14:45  UT.  This  interval  of  strong  poleward  flow  is  associated 
with  a  large  southward  IMF.  The  simulation  reproduces  the  observed  poleward  flow  at  high  latitudes 
suggesting  that  strong  DPZ  currents  similar  to  those  in  the  model  exist  during  strong  southward 
IMF  conditions.  At  lower  latitudes,  around  75  degrees  invariant,  the  simulation  disagrees  with  the 
observed  convection  strength  and  direction.  This  suggests  a  different  configuration  of  the  region 
1  and  2  currents  in  this  region  under  these  circumstances.  In  addition,  the  interval  following  the 
strong  southward  IMF  is  an  interval  of  strong  northward  IMF  and  this,  also,  is  not  modeled  well. 

The  third  interval,  July  24,  1983,  shows  a  variety  of  complex  temporal  variations.  Firom  12:30 
UT  through  13:30  UT  a  region  of  strong  poleward  flow  associated  with  strong  southward  IMF 
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FIELD  ALIGNED  CURRENT 


Figure  4:  Polar  plot  of  the  distribution  of  field-aligned  current  used  in  the  computer  model. 
Three  regions  of  current  are  shown:  (1)  the  region  1  and  region  2  currents  around  the  auroral 
oval  with  inward  current  denoted  by  shading,  (2)  the  DPY  sheet  currents  centered  at  local 
noon  and  lying  on  and  poleward  of  the  region  l  currents,  and  (3)  the  DPZ  line  currents  at 
the  longitudinal  ends  of  the  DPY  currents.  The  polarities  of  the  currents  in  this  illustration 
correspond  to  the  case  of  IMF  By  >  0  and  B,  >  0. 
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is  observed  and  is  similar  to  the  previous  example.  However,  examination  of  the  timing  of  the 
southward  variation  and  the  onset  of  poleward  flow  reveals  a  problem. 

Clauer  et  al.t  [1984]  and  Clauer  and  Banks,  [1986]  show  that  the  dayside  high  latitude  ionospheric 
electric  field  generally  responds  to  polarity  changes  in  the  IMF  By  and  Bt  components  within  about 
10  to  15  minutes.  On  July  24,  1983,  however,  the  poleward  flow  observed  at  12:30  UT  is  associated 
with  a  southward  turning  of  the  IMF  Bz  component  which  occurs  at  11:15  UT.  The  expected 
delay  time  between  IMP-8  and  the  magnetopause  is  about  8  to  12  minutes.  Thus,  the  time  delay 
from  when  the  southward  fluctuation  is  expected  to  encounter  the  magnetopause  and  the  onset  of 
observed  poleward  flow  in  the  ionosphere  is  over  50  minutes. 

We  propose  two  possible  explanations  for  this  observation.  Examination  of  the  IMF  observations 
indicate  that  By  was  large  and  positive  at  the  time  Bz  turned  southward.  By  became  nearly  0  at 
about  13:00  UT.  This  suggests  the  possibility  that  the  relative  magnitudes  of  By  and  Bz  may  be 
important  in  determining  the  ionospheric  effects  of  the  DPY  and  DPZ  currents.  Alternatively, 
perhaps  the  effects  of  the  DPZ  current  were  not  observed  until  around  12:30  UT  because  of  a  shift 
in  the  spatial  distribution  of  the  DPZ  currents.  It  is  interesting,  however,  that  the  simulation  of 
ion  velocity  measurements  agrees  very  well  with  the  observations.  The  interval  of  poleward  flow 
between  1230  and  1330  UT  is  reproduced.  Since  the  locations  of  the  currents  are  constant  in  the 
simulation,  this  suggests  that  adjusting  the  relative  strengths  of  the  DPY  and  DPZ  currents  is 
sufficient  to  explain  the  observations.  We  also  note  that  following  1330  UT  the  simulated  and 
observed  convection  velocities  do  not  agree.  This  is  during  an  interval  of  strongly  northward  IMF. 
While  the  results  of  the  Clauer  and  Banks,  [1984]  paper  show  that  a  remarkably  simple  model  can 
account  for  much  of  the  observed  high  latitude  dayside  convection  observations,  further  refinement 
is  required  to  better  address  the  situation  for  IMF  Bz  northward  and  for  the  auroral  electrojet 
region  near  local  magnetic  noon.  Therefore,  to  address  more  carefully  the  situation  for  both  IMF 
Bz  northward  and  southward,  Clauer  and  Friis-Christensen,  [1988]  have  refined  the  field- aligned 
current  distributions.  The  new  distributions  of  field-aligned  current  proposed  by  Clauer  and  Friis- 
Christensen  are  shown  in  Figure  5  and  their  corresponding  ionospheric  electric  potential  patterns 
are  shown  in  Figure  6. 

The  modified  distributions  of  current  better  reflect  our  understanding  of  the  dayside  high  latitude 
currents  and  our  ideas  regarding  the  origins  of  these  current  systems.  We  hypothesize  that  the  region 
1  currents  map  to  the  vicinity  of  the  inner  edge  of  the  low  latitude  boundary  layer  and  the  region  2 
currents  map  to  the  ring  current  region  of  the  inner  magnetosphere.  Since  the  field-aligned  current 
is  observed  to  be  continuous  in  the  observations  of  low  altitude  polar  satellites,  this  implies  that  the 
current  is  continuous  throughout  the  region  of  the  magnetosphere  between  the  ring  current  and  the 
low  latitude  boundary  layer.  During  periods  of  northward  IMF,  the  region  1  and  region  2  currents 
diminish  substantially  in  strength  and  are  located  at  higher  latitudes.  Near  local  noon,  where 
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Ionospheric  Distribution  of  Field-Aligned  Current 
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Figure  5:  Polar  plot  of  the  distribution  of  field-aligned  current  used  in  the  revised  computer 
model.  The  shaded  region  with  solid  contours  indicates  current  into  the  ionosphere  and 
the  unshaded  region  with  dashed  contours  indicates  current  out  of  the  ionosphere  in  the 
northern  hemisphere.  The  contour  interval  is  a  0.5  microamps  per  square  meter. 
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Ionospheric  Electric  Potential  Pattern 


By<° 


By  =  0 


By  >  « 


POTENTIAL  BYNBZP 


POTENTIAL  BYOBZP 


POTENTIAL  BYPBZP 


POTENTIAL  BYNBZO 


POTENTIAL  BYOBZO 


POTENTIAL  BYPBZO 


/  //,/ 


IWI 


POTENTIAL  BYNBZN 


POTENTIAL  BYOBZN 


POTENTIAL  BYPBZN 


w! 


» 


Figure  6:  Polar  plot  of  the  ionospheric  electric  potential  distribution  associated  with  each 
of  the  current  distributions  shown  in  Figure  5.  The  contour  interval  is  8  kilovolts.  The  zero 
potential  curve  crosses  the  center  of  the  polar  cap  and  solid  contours  indicate  positive  values 
and  dashed  contours  indicate  negative  values. 


the  low  latitude  boundary  layer  is  quite  thin  and  the  boundary  layer  flow  is  small,  the  primary 
contribution  to  the  region  1  current  system  are  the  DPY  and  DPZ  currents  which  result  from  a 
direct  electrical  connection  between  the  solar  wind  and  the  magnetosphere.  We  do  not  address  how 
this  electrical  connection  is  accomplished  but  feel  that  it  is  consistent  with  the  variety  of  magnetic 
merging  models  discussed  in  the  literature. 

The  DPY  and  DPZ  field-aligned  currents  are  each  separately  controlled  by  the  IMF  By  and  Bz 
components  (or  equivalently  by  the  interplanetary  electric  field  Ez  and  Ev  components  as  discussed 
by  Banks  et  al.,  [1984]  and  D’Angelo,  [1980]).  The  DPY  ground  magnetic  disturbance  is  produced 
by  the  ionospheric  Hall  current  flowing  between  a  pair  of  oppositely  directed  field-aligned  current 
sheets  oriented  longitudinally  across  and  centered  near  the  noon  meridian.  One  of  the  current 
sheets  is  collocated  with  the  region  1  currents  and  the  other  current  sheet  is  located  poleward  of 
the  region  1  currents.  During  northward  IMF  the  DPY  currents  are  located  at  higher  latitudes  but 
maintain  their  relationship  relative  to  the  region  1  currents.  The  poleward  sheet  lies  on  field  lines 
near  the  poleward  boundary  of  the  polar  cleft  and  the  equatorward  sheet  which  partially  forms  the 
noon  sector  region  1  currents  lies  near  the  equatorward  boundary  of  the  cleft.  Associated  with  the 
IMF  By  component  is  an  Ez  electric  field  which  drives  the  DPY  currents. 

For  southward  IMF  (negative  Bz)  the  DPZ  current  system  is  a  field-aligned  current  wedge  with 
current  flowing  from  the  vicinity  of  the  prenoon  equatorial  magnetopause  into  the  ionosphere  at  the 
latitude  of  the  region  I  current  and  out  of  the  ionosphere  in  the  postnoon  region  1  current  along 
field  lines  which  extend  to  the  vicinity  of  the  postnoon  equatorial  magnetopause.  During  northward 
IMF  the  DPZ  system  (or  NBZ  system  as  discussed  by  some  researchers)  is  a  similar  current  wedge 
with  opposite  current  direction  but  located  on  field  lines  near  the  poleward  boundary  of  the  polar 
cleft  and  thus  collocated  with  the  poleward  sheet  of  the  DPY  current  system.  Associated  with  the 
IMF  Bz  component  is  and  Ey  electric  field  which  drives  the  DPZ  currents  at  the  longitudinal  ends 
of  the  cleft. 

The  latitudinal  changes  in  the  currents  agree  better  with  observations  and  are  more  consistent 
with  the  theoretical  notions  regarding  merging  on  lobe  field  lines  poleward  of  the  cleft  for  northward 
IMF  and  merging  on  field  lines  equatorward  of  the  cleft  for  southward  IMF  [Dungey,  1963,  Crooker, 
1979;  Reiff,  1982;  Luhmann  et  al.,  1984;  Chiu  et  al.,  1985,  and  references  therein]. 

Clauer  and  Friis-Christensen,  [1988]  examine  the  consequences  of  the  new  distributions  of  current 
and  compare  with  observations  on  July  23,  1983.  During  this  day,  the  IMF  turns  strongly  northward 
for  several  hours.  The  observations  on  this  day  represent  an  extreme  case  and  therefore  are  very 
interesting  for  investigation. 

Perhaps  one  of  the  most  noteworthy  observations  presented  in  the  Clauer  and  Friis-Christensen 
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paper  is  the  rapid,  almost  immediate,  response  of  the  dayside  high  latitude  ionospheric  currents  to 
the  northward  turning  of  the  IMF.  The  IMF  (Figure  7)  was  unusually  large  on  July  23,  1982,  having 
total  field  values  of  about  22  nT  or  larger  for  much  of  the  day.  At  10:15  UT  the  field  measured  at 
IMP-8  turns  strongly  northward  with  Bz  becoming  about  20  nT.  It  remains  northward  for  several 
hours.  The  Bv  component  of  the  field  is  much  smaller  than  Bz  from  10:18  until  about  13:10  UT. 
The  magnitude  of  By  becomes  larger  than  Bz  at  13:30  UT.  The  time  delay  between  IMP-8  and  the 
magnetopause  is  estimated  to  be  8  minutes  indicating  that  the  northward  Bz  polarity  change  would 
encounter  the  magnetopause  at  10:23  UT.  The  change  in  the  high  latitude  currents  is  observed  to 
occur  at  the  Greenland  high  latitude  magnetic  observatories  Thule  (THL),  Kuvdlorssuaq  (KUV) 
and  Upernarvik  (UPN).  Station  locations  are  shown  in  Table  1  in  Clauer  and  Friis-Christensen, 
[1988]  in  Appendix  1.  The  D  component  of  the  magnetic  field  measurements  are  shown  in  Figure  8. 
The  response  of  the  currents  is  seen  to  begin  3  minutes  following  the  northward  IMF  encountering 
the  magnetopause  and  to  develop  over  a  period  of  about  22  minutes. 

The  observed  ionospheric  convection  pattern  is  shown  in  Figure  9.  To  create  this  figure  we  have 
used  the  radar  measured  ionospheric  convection  data  plus  the  magnetic  data  from  the  three  mag¬ 
netic  observatories  THL,  KUV,  and  UPN.  The  magnetic  data  are  15  minute  averages  and  the  total 
horizontal  field  component  is  rotated  90  degrees  counter  clockwise  to  give  the  equivalent  convection 
direction.  The  magnetic  measurements  are  not  calibrated  with  the  radar  velocity  measurements  so 
the  magnitude  of  the  vectors  is  only  roughly  proportional  to  the  velocity. 

From  0800  to  1300  UT,  the  convection  reversal  boundary  appears  to  be  located  between  73 
and  75  degrees  invariant  latitude.  The  anti-sunward  convection  observed  at  latitudes  above  this 
boundary  is  thought  to  be  tailward  convection  on  polar  cap  field  lines.  The  low  latitude  of  the 
boundary  indicates  an  expanded  polar  cap  and  is  consistent  with  the  several  intervals  of  large 
southward  IMF  prior  to  1015  UT.  Extending  the  convection  measurements  further  poleward  using 
magnetic  data,  we  see  a  sunward  rotation  of  the  convection  following  10:15  UT.  The  sunward 
convection  is  restricted  to  very  high  latitudes. 

Figure  10  shows  the  model  simulation  based  upon  the  current  distributions  shown  in  Figure  9. 
Figure  10  can  be  directly  compared  with  Figure  9.  The  convection  reversal  boundary  is  located  at 
about  the  same  location  in  both  figures.  The  simulated  high  latitude  vectors  at  the  magnetometer 
measurement  points  at  dawn  and  dusk  are  aligned  more  east-  west  than  the  measured  vectors, 
however,  this  is  likely  due  to  skewing  of  the  magnetometer  vectors  due  to  conductivity  gradients 
and  the  associated  field-aligned  currents.  The  highest  latitude  vectors  measured  by  THL  are  well 
reproduced  by  the  simulation.  The  rotation  is  seen  at  1030  UT  in  both  sets.  The  measurements  at 
KUV,  the  next  highest  latitude  vectors,  are  not  as  well  simulated  at  1030  UT  but  are  fairly  well 
modeled  after  1100  UT.  Thus,  the  spatial  arrangement  of  the  high  latitude  currents  only  approxi¬ 
mates  the  actual  distribution,  but  the  agreement  is  quite  reasonable  and  the  time  variations  match 
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Figure  7:  Interplanetary  magnetic  field  measured  by  the  IMP-8  spacecraft  during  July  23, 
1983  shown  in  GSM  coordinates. 
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Figure  8:  Horizontal  magnetic  field  declination  variations  measured  relative  to  the  average 
declination  at  Thule  (top),  Kuvdlorssuaq  (middle)  and  Upernavik  (bottom).  At  10:23  UT 
the  northward  IMF  is  estimated  to  encounter  the  magnetopause  and  at  10:26  UT  the  field 
at  THL  and  KUV  begins  a  120° rotation. 
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MEASURED  ION  VELOCITY 


Figure  9:  Polar  plot  showing  F-region  ion  convection  velocity  vectors  as  a  function  of  latitude 
and  local  time  for  July  23,  1983.  Local  time  is  indicated  around  the  outer  circumference  of 
the  plot  and  universal  time  is  indicated  around  the  inner  circumference  of  the  plot.  Local 
magnetic  noon  is  at  the  top  of  the  plot.  The  three  complete  rows  of  highest  latitude  vectors 
are  convection  velocities  inferred  from  ground  magnetometer  data.  The  other  vectors  in  the 
plot  are  obtained  using  the  Sondre  Stromfjord,  Greenland  incoherent  scatter  radar. 
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Figure  10:  Simulated  F-region  ion  convection  velocity  vectors  for  July  23,  1983. 


well.  Between  12  and  13  UT  the  magnitudes  of  the  vectors  equatorward  of  the  convection  reversal 
boundary  are  smaller  than  observed.  This  has  been  a  consistent  feature  in  most  of  the  simula¬ 
tions  shown  here  and  those  that  have  not  been  published.  We  intend  to  examine  the  conductivity 
measurements  in  this  region  to  determine  if  that  is  the  source  of  the  disagreement. 

A  more  thorough  discussion  of  these  data  and  the  results  is  available  in  the  paper  by  Clauer  and 
Friis-Christensen  attached  in  Appendix  1. 

COORDINATED  ANALYSIS  OF  SATELLITE 
AND  GROUND  BASED  MEASUREMENTS 

The  second  approach  to  our  investigation  of  the  large  scale  dayside  high  latitude  current  systems 
uses  the  coordinated  analysis  of  simultaneous  measurements  of  ionospheric  convection  obtained  by 
the  Sondre  Stromftord  radar  and  in  situ  observations  of  field-aligned  currents,  particle  precipitation 
and  plasma  drift  obtained  by  low  altitude  polar  satellites  nearly  conjugate  with  the  radar  measure¬ 
ments.  This  work  has  been  discussed  in  our  quarterly  progress  reports,  has  been  presented  in  papers 
at  some  professional  meetings,  and  is  still  in  progress.  Papers  are  presently  in  preparation  and  we 
will  discuss  our  results  as  of  the  completion  of  our  contract.  Results  from  the  analysis  of  these 
coordinated  simultaneous  measurements  should  aid  us  in  setting  parameters  more  accurately  in  the 
model  current  systems  discussed  above.  A  variety  of  coordinated  data  are  available  for  analysis  of 
dayside  high  latitude  electrodynamics  during  March  7,  16  and  20,  1984.  We  have  conducted  an 
number  of  studies  of  the  coordinated  electrodynamic  parameters  measured  by  the  Sondre  Strom- 
(jord,  Greenland  incoherent  scatter  radar  together  with  in  situ  measurements  from  low  altitude 
polar  satellites.  Measurements  in  the  pre-noon  sector  have  utilized  the  DMSP  F7  satellite  and 
measurements  in  the  post-noon  sector  have  utilized  the  HILAT  satellite.  During  these  days,  we 
have  a  variety  of  different  plasma  convection  patterns  corresponding  to  IMF  Bz  >  0  and  <  0  and 
IMF  By  >  0  and  <  0. 

MARCH  7,  1984  '  Pre  Noon  hours. 

A  DMSP  satellite  pass  occurred  at  1130  UT  which  was  nearly  conjugate  with  Sondrestrom. 
DMSP  magnetometer  data  show  2  sheets  of  current  with  the  equatorward  sheet  being  outward  and 
the  poleward  sheet  directed  inward.  Figure  11  shows  the  DMSP  magnetometer  data  together  with 
the  derived  current  density.  Sketched  along  the  top  of  the  plot  is  a  schematic  indication  of  the  type 
of  particle  precipitation  observed  by  the  plasma  instrument.  In  Figure  12  we  show  the  ionospheric 
convection  measured  by  the  Sondrestrom  radar  at  the  time  of  the  DMSP  pass.  We  have  also 
indicated  schematically  the  locations  of  the  DMSP  observed  currents  and  particle  precipitation. 
The  region  marked  “oval”  represents  the  region  of  hard  auroral  precipitation  consisting  of  energies 
of  about  1  KeV.  The  “cusp”  region  consists  of  high  fluxes  of  low  energy  (50  eV)  particles,  and  the 
“polar  cap”  region  consists  of  very  small  fluxes.  We  note  that  the  cusp  plasma  is  observed  in  a 
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region  of  sunward  flowing  convection  and  the  electric  field  is  rather  strong  (long  vectors). 

Figures  13  and  14  show  data  obtained  from  an  elevation  scan  made  by  the  Sondrestrom  radar. 
Figure  13a  shows  the  electron  density  and  Figure  13b  shows  the  ion  temperature.  In  Figure  13a 
there  is  a  density  enhancement  located  at  about  100  km  north  of  the  radar  at  150  km  altitude.  This 
E-region  enhancement  is  due  to  the  1  KeV  oval  precipitation.  Poleward  of  this,  in  the  region  of  the 
large  fluxes  of  low  energy  cusp  precipitation,  the  ionospheric  density  drops  to  very  low  values.  The 
cusp  precipitation  is  not  producing  am  ionospheric  density  enhancement  in  the  F-region  as  might 
be  expected.  Figure  13b  provides  some  insight  as  to  why.  Very  high  ion  temperatures  are  measured 
in  the  cusp  precipitation  region.  These  high  temperatures  are  the  result  of  joule  heating  due  to 
the  high  electric  fields  and  large  convection  velocities.  The  consequence  of  the  high  temperatures 
is  that  the  recombination  rate  is  increased  and  the  ionization  is  eliminated. 

These  observations  occur  during  an  interval  of  IMF  By  <  0  and  the  IMF  observations  measured 
by  IMP-8  are  shown  in  Figure  14.  A  summary  of  the  observations  superposed  upon  a  schematic 
of  the  global  convection  pattern  assumed  for  IMF  By  <  0  is  shown  in  Figure  15.  The  features 
which  we  show  on  the  figure  include,  the  cusp  plasma  observed  primarily  in  the  poleward  region 
of  the  upward  (region  2)  but  spreading  slightly  into  the  inward  region  1  current.  The  hard  oval 
precipitation  lies  in  the  equatorward  portion  of  the  region  2  current.  The  region  1  current  lies  in 
the  region  of  primarily  polar  cap  precipitation. 

The  lack  of  F-region  ionization  observed  in  the  region  of  cusp  precipitation  is  a  curious  ob¬ 
servation  since  ionosonde  observations  have  often  been  utilized  to  identify  the  polar  cusp.  If  this 
observation  it  a  general  feature,  then  those  ionosonde  observations  which  have  claimed  to  identify 
the  cusp  could  really  be  oblique  reflections  from  the  strong  E-region  enhancement  at  lower  latitudes. 
This  is  something  which  should  be  investigated  more  thoroughly  through  the  analysis  of  additional 
data. 

MARCH  16,  1984'  P™  Noon  hours. 

On  March  16,  1984,  the  DMSP  pass  at  1T.50  UT  is  nearly  conjugate  with  the  observations 
of  the  Sondrestrom  radar.  The  DMSP  magne'ometer  on  this  pass  observed  3  sheets  of  current; 
the  equatorward  and  poleward  sheets  being  outward  and  the  middle  sheet  being  inward.  This  is 
shown  in  Figure  16  using  the  same  format  as  Figure  11.  A  summary  of  the  particle  observations 
is  schematically  shown  along  the  top  of  Figure  16.  In  contrast  with  the  observations  on  March  7, 
1984,  the  region  2  current  is  entirely  associated  with  1  KeV  oval  precipitation  and  the  soft  cusp 
precipitation  occurs  across  the  entire  region  1  current  and  into  the  poleward  cusp  current  sheet.  The 
boundary  between  sunward  convection  and  anti-sunward  convection  lies  at  the  boundary  between 
the  region  1  and  region  2  currents  and  is  shown  in  the  Sondrestrom  convection  measurements  made 
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Figure  11:  Magnetic  field  measurements  and  derived  field-aligned  current  density  (shaded 
curve)  obtained  by  the  DMSP  polar  satellite  during  a  northern  hemisphere  pass  beginning 
at  11:30  UT  on  March  7,  1984.  UT,  Geomagnetic  Latitude  and  Magnetic  Local  Time 
are  shown  along  the  horizontal  axis,  magnetic  field  strength  on  the  left  axis,  and  current 
density  in  pA/m2  on  the  right  axis.  Along  the  top  of  the  plot  is  a  schematic  representation 
of  the  particle  precipitation  measured  by  the  plasma  instruments  on  the  satellite.  The  oval 
represents  fluxes  of  KEV  energy  particles,  the  cusp  is  a  region  of  enhanced  fluxes  of  low 
(few  Ev)  energy  particles,  and  the  polar  cap  is  a  region  of  low  precipitation  fluxes.  Current 
directed  into  the  ionosphere  is  positive,  outward  current  is  negative. 
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Figure  12:  Measurements  of  ionosphpric  plasma  convection  obtained  by  the  Sondre  Strom- 
fjord,  Greenland  radar  during  the  period  of  the  March  7,  1984,  11:30  UT  DMSP  satellite 
pass.  The  location  of  the  DMSP  track  is  indicated  on  the  plot  along  with  a  schematic 
representation  of  the  field-aligned  current  direction  and  plasma  precipitation  regions. 
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Figure  13  a)  Measurements  of  the  ionospheric  electron  density  along  the  Sondrestrom  mag¬ 
netic  meridian  obtained  during  the  interval  11:32  -  11:37  UT  on  March  7,  1984.  A  density 
enhancement  of  2  •  105  electrons/cc  is  observed  at  E-region  altitudes  about  100  km  north  of 
the  radar,  b)  ion  temperature  measurements  obtained  during  the  same  meridian  scan.  A  F- 
region  temperature  enhancement  is  observed  poleward  of  the  density  enhancement  measured 
in  panel  a). 
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Figure  14:  Interplanetary  magnetic  field  measurements  on  07  March  1984. 
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at  the  time  of  the  DMSP  satellite  pass  and  shown  in  Figure  17. 


These  observations  occur  during  an  interval  of  IMF  By  >  0.  The  IMF  measurements  obtained 
by  IMP-8  are  shown  in  Figure  18.  A  summary  of  the  observations  superposed  upon  a  schematic  of 
the  global  convection  pattern  assumed  for  IMF  By  >  0  is  shown  in  Figure  19. 

MARCH  20,  1984:  Pn  Noon  hours. 

The  prenoon  observations  on  March  20,  1984  are  very  similar  to  those  for  March  16,  1984.  The 
IMF  observations  are  shown  in  Figure  6.  Several  DMSP  passes  have  been  examined,  and  all  show 
current  structures  similar  to  those  observed  on  March  16,  1984.  DMSP  passes  occur  at  10:28  UT, 
12:08  UT  and  13:50  UT.  The  12:08  pass  is  measuring  nearly  the  same  meridian  as  the  Sondrestrom 
radar.  The  DMSP  data  are  shown  in  Figure  20  for  the  three  passes.  In  Figure  21  we  show  as 
schematic  representation  of  the  DMSP  observations  together  with  the  plasma  convection  measured 
by  the  Sondrestrom  radar.  The  consistency  of  the  field-aligned  currents  with  the  longitudinally 
separated  measurements  of  the  convection  at  10:28  UT  and  13:30  UT  suggest  that  the  triple  current 
sheet  structure  observed  in  the  pre-noon  hours  during  IMF  By  >  0  extends  longitudinally  through 
a  region  of  several  hours  of  local  time. 

In  contrast  with  the  observations  on  March  7,  1984,  the  region  2  current  is  entirely  associated 
with  1  KeV  oval  precipitation  and  the  soft  cusp  precipitation  occurs  across  the  entire  region  1 
current  and  into  the  poleward  cusp  current  sheet.  The  boundary  between  sunward  convection  and 
anti-sunward  convection  lies  at  the  boundary  between  the  region  1  and  region  2  currents  and  is 
shown  in  the  Sondrestrom  convection  measurements  made  at  the  time  of  the  DMSP  satellite  pass 
and  shown  in  Figure  17. 

These  observations  occur  during  an  interval  of  IMF  By  >  0.  The  IMF  measurements  obtained 
by  IMP-8  are  shown  in  Figure  18.  A  summary  of  the  observations  superposed  upon  a  schematic  of 
the  global  convection  pattern  assumed  for  IMF  By  >  0  is  shown  in  Figure  19. 

MARCH  7,  1984:  Post  Noon  Hours 

Some  preliminary  analysis  has  been  completed  of  coordinated  HILAT  magnetometer  data  and 
Sondrestrom  radar  data  during  the  post  noon  hours.  These  obserbvations  occured  during  intervals 
of  variable  IMF  Bz  and  Bz  «  0  while  By  <  0  for  both  passes.  The  IMF  observations  are  shown  in 
Figure  14  and  the  HILAT  magnetic  observations  are  shown  in  Figure  22.  The  HILAT  data  show 
considerable  structure  in  the  field-aligned  currents  with  multiple  current  sheets  observed.  The 
top  panel  of  Figure  22  shows  the  HILAT  pass  beginning  at  15:24  UT.  There  are  5  current  sheets 
observed.  The  observations  for  the  17:04  pass  (bottom  panel)  show  considerably  more  structure, 
however,  there  are  basicly  3  major  current  sheets.  Figure  23  shows  the  Sondrestrom  convection 


—  1 


25 


it 


Observation  Summary  for 
March  7, 1 984 


Oval 


DMSP  F 7  orbit  track 
.  1 1 :30  UT 


Oust 


1  KeV  precipitation 
producing  ionization 
\  at  150  Km. 


r 

Polar  cap  y 

Vi 

/  /  ' 

^85 

?  so  / 

75 

Figure  15:  Schematic  summary  of  March  7.  1984,  11:30  UT  DMSP  observations  superim¬ 
posed  upon  the  large  scale  convection  pattern  assumed  for  IMF  Dv  <  0. 
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Figure  17:  Measurements  of  ionospheric  plasma  convection  obtained  by  the  Sondre  Stronv 
fjord,  Greenland  radar  during  the  period  of  the  March  16,  1984,  11:50  UT  DMSP  satellite 
pass.  The  location  of  the  DMSP  track  is  indicated  on  the  plot  along  with  a  schematic 
representation  of  the  field-aligned  current  direction  and  plasma  precipitation  regions. 
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Figure  18:  Interplanetary  magnetic  field  measurements  on  16  March  1984. 
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Observation  Summary  for 
March  16, 1984 


Figure  19:  Schematic  summary  of  March  16,  1981.  11  50  UT  DMSP  observations  superim¬ 
posed  upon  the  large  scale  convection  pattern  assumed  for  IMF  Dv  >  0. 
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Figure  21:  Ionospheric  convection  measured  by  the  Sondrestrom  radar  during  each  of  the 
DMSP  satellite  passes,  10:30  UT  (a),  12:13  UT  (b)  and  13:54  UT  (c),  shown  in  Figure  20. 
Also  shown  is  the  track  of  the  satellite  pass  and  a  schematic  representation  of  the  current 
regions. 
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Figure  22:  Magnetic  field  and  computed  field-aligned  current  density  for  two  consecutive 
HILAT  satellite  passes  on  March  7,  1984  at  15:24  UT  (top)  and  17.04  UT  (bottom).  These 
passes  are  in  the  13  -  14  local  time  sector.  Current  directed  inward  is  positive,  outward 
directed  current  is  negative. 
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Figure  23:  Ionospheric  convection  measured  by  the  Sondrestrom  radar  during  each  of  the 
HILAT  satellite  passes,  15:24  UT  (top)  and  17:04  UT  (bottom)  shown  in  Figure  22.  Also 
shown  is  the  track  of  the  satellite  pass  and  a  schematic  representation  of  the  observed  current 
regions. 
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Figure  24:  Magnetic  field  and  computed  field-aligned  current  density  for  two  consecutive 
HILAT  satellite  passes  on  March  16,  1984  at  14:52  UT  (top)  and  16:35  UT  (bottom).  These 
passes  are  in  the  12  -  14  local  time  sector.  Current  directed  into  the  ionosphere  is  positive, 
outward  current  is  negative 
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Figure  25:  Ionospheric  convection  measured  by  the  Sondrestrom  radar  during  each  of  the 
HILAT  satellite  passes,  14:52  UT  (top)  and  16:35  UT  (bottom)  shown  in  Figure  24.  Also 
shown  is  the  track  of  the  satellite  pass  and  a  schematic  representation  of  the  observed  current 
regions. 


observations  for  each  of  the  HILAT  satellite  passes. 

MARCH  16,  1984:  Post  Noon  Hours 

Preliminary  analysis  of  the  March  16,  1984  HILAT  magnetic  data  and  correlative  Sondrestrom 
radar  data  have  also  been  completed.  The  IMF  observations  on  March  16  are  shown  in  Figure  18. 
The  HILAT  magnetic  data  are  shown  in  Figure  24.  The  HILAT  pass  beginning  at  14:50  (top  panel) 
occurs  during  and  interval  of  IMF  By  <  0  and  Bz  <  0.  The  HILAT  pass  beginning  at  16:30  UT 
(bottom  panel)  occurs  during  a  data  gap  in  the  IMP-8  data,  however,  it  seems  possible  to  guess 
that  during  the  gap  Bz  >  0  and  By  >  0.  The  14:52  UT  pass  shows  5  current  sheets  while  the  16:35 
UT  pass  shows  considerable  structure  and  multiple  current  sheets,  though  possibly  the  sheets  form 
a  basic  3  sheet  pattern 

SUMMARY 

These  sets  of  coordinated  measurements  provide  3  variety  of  conditions  where  explicit  knowledge 
of  the  field-aligned  current  structure  and  ionospheric  convection  are  known.  In  continuing  investi¬ 
gations  supported  through  other  grants,  we  will  incorporate  these  current  configurations  into  our 
modeling  effort  to  better  understand  the  complex  current  distributions  in  the  dayside  cleft  region, 
their  relationship  to  high  latitude  ionospheric  convection  and  their  control  by  solar  wind  parame¬ 
ters.  Figure  25  shows  the  Sondrestrom  ionospheric  convection  measurements  obtained  during  the 
HILAT  pass  intervals  shown  in  Figure  24. 

SMALL  SCALE  COUPLING  PHENOMENA 

Our  investigations  of  small  scale  or  localized  solar  wind-magnetosphere-ionosphere  coupling  has 
concentrated  in  three  complimentary  areas  involving  small  scale  field-aligned  currents  in  the  vicin¬ 
ity  of  the  polar  cleft.  (1)  Theoretical  modeling  of  the  ground  magnetic  signatures  and  ionospheric 
consequences  of  hypothetical  flux  transfer  event  current  systems;  (2)  An  alysis,  in  collaboration  with 
Dr.  Eigil  Friis-Christensen,  of  magnetometer  data  from  the  Greenland  chain  and  particularly  the 
high  time  resolution  data  from  the  array  of  temporary  stations  operated  around  Sondrestrom  during 
the  summer  of  1986;  and  (3)  analysis  of  data  from  radar  experiments  which  focus  on  rapid  mea¬ 
surements  of  the  ionospheric  convection  around  the  convection  reversal  boundary.  Progress  in  these 
areas  has  been  very  good  with  a  paper  published  describing  our  modeling  techniques  [McHenry  and 
Clauer,  1987]  and  the  publication  of  our  discovery  this  summer  of  ionospheric  traveling  convection 
vortices  (TCYs)  in  the  Greenland  magnetometer  data  [Friis-Christensen  et  al . ,  1988].  Both  of  these 
papers  are  attached  to  Appendix  1. 

Our  rapid  measurements  of  the  ionospheric  convection  near  the  reversal  boundary  have  revealed 
that  the  boundary  can  be  either  steady  or  very  dynamic.  In  the  dynamic  condition,  it  can  move 
several  hundred  kilometers  in  a  few  minutes.  Our  initial  efforts  to  understand  the  cause  of  these 
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Figure  26:  Data  from  the  Greenland  magnetometer  chain  showing  the  II  perturbation  in  the 
left  panel  and  the  Z  perturbation  in  the  right  panel.  All  of  the  stations  are  along  the  west 
coast  except  for  the  bottom  4  stations  in  each  panel  which  are  located  on  the  east  coast. 
The  large  single  cycle  pulsation  at  10:10  UT  is  the  TCV  event  under  consideration. 
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several  hundred  kilometer  movements  of  the  boundary  were  inconclusive.  However,  we  are  now 
investigating  the  association  of  TCVs  with  the  observed  deflections  of  the  boundary  and  there 
appears  to  be  a  relationship.  The  combined  radar  and  ground  magnetometer  analysis  promises  to 
reveal  new  information  about  the  dynamics  of  the  convection  reversal  boundary. 

Theoretical  modeling  efforts  have  concentrated  on  the  ground  magnetic  signature  to  be  expected 
from  small  scale  field-aligned  current  systems  as  hypothesized  for  the  so  called  flux  transfer  events 
(FTEs).  FTEs  are  the  basis  for  a  popular  theory  to  explain  the  observed  transfer  of  magnetic 
flux  from  closed  to  open  field  lines  through  impulsive,  localized  magnetic  reconnection  [Russell  and 
Elphic,  1979].  FTEs  have  been  observed  by  satellites  close  to  the  magnetopause  using  a  repeatable 
in  situ  magnetic  perturbation  caused  by  strong  field-aligned  currents.  This  magnetic  signature  is 
associated  also  with  observations  of  plasma  inside  the  flux  tube  having  a  mixture  of  magnetospheric 
and  magnetosheath  sources. 

One  of  our  goals  has  been  to  detect  the  signature  of  these  events  on  the  ionosphere.  Because  the 
plasma  connecting  the  ionosphere  and  the  area  where  FTEs  have  been  detected  obeys  the  “frozen 
in”  MHD  approximation,  current  densities  and  plasma  velocities  in  the  magnetosphere  near  the 
magnetopause  and  in  the  ionosphere  should  be  very  similar  except  for  a  scale  factor  depending 
upon  the  background  magnetic  field  strength  which  is  fairly  well  known. 

Our  theoretical  effort  has  been  to  use  the  hypothesized  field-aligned  current  structure  of  FTEs, 
assume  that  it  is  the  same  at  the  ionosphere  except  for  the  magnetic  scale  factor  and  solve  for  the  I 

electric  field  in  the  ionosphere.  Because  the  dayside  ionosphere  is  almost  uniform,  it  is  straight 
forward  to  solve  for  the  ionospheric  currents  and  then  the  magnetic  field  on  the  ground.  This  is  pre¬ 
sented  in  the  paper  by  McHenry  and  Clauer  [1987].  A  copy  of  this  paper  showing  the  mathematical 
details  and  the  resulting  magnetic  field  strengths  and  patterns  found  for  the  FTE  current  systems  | 

is  attached  in  Appendix  1.  Simple  estimates  for  the  rate  at  which  ground  magnetometers  would 
detect  FTEs  are  also  computed.  Using  FTE  detection  statistics  at  the  magnetopause,  we  estimate 
the  detection  rate  of  FTEs  at  a  ground  station  under  the  dayside  convection  reversal  boundary  to  be 
between  12  and  60  minutes  per  sighting  depending  upon  which  model  current  system  is  considered. 

We  have  also  examined  the  ground  induction  effects  caused  by  moving  current  systems.  Results 
similar  to  the  case  for  time  varying  currents  fixed  in  location  were  found.  At  high  ionospheric 
current  speeds,  the  vertical  magnetic  field  on  the  ground  will  diminsh  and  the  horizontal  field  will 
double.  For  translation  speeds  of  convection  vortices  in  the  ionosphere  of  10  km/sec  and  greater,  I 

induction  effects  are  very  strong.  For  speeds  much  below  10  km/sec,  induction  effects  are  weak. 

Thus,  for  the  observed  speeds  which  are  between  1  and  4  km/sec,  induction  effects  will  be  small 
and  we  neglect  them.  The  large  change  in  ground  conductivity  known  as  the  ’coast  effect’  was  also 
examined.  Because  the  induction  currents  penetrate  to  depths  of  several  tens  of  kilometers,  i.e. 
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much  more  then  the  depth  of  the  ocean,  the  large  jump  in  surface  conductivity  between  the  sea 
and  land  should  not  effect  geomagnetic  observations  of  ionospheric  currents  moving  less  then  10 
km/sec.  This  implies  that  vortex  like  disturbances  in  the  ionosphere,  which  create  a  radial  ground 
magnetic  field,  will  still  generate  a  radial  pattern  on  the  ground  even  near  the  sea  coast  where  most 
of  the  Greenland  high  latitude  observatories  are  located. 

Guided  by  our  modeling  efforts  and  in  collaboration  with  Dr.  Eigil  Friis-Christensen,  we  have 
focused  our  data  analysis  on  a  set  of  pulsation  bursts  observed  in  the  Greenland  magnetometer 
data.  First  results  from  this  investigation  showing  the  discovery  of  traveling  convection  vortices 
moving  tailward  near  the  ionospheric  convection  reversal  boundary  have  been  published  by  Friis- 
Christensen  et  al.,  [1988]  attached  in  Appendix  1.  Figure  26  shows  such  a  pulsation  recorded  by  the 
Greenland  magnetometer  chain  with  1-minute  data  samples.  During  the  summer  of  1986,  an  array 
of  magnetometer  stations  was  operated  around  Sondre  Stromfjord,  Greenland  collecting  data  at  20- 
second  resolution.  Figure  27  shows  a  series  of  plots  from  this  array  showing  the  horizontal  magnetic 
perturbation  vectors  plotted  at  the  geographic  location  of  each  station.  The  vectors  appear  to  be 
pointing  radially  from  a  source  which  is  at  a  different  location  in  each  frame,  ie.  moving  across 
the  array  of  stations.  The  source  of  this  radial  magnetic  magnetic  perturbation  is  thought  to  be 
a  moving  field-aligned  current  filament.  These  six  frames  of  data  represent  only  2-minutes  of  the 
20-second  data.  In  Figure  28  we  plot  the  entire  event  in  a  new  format  in  a  single  plot. 

If  we  assume  that  the  magnetic  perturbations  which  are  observed  by  the  stations  result  from 
a  steady  current  system  moving  across  the  array  of  stations,  we  can  equally  well  represent  this 
as  moving  the  array  of  stations  across  a  fixed,  steady  current  system.  In  Figure  28  we  plot  the 
horizontal  perturbation  vectors  (which  have  been  rotated  by  -90°to  approximate  the  direction 
of  the  overhead  convection)  from  the  array  and  offset  each  set  of  measurements  by  a  distance 
which  corresponds  to  the  velocity  of  the  moving  current  structure.  In  this  example,  the  motion 
of  the  current  is  4  km/sec.  This  velocity  can  be  judged  fairly  precisely  since  the  array  of  stations 
contains  one  site  located  east  of  the  meridional  chain.  The  result  is  a  twin  vortex  pattern  which 
moves  westward  (tailward)  over  the  stations.  (The  stations  are  located  at  about  10  Local  Time). 
Numerous  events  or  traveling  convection  vortices  (TCVs)  such  as  this  have  been  identified  in  the 
data  and  are  presently  under  investigation.  At  least  two  classes  of  TCVs  have  been  identified. 
One  class  is  associated  with  large  single  cycle  pulsations  such  as  shown  in  this  example.  Another 
class  consists  of  large  intervals  (hours)  of  irregular  pulsations  which  can  be  explained  as  series  of 
oppositely  directed  field-aligned  currents  moving  steadily  across  the  array  of  stations. 

Least  squares  fitting  the  modeled  ground  magnetic  field  of  a  single  field-aligned  current  filament 
to  Greenland  magnetometer  data,  traveling  vortices  in  the  dayside  ionospheric  flow  have  been 
detected  and  tracked.  They  are  observed  during  quiet  periods  for  several  hours  at  a  time  and  move 
generally  anti-sunward.  Assuming  each  vortex  to  be  the  convection  pattern  produced  by  a  small 
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Figure  27:  Magnetic  perturbations  are  plotted  in  successive  panels  at  the  location  of  the 
corresponding  station  for  the  temporary  array  of  observatories  operated  around  Sondre 
Stromfjord,  Greenland.  This  data  is  obtained  at  20-second  resolution  so  the  time  difference 
between  each  panel  is  20  seconds.  The  total  horizontal  components  of  the  magnetic  pertur¬ 
bation  are  shown  with  a  solid  vector  while  the  vertical  component  perturbations  are  shown 
with  a  vertical  dashed  line.  At  10:08  UT,  the  disturbance  approached  the  stations  from  the 
east  causing  a  horizontal  perturbation  which  points  radially  from  a  source  position  as  well 
as  a  large  vertical  (dashed  lines)  perturbation.  At  10:10  UT  the  disturbance  is  centered  over 
the  array  of  stations. 
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Figure  28:  The  magnetic  perturbations  in  the  vertical  component  (top)  and  horizontal 
component  (bottom)  from  the  array  of  20-second  stations  for  the  interval  10:06  UT  to  10:21 
UT.  The  horizontal  component  has  been  rotated  90  degrees  counter  clockwise  to  show  the 
equivalent  ionospheric  convection  direction.  We  ofTset  the  station  locations  by  80  km  for 
each  set  of  20  second  data,  corresponding  to  a  velocity  of  about  4  km/sec.  The  resulting  twin 
vortex  pattern  is  indicative  of  a  pair  of  oppositely  directed  field-aligned  current  filaments 
moving  rapidly  over  the  array  of  stations. 
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field  aligned  current  moving  across  the  ionosphere,  the  amount  of  field  aligned  current  is  found  to 
be  10  to  50  kamps  and  they  move  at  1-3  km/sec  across  the  ionosphere.  During  some  periods,  the 
vortices  are  aligned  into  ’trains’  and  the  calculated  field  aligned  current  is  seen  to  be  steady  for  each 
vortex  and  neighboring  vortices  have  currents  of  opposite  sign.  Other  intervals  have  vortices  which 
again  move  generally  anti-sunward,  have  alternating  field  aligned  current  but  the  vortex  trajectories 
and  amount  of  field  aligned  current  are  irregular.  These  different  situations  indicate  that  ‘here  are 
different  classes  of  small  scale  convection  near  the  convection  reversal.  Simultaneous  measurements 
using  the  Sondre  Stromflord  Radar  show  the  vortices  are  traveling  along  the  convection  reversal 
boundary.  Because  the  vortices  are  conjugate  to  the  magnetopause  shear  region  and  repeat  in  a 
regular  fashion,  we  suggest  that  they  probably  results  from  the  Kelvin-Helmholtz  instability  at  the 
magnetopause.  This  is  supported  by  a  correlation  between  field  aligned  current  strength  and  solar 
wind  speed.  Little  correlation  is  found  with  IMF  measurements. 

Some  time  periods  have  been  found  where  significant  pulsations  of  the  ground  magnetic  field 
exist  but  the  data  can  not  be  fit  with  the  modeled  field  of  the  simple  vortex.  Examining  the 
cotemporal  IMF  in  detail  of  these  cases,  we  have  found  that  the  IMF  has  a  high  variance.  It 
appears  that  these  observed  pulsations  are  directly  related  to  the  variance  in  the  IMF.  Most  likely, 
the  pulsations  are  caused  by  the  large  scale  high  latitude  field  aligned  currents  being  modulated  by 
the  driving  IMF.  Since  the  magnetic  field  of  the  large  scale  currents  is  quite  large,  modulating  them 
causes  irregular  pulsations  that  are  larger  than  the  ground  magnetic  field  caused  by  any  traveling 
vortices. 
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On  July  23,  1983,  the  Interplanetary  Magnetic  Field  turned  strongly  northward,  becoming 
about  22  nT  for  several  hours.  Using  a  combined  data  set  of  ionospheric  convection  measurements 
made  by  the  Sondre  Stromfjord  incoherent  scatter  radar  and  convection  inferred  from  Greenland 
magnetometer  measurements,  we  observe  the  onset  of  the  reconfiguration  of  the  high-latitude 
ionospheric  currents  to  occur  about  3  min  following  the  northward  IMF  encountering  the  magne¬ 
topause  The  large-scale  reconfiguration  of  currents,  however,  appears  to  evolve  over  a  period  of 
about  22  min.  Using  a  computer  model  in  which  the  distribution  of  field-aligned  current  in  the 
polar  cleft  is  directly  determined  by  the  strength  and  orientation  of  the  interplanetary  electric 
field,  we  are  able  to  simulate  the  ime-varying  pattern  of  ionospheric  convection,  including  the 
onset  of  high-latitude  "reversed  convection"  cells  observed  to  form  during  the  interval  of  strong 
northward  IMF.  These  observations  and  the  simulation  results  indicate  that  the  dayside  polar 
cap  electric  field  observed  during  strong  northward  IMF  is  produced  by  a  direct  electrical  current 
coupling  with  the  solar  wind. 


1.  Introduction 

It  is  well  established  that  magretosphcric  convection  and 
the  associated  ionospheric  plasma  convection  are  strongly 
controlled  by  the  strength  and  orientation  of  the  interplan¬ 
etary  magnetic  field  (IMF)  observed  at  the  magnetopause 
[ Heelis ,  1984;  Zanetti  et  al.,  1984:  Clauer  et  al.,  1984;  Re  iff 
and  Burch ,  1985;  Frits- Christensen  et  al-,  1985;  Clauer  and 
Banks,  1986;  and  references  therein).  A  simple  model  to 
explain  this  observed  control  has  been  suggested  by  Bunks 
et  al.  [1984].  In  this  model,  the  high-latitude  ionospheric 
electric  field  is  a  consequence  of  direct  electrical  coupling  be¬ 
tween  the  interplanetary  electric  field  and  the  high  latitude 
ionosphere  via  field-aligned  currents. 

This  model  was  adapted  by  Clauer  and  Banks  [1986[  to 
perform  a  series  of  electrostatic  computations  of  the  iono¬ 
spheric  electric  field  based  upon  measured  IMF  values  ob¬ 
tained  by  the  IMP  8  spacecraft.  Using  these  computa¬ 
tions,  Clauer  and  Banks  simulated  the  ion  convection  ve¬ 
locity  which  would  be  observed  by  the  Sondre  Stromfjord, 
Greenland,  incoherent  scatter  radar  under  the  conditions 
described  by  the  model.  The  results  of  the  simulations  com¬ 
pared  favorably,  in  general,  with  the  ionospheric  plasma  con¬ 
vection  measured  by  the  radar. 

The  time  intervals  which  were  studied  by  Clauer  and 
Banks  [1986)  consisted  primarily  of  intervals  having  south¬ 
ward  IMF  and  occasional  brief  intervals  of  northward  IMF. 
The  simulation  of  the  northward  IMF  intervals  did  not  agree 
well  with  the  observations.  Therefore,  in  this  investigation 
we  address  more  carefully  the  situation  for  both  IMF  Bz 
northward  and  southward.  The  current  distributions  used 
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by  Clauer  and  Banks  [1986)  have  been  refined  and  are  now 
more  consistent  with  the  field-aligned  current  systems  which 
we  hypothesize  to  connect  the  solar  wind  emf  with  the  day- 
side  high-latitude  ionosphere. 

The  model  follows  a  simple  decomposition  of  the  dayside 
high-latitude  field-aligned  rurrents  into  systems  which  are 
separately  controlled  by  the  IMF  Bv  and  Bz  components. 
These  field-aligned  current  systems  are  rc^rred  to  as  the 
"field-aligned”  DRY  and  DPZ  current  sysu..  s  as  they  are 
associated  with  the  ionospheric  systems  which  have  been 
originally  labeled  with  these  names  [Fms-Chnstensen  and 
Wtlhjelm,  1975).  Here,  however,  we  modify  slightly  the  defi¬ 
nition  of  the  DPZ  disturbance.  In  the  original  definition,  the 
DPZ  disturbance  included  the  entire  two-cell  current  system 
also  known  as  the  DPT.  system.  We  now  limit  the  definition 
of  the  DPZ  system  to  the  more  local  two-cell  pattern  ob¬ 
served  on  the  dayside  at  times  when  IMF  |Bj|  >>  \BV\. 
This  system  corresponds  to  the  cusp  currents  statistically 
described  by  /tjtma  and  Potemra  [1976]  for  IMF  Bz  south¬ 
ward  and  to  the  system  observed  by  Araki  et  al.  [1984] 
and  termed  NBZ  currents  by  Itjtma  et  al.  [1984],  In  our 
model,  the  field-aligned  DPZ  system  is  contained  entirely 
on  the  dayside  of  the  polar  cap  and  does  not  extend  into  the 
nightside  as  described  by  Iijtma  et  al.  [1984],  Zanetti  et  al. 
(1984),  and  Itjtma  and  Shtbajt  [1987].  The  DPY  system  is  a 
single  ionospheric  system  consisting  of  a  Hall  current  flowing 
between  a  pair  of  oppositely  directed  field-aligned  current 
sheets  oriented  longitudinally  across  and  centered  near  the 
noon  meridian,  as  observed  by  Wtlhjelm  et  al.  [1978). 

Refinements  to  the  field-aligned  DPY  and  DPZ  current 
systems  in  our  simulation  now  include  latitudinal  changes 
in  the  current  systems  when  IMF  Bz  changes  sign.  The 
latitudinal  changes  in  the  currents  agree  better  with  obser¬ 
vations  and  are  more  consistent  with  the  theoretical  notions 
regarding  merging  on  lobe  field  lines  poleward  of  the  cleft  for 
northward  IMF  and  merging  on  field  lines  equatorward  of 
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Fig.  1.  Interplanetary  magnetic  field  measurements  in  GSM  coordinates  during  July  23,  1983 


the  cleft  for  southward  IMF  [ Dvngey .  1963:  Crooker.  1979; 
Re  iff.  1982;  Luhmann  et  ai,  1984:  Chiu  et  ai.  1985  and 
references  therein]. 

To  examine  the  consequences  of  the  new  distributions  of 
currents  and  to  compare  with  observations,  we  are  fortu¬ 
nate  to  have  both  radar  convection  observations  and  IMF 
measurements  on  July  23,  1983,  during  an  interval  when  the 
IMF  turns  strongly  northward  for  several  hours.  This  day 
represents  an  extreme  case  and  therefore  is  very  interesting 
for  investigation. 

2.  Results 

Observations 

Figure  1  shows  measurements  of  the  interplanetary  mag¬ 
netic  field  in  GSM  coordinates  obtained  by  the  IMP  8  space¬ 
craft  located  at  Xgie  =  28  Re  and  V'98p  =  - 10  /?£.  The 
solar  wind  velocity  is  about  400  km/s  and  the  proton  num¬ 
ber  density  is  about  10/cm2.  Assuming  that  the  solar  wind 
dynamic  pressure  is  balanced  by  the  magnetic  pressure  of 
the  Earth’s  magnetic  field,  the  distance  from  the  center  of 
the  Earth  to  the  subsolar  magnetopause  is  about  7.5  Re 
Using  this  information,  we  estimate  that  the  transit  time 
between  IMP  8  and  the  magnetopause  is  about  7  to  9  min 
(see  Clauer  and  Banks  [1986]  for  a  description  of  the  esti¬ 
mation  procedure). 

The  IMF  was  unusually  large  on  July  23,  1983.  having 
total  field  values  of  about  22  nT  or  larger  for  much  of  the 
day.  At  1015  UT  the  field  turns  strongly  northward  with  Bs 
becoming  about  20  nT.  It  remains  strongly  northward  for 
many  hours.  The  By  component  of  the  field  is  much  smaller 
than  the  strong  B2  component  from  1018  until  about  1310 
UT.  The  magnitude  of  By  becomes  larger  than  Bt  at  1330 
UT. 

On  July  23.  1983.  the  Sondre  Stromfjord  radar  made  iono¬ 


spheric  convection  measurements  in  the  invariant  latitude 
range  between  70°  and  80°  with  about  25-min  time  resolu¬ 
tion.  In  addition  to  the  operations  of  the  radar,  the  Dan¬ 
ish  Meteorological  Institute  Greenland  magnetometer  chain 
operates  continuously  and  is  roughly  aligned  along  the  Son¬ 
dre  Stromfjord  magnetic  meridian.  In  Figure  2  the  convec¬ 
tion  observations  obtained  from  both  the  Sondre  Stromfjord 
radar  and  from  Greenland  magnetometer  measurements  are 
shown  in  a  clock  dial  format.  Local  noon  is  at  the  top  of 
the  plot,  and  local  time  is  shown  around  the  outer  circum¬ 
ference  of  the  plot.  UT  is  shown  around  the  interior  of  the 
plot,  and  the  radial  dimension  indicates  invariant  latitude. 
Velocity  vectors  are  shown  with  a  point  at  the  measurement 
position  and  a  line  in  the  direction  of  the  flow.  The  length 
of  the  line  is  proportional  to  the  flow  speed. 

In  general,  the  velocity  vectors  determined  from  the  Son¬ 
dre  Stromfjord  radar  data  lie  below  79.5°.  There  are  occa¬ 
sional  vectors,  however,  above  80.0°  obtained  when  signifi¬ 
cant  signal  is  available  in  the  most  distant  range  gates.  De¬ 
tails  of  the  operating  mode  and  data  reduction  are  described 
by  Clauer  and  Banks  [1986],  In  addition,  we  show  three  sets 
of  velocity  vectors  at  invariant  latitudes  79.9°,  82.6°.  and 
85.5°  which  were  inferred  from  the  magnetic  measurements 
of  the  three  northernmost  stations  in  the  west  Greenland 
magnetometer  chain.  The  station  names  and  coordinates 
are  given  in  Table  1. 

These  equivalent  velocity  vectors  were  obtained  by  ro¬ 
tating  the  horizontal  magnetic  perturbation  vector  by  90° 
to  give  the  equivalent  current  direction,  then  rotating  an¬ 
other  180°  to  give  the  convection  velocity  direction.  The 
magnetic  measurements  are  not  calibrated  with  the  radar 
velocity  measurements,  so  the  magnitude  of  the  vectors  is 
only  roughly  proportional  to  the  velocity.  The  direction  in¬ 
dicated  bv  the  magnetirally  derived  equivalent  convection 
velocity  vectors,  however,  is  generally  consistent  with  the 
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Table  1.  Greenland  Magnetic  Observatories 


Station 

Symbol 

Geographic 

North 

Geographic 

East 

Invariant 

Latitude 

Declination 

Thule 

THL 

77.48° 

290.83° 

85.5° 

S 

CO 

o 

Kuvdlorssuaq 

KUV 

74.57° 

302.82° 

81.6° 

-57° 

Upernavik 

UPN 

72.78° 

303.85° 

79.9° 

-54° 

actual  convection  direction.  Discrepancies,  nevertheless,  do 
occur  in  regions  where  there  are  conductivity  gradients,  and 
thus  the  equivalent  velocity  vectors  near  dawn  and  dusk  may 
be  rotated  at  some  angle  to  the  true  convection.  For  exam¬ 
ple,  this  may  explain  the  rotation  of  the  magnetic  vectors 
relative  to  the  highest-latitude  radar  observations  near  dusk. 
During  the  northern  hemisphere  summer,  the  dayside  high- 
latitude  ionosphere  has  a  fairly  uniform  conductivity  except 
near  the  dayside  terminator.  The  equivalent  convection  vec¬ 
tors  should  nevertheless  respond  in  a  consistent  manner  to 
changes  in  convection  produced  by  changes  in  the  external 
drivers  of  the  convection. 

Confining  our  attention  to  the  radar  data  which  lie  below 
79°,  apart  from  a  small  variation  in  the  latitudinal  loca¬ 
tion  of  the  convection  reversal  boundary  between  1000  UT 
and  1030  UT  we  see  no  change  in  convection  in  response  to 
the  northward  turning  of  IMF  Bt  at  1015  UT.  From  0800 
to  1300  UT,  the  convection  reversal  boundary  is  located  be¬ 
tween  71°  and  75°.  The  antisunward  convection  observed  at 
latitudes  above  this  boundary  is  thought  to  be  tailward  con¬ 
vection  on  polar  cap  field  lines.  The  low  latitude  of  the  con¬ 
vection  reversal  boundary  indicates  an  expanded  polar  cap 
and  is  consistent  with  the  several  intervals  of  large  south¬ 
ward  IMF  prior  to  1015  UT.  The  convection  is  strong,  being 
between  2  and  3  km/s.  Note  that  the  IMF  By  component 
is  extremely  large  during  this  interval  as  well. 

It  is  interesting  that  there  is  no  apparent  change  in  the 


radar-measured  convection  at  around  1015  UT  when  the 
IMF  B i  turned  northward.  The  radar  measurements  extend 
well  into  the  polar  cap  antisunward  convection.  Extending 
the  convection  measurements  further  poleward  by  use  of  the 
magnetometer  data  does,  however,  show  a  change  in  convec¬ 
tion  at  very  high  latitudes.  Both  KUV  at  A  =  81.6°  and 
THL  at  A  =  85.5°  observe  a  sunward  rotation  of  the  con¬ 
vection  velocity.  The  sunward  convection  is  restricted  to 
very  high  latitudes.  It  is  also  noteworthy  that  the  change  in 
convection  observed  at  KUV  and  THL  follows  the  polarity 
change  in  Bz  very  directly. 

Allowing  an  8-min  transit  time  from  IMP  8  to  the  mag¬ 
netopause,  the  polarity  change  would  encounter  the  mag¬ 
netopause  at  1023  UT,  and  the  ionospheric  response  ob¬ 
served  in  Figure  2  is  seen  at  1045.  In  Figure  3  we  have 
plotted  the  SD  component  of  the  magnetic  field  at  TLH, 
KUV  and  UPN  using  1-min  resolution  magnetometer  data. 
The  quantity  6D  is  the  perturbation  field  declination  an¬ 
gle  measured  relative  to  the  average  magnetic  declination 
at  each  station.  Examination  of  the  1-min  magnetometer 
data  shown  in  Figure  3  indicates  that  the  field  at  THL  and 
KUV  begins  a  slow  rotation  at  1026  UT  which  continues 
through  1058  UT.  Similar  timing  is  seen  at  UPN.  The  plot 
shows  that  at  THL  and  KUV  the  magnetic  perturbation 
vector  rotates  from  6D  =  30°  to  6D  =  150°  between  1026 
and  1058  UT.  Thus,  changes  in  the  ionospheric  currents  as¬ 
sociated  with  the  northward  turning  of  Bz  begin  about  3 


MEASURED  ION  VELOCITY 


Fig.  2.  Polar  plot  showing  F  region  ion  convection  velocity  vectors  as  a  function  of  latitude  and  local  time  for 
July  23,  1983.  Local  time  is  indicated  around  the  outer  circumference  of  the  plot,  and  universal  time  is  indicated 
around  the  inner  circumference.  Local  magnetic  noon  ( 1400  UT)  is  at  the  top  of  the  plot  The  three  complete  rows 
of  highest-latitude  vectors  are  convection  velocities  inferred  from  ground  magnetometer  data.  The  other  vectors 
in  the  plot  were  obtained  using  the  Sondre  Stromfjord.  Greenland,  incoherent  scatter  radar 
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Greenland  Magnetometer  Data 
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Fig.  3.  Horizontal  magnetic  field  declination  variations  measured  relative  to  the  average  declination  at  Thiile  (top), 
Ktivdlorssuaq  (middle),  and  Upernavik  (bottom).  At  1023  I  T  the  northward  IMF  is  estimated  to  encounter  the 
magnetopause,  and  at  1026  U'T  the  field  at  THL  and  KUV  begins  a  120°  rotation 


min  following  the  estimated  time  that  the  northward  field 
encounters  the  magnetopause.  This  delay  time  is  compara¬ 
ble  with  the  estimated  propagation  time  for  an  Alfven  wave 
to  travel  from  the  dayside  magnetopause  to  the  ionosphere. 
We  also  notice  that  during  the  rotation  of  the  field,  there 
is  a  periodic  variation  in  the  field  rotation  with  a  period  of 
about  6  min.  Perhaps  this  is  due  to  multiple  reflections  of 
the  incident  Alfven  wave  which  initiates  the  current  recon¬ 
figuration.  The  large-scale  reconfiguration  of  the  currents 
appears  to  evolve  over  the  22  min  during  which  the  field 
rotation  at  THL  and  KUV  takes  place. 

Model  Simulation 

Using  the  technique  described  by  Claut  r  and  Banks  [1986] 
we  have  used  a  computer  model  to  simulate  the  convection 
pattern  observed  on  July  23,  1983.  The  model  assumes  that 
the  distribution  and  strength  of  field-aligned  currents  in  the 
dayside  high-latitude  regions  are  controlled  by  the  direction 
and  magnitude  of  the  IMF.  For  a  given  conductivity  dis¬ 
tribution  and  field-aligned  current  distribution,  the  model 
computes  the  global  ionospheric  electrostatic  potential  dis¬ 
tribution.  Here  we  have  adopted  the  realistic  but  simple 
conductivity  model  from  Kamide  and  Matsushita  [1979]  with 
the  additional  inclusion  of  a  conductivity  enhancement  in 
the  auroral  oval.  From  the  computed  electric  potential,  the 
E  x  B  F-region  plasma  drift  can  be  determined  Thus,  us¬ 
ing  the  measurements  of  the  IMF  recorded  by  IMP  8  to 
determine  the  distribution  and  magnitude  of  the  dayside 
high-latitude  field-aligned  currents,  we  compute  the  time- 
varying  ionospheric  electrostatic  field  and  the  associated 
plasma  drift  which  would  be  measured  by  the  Sondre  Stro m- 
fjord  radar  and  high-latitude  Greenland  magnetometers. 


We  have  now  modified  the  distribution  of  current  to  bet¬ 
ter  reflect  our  understanding  of  the  dayside  high-latitude 
currents  and  our  ideas  regarding  the  origins  of  these  current 
systems.  The  distributions  of  current  in  our  model  have 
been  guided  by  the  following  observational  and  theoretical 
considerations.  We  hypothesize  that  the  region  1  currents 
map  to  the  vicinity  of  the  low-latitude  boundary  layer,  and 
the  region  2  currents  map  to  the  ring  current  region  of  the  in¬ 
ner  magnetosphere.  Since  the  current  is  observed  to  be  con¬ 
tinuous  in  the  observations  of  low-altitude  polar  satellites, 
we  assume  that  the  current  is  continuous  throughout  the 
region  in  the  magnetosphere  between  the  ring  current  and 
the  low-latitude  boundary  layer.  During  periods  of  north¬ 
ward  IMF.  the  region  1  and  2  currents  diminish  substantially 
in  strength  and  are  located  at  higher  latitudes.  Near  local 
noon,  where  the  low-latitude  boundary  layer  is  thin  and  the 
boundary  layer  flow  is  small,  the  primary  contributions  to 
the  region  1  current  system  are  the  field-aligned  DPY  and 
Z?  PZcurrcnts  which  result  from  a  direct  electrical  connection 
between  the  solar  wind  and  the  magnetosphere.  We  do  not 
address  how  this  eleetriral  connection  is  accomplished  but 
feel  that  it  is  consistent  with  the  variety  of  magnetic  merg¬ 
ing  models  discussed  in  the  literature.  Also,  since  we  are 
primarily  concerned  with  dayside  observations,  these  con¬ 
siderations  may  oversimplify  the  situation  in  the  nightside 
magnetosphere. 

In  our  model,  we  consider  a  set  of  independent  dayside 
high-latitude  field-aligned  current  systems  each  separately 
controlled  by  the  IMF  By  ami  Bz  components  (or  equiva¬ 
lently  by  the  interplanetary  electric  field  Ez  and  Ey  com- 
po’n  nts  as  discussed  by  Banks  et  al.  [1984]  and  D'Angelo 
[19X11,;  In  the  model,  the  DPY  ground  magnetic  distur¬ 
bance  is  produced  by  the  ionospheric  Hall  current  flowing 
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Ionospheric  Distribution  of  Field- Aligned  Current 


By  <  0 


By  >  0 


Bz  >  0 


Bz  =  0 


Bz  <  0 


Current  into  the  ionosphere 


□  Current  out  of  the  ionosphere 

Fig  4.  Polar  plot  of  the  distribution  of  field-aligned  current  used  in  the  computer  model.  The  shaded  region 
with  solid  contours  indicates  current  into  the  ionosphere  and  the  unshaded  region  with  dashed  contours  indicates 
current  out  of  the  ionosphere  in  the  northern  hemisphere.  The-  contour  interval  is  at  0  5  p.4/m2. 


between  a  pair  of  oppositely  directed  field-aligned  current  currents  are  located  at  higher  latitudes  but  maintain  their 
sheets  oriented  longitudinally  across  and  centered  near  the  relationship  relative  to  the  region  1  currents.  The  mapping 
noon  meridian  as  observed  by  Wilhjelm  et  al  [l'J78j.  One  of  these  currents  to  the  magnetopause  is  that  suggested  by 
of  the  current  sheets  is  collocated  with  the  region  1  cur-  Banka  et  al.  [1984]  and  shown  schematically  in  their  Figure 
rents,  and  the  other  current  sheet  is  located  poleward  of  2  The  poleward  current  sheet  lies  on  field  lines  near  the 
the  region  1  currents.  During  northward  IMF,  the  DPY  poleward  boundary  of  the  polar  cleft,  and  the  rquatorward 
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Ionospheric  Electric  Potential  Pattern 
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Fig.  5  Polar  plot  of  the  ionospheric  potential  distribution  assoi  uite'l  with  each  of  the  <  urrent  distributions  shown 
in  Figure  4.  The  contour  interval  is  8  kV.  The  zero  potential  curve  t  rosses  the  renter  of  t Fir-  polar  cap;  solid 
contours  indicate  positive  values,  and  dashed  contours  indicate  negative  values 


sheet  which  partially  forms  the  noon  sector  region  i  system 
lies  near  the  equatorward  boundary  of  the  cleft .  Associated 
with  I  lie  IMF  By  component  is  an  Ez  interplanetary  electric 
field  which  drives  the  DPI' current  sheets. 

For  southward  IMF  (negative  Bz)  the  field-aligned  DPZ 
current  system  is  a  current  wedge  with  current  flowing  from 
the  vicinity  of  the  prenoon  equatorial  magnetopause  into  the 
ionosphere  at  the  latitude  of  the  region  1  current  and  out  of 
the  ionosphere  in  the  postnoon  region  1  current  along  held 
lines  which  extend  to  the  vicinity  of  the  postnoon  equato¬ 
rial  magnetopause  (also  shown  schematically  in  Bank >  */ 


al.  [1984]  Figure  2).  During  northward  IMF  the  field- 
aligned  DPZ  system  (or  NBZ  system  as  discussed  by  some 
researchers  [o.g..  hjimu  rt  al..  1984])  is  a  similar  current 
wedge  with  opposite  current  direction  polarity,  but  located 
on  held  lines  near  the  poleward  boundary  of  the  polar  cleft 
-rid  thus  collocated  with  the  poleward  sheet  of  the  field- 
aligned  DPY  current  system.  Associated  with  the  IMF  Bz 
component  is  an  Ey  interplanetary  electric  field  which  drives 
the  DPZ  currents  at  the  longitudinal  ends  of  the  cleft. 

The  nine  plots  in  Figure  4  show  the  variability  with  re¬ 
spect  to  the  IMF  of  the  current  distributions  used  in  our 
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Fig.  6.  Simulated  F  region  ion  convection  velocity  vectors  for  July  23,  1983,  shown  in  the  same  format  as  Figure 
2. 


present  model.  The  pattern  in  the  center  shows  the  baseline 
Bx  =  By  =  0  configuration.  The  top  row  shows  patterns 
for  IMF  Bx  >  8  nT,  and  the  bottom  row  shows  patterns 
for  IMF  B,  <  -8  nT.  The  right  column  shows  patterns 
for  IMF  By  >  8  nT,  and  the  left  column  shows  patterns  for 
IMF  By  <  -8  nT.  For  larger  values  of  the  field  components, 
the  currents  remain  saturated  at  these  levels.  The  currents 
vary  linearly  in  strength  for  IMF  values  between  these  sat¬ 
uration  levels,  and  the  positions  change  abruptly  when  the 
Bx  component  changes  sign. 

For  positive  Bx  the  total  inward  or  outward  region  1  cur¬ 
rent  is  about  0.9  x  10®  A,  and  the  total  inward  or  outward 
region  2  current  is  about  0.7  x  106  A  while  for  large  neg¬ 
ative  Bz  the  total  inward  or  outward  region  1  current  is 
about  5.5  x  106  A  and  the  total  inward  or  outward  region  2 
current  is  about  4.4  x  106  A.  The  total  inward  or  outward 
field-aligned  current  in  the  DPY  system  is  about  1.5  x  106 
A  for  Bx  >  0  and  about  2.  x  106  A  for  Bz  <  0.  The  larger 
total  DPY  current  for  negative  Bz  is  primarily  a  result  of 
the  larger  spatial  area  of  the  currents  spanning  the  same 
longitudinal  range  but  at  lower  latitudes.  The  DPY  cur¬ 
rent  density  is  the  same  for  positive  or  negative  Bz  with 
1.9  pA/m 2  in  the  poleward  sheet  and  2.4  pA/m2  in  the 
equatorward  sheet.  The  current  density  in  the  field-aligned 
DPZ  currents  is  1.7  pA/m2.  The  total  inward  or  outward 
DPZ  current  is  about  1.  x  10®  A  for  negative  Bz  and  about 
0.6  x  106  A  for  positive  Bz. 

In  Figure  5  we  show  the  corresponding  electrical  potential 
patterns  for  these  nine  current  distributions.  These  poten¬ 
tial  patterns  are  derived  using  the  modeling  code  described 
by  Banks  et  of.  [1984]. 

In  order  to  simulate  the  convection  pattern  for  July  23, 
1983,  we  perform  a  series  of  steady  state  calculations  using 
distributions  of  currents  determined  by  the  IMF  measured 
using  the  IMP  8  satellite  and  corrected  for  propagation  de¬ 
lay.  A  20-min  offset  is  used  between  the  measurements  at 
IMP  8  and  the  simulated  velocity  measurements.  This  off¬ 
set  is  composed  of  the  8-min  transit  time  from  IMP  8  to 
the  magnetopause  and  the  12-min  delay  for  the  large-scale 
currents  to  develop  which  was  determined  by  Clauer  et  al. 


[1984]  and  Ctauer  and  Banks  [1986].  From  the  computed 
electric  field  we  calculate  the  ionospheric  E  x  B  drift  veloc¬ 
ity  at  the  radar  observation  points  and  display  them  in  a 
clock  dial  plot.  The  simulated  convection  pattern  for  July 
23,  1983,  is  shown  in  Figure  6.  The  simulated  velocity  vec¬ 
tors  are  only  determined  over  the  interval  where  IMF  mea¬ 
surements  are  available  from  IMP  8.  As  is  apparent  from 
Figure  1,  there  is  a  data  gap  from  1447  to  1929  UT. 

Figure  6  can  be  directly  compared  with  the  measurements 
shown  in  Figure  2.  The  convection  reversal  boundary  is 
located  at  about  the  same  location  in  both  figures.  The 
simulated  high-latitude  vectors  at  the  magnetometer  mea¬ 
surement  points  at  dawn  and  dusk  are  aligned  more  east- 
west  than  the  measured  vectors;  however,  this  is  likely  due 
to  skewing  of  the  magnetometer  vectors  due  to  conductiv¬ 
ity  gradients  and  the  associated  field-aligned  currents.  The 
highest-latitude  vectors  measured  by  THL  are  well  repro¬ 
duced  by  the  simulation.  The  rotation  is  seen  at  1030  UT 
in  both  sets.  The  measurements  at  KUV,  the  next  highest 
latitude  vectors,  are  not  as  well  simulated  at  1030  UT  but 
are  fairly  well  modeled  after  1200  UT.  Thus  the  spatial  ar¬ 
rangement  of  the  high-latitude  currents  only  approximates 
the  actual  distribution,  but  the  agreement  is  quite  reason¬ 
able  and  the  time  variations  match  well.  Between  1400  and 
1500  UT  the  magnitudes  of  the  vectors  equatorward  of  the 
convection  reversal  are  smaller  than  observed.  This  has  been 
a  consistent  feature  in  the  simulations  shown  by  Claner  and 
Banks  [1986J  and  here  also.  In  future  work  we  will  examine 
the  conductivity  measurements  in  this  region  to  determine 
if  that  is  the  source  of  the  disagreement. 

The  success  of  the  model  in  simulating  the  high-latitude 
time-varying  convection  during  a  transition  from  strongly 
negative  IMF  Bz  to  strongly  positive  IMF  Bz  provides  addi¬ 
tional  support  for  the  suggestion  of  direct  electrical  coupling 
between  the  solar  wind  and  high- latitude  ionosphere  during 
both  IMF  orientations.  The  poleward  motion  of  the  DPY 
and  DPZ  currents  is  consistent  with  motion  of  the  coupling 
region  from  field  lines  below  the  cleft  to  lobe  field  lines  above 
the  cleft  during  periods  of  strongly  northward  Br.  In  this 
instance,  this  occurs  during  an  active  interval  having  an  ex- 
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panded  polar  cap.  The  result  is  the  apparent  development 
of  “reversed"  convection  cells  within  the  polar  cap. 

3.  Discussion  and  Conclusions 

Perhaps  the  most  noteworthy  observation  presented  with¬ 
in  this  paper  is  the  rapid,  almost  immediate,  response  of  the 
dayside  high-latitude  ionospheric  currents  to  the  northward 
IMF.  We  are  fortunate  that  the  location  of  the  IMP  8  space¬ 
craft  yields  little  difference  in  travel  time  for  orientations  of 
the  IMF  discontinuity  which  lie  between  the  Parker  spiral 
angle  and  the  plane  orthogonal  to  the  Earth-Sun  line.  The 
most  likely  orientation  of  the  IMF  discontinuity  is  assumed 
to  lie  within  this  range  of  angles. 

Observing  an  ionospheric  response  within  a  few  minutes 
of  the  time  that  the  northward  IMF  is  expected  to  encounter 
the  magnetopause  suggests  very  strongly  that  a  direct  elec¬ 
trical  coupling  between  the  solar  wind  and  the  dayside  high- 
latitude  ionosphere  occurs  during  northward  IMF.  The  rapid 
response  is  also  somewhat  in  conflict  with  the  notion  that 
northward  IMF  corresponds  to  a  “turning  off’  or  “turning 
down”  of  the  reconnection  process  and  that  it  requires  sev¬ 
eral  hours  for  the  ionospheric  polar  cap  potential  to  respond 
[Wygant  et  al.  [1983].  This  may  be  true  when  Bx  becomes 
near  0  or  slightly  positive,  but  for  large  values  of  northward 
IMF.  the  polar  cap  electric  field  appears  to  respond  rapidly. 
Following  the  initial  response  to  the  northward  IMF,  the  cur¬ 
rents  are  observed  to  evolve  for  about  22  min  before  achiev¬ 
ing  their  final  configuration. 

To  explain  the  apparent  control  of  the  observed  high- 
latitude  ionospheric  electric  field  configuration  we  have  used 
a  simple  model  in  which  the  dayside  high-latitude  field- 
aligned  current  distribution  is  controlled  by  the  strength 
and  orientation  of  the  IMF.  The  dayside  field-aligned  cur¬ 
rents  are  thought  to  connect  to  the  solar  wind  at  the  tops 
and  flanks  of  the  magnetosphere.  We  speculate  that  they 
connect  the  ends  of  magnetic  merging  lines  on  the  magne¬ 
topause  to  the  high-latitude  ionosphere.  The  DPY  and  DPZ 
systems  are  similar  to  that  proposed  by  Banks  et  al.  [1984] . 
The  DPY  system  is  thought  to  be  driven  by  the  north-south 
component  of  the  interplanetary  electric  field  (IEF)  shunted 
through  the  cleft  ionosphere  via  sheets  of  field- aligned  cur¬ 
rent  while  the  DPZ  system  is  driven  by  a  similar  shunting  of 
the  east-west  IEF  component  through  the  ionosphere.  The 
electric  connection  is  made  in  the  merging  region  poleward 
of  the  cleft  for  positive  Bx  and  equatorward  of  the  cleft  for 
negative  Bx . 

We  have  simplified  the  behavior  of  these  currents  by  fixing 
the  longitudinal  position  of  the  currents  and  varying  only 
the  latitudinal  position,  strength,  and  polarity  of  the  cur¬ 
rents  as  determined  by  the  time- varying  IMF.  Nevertheless, 
this  simple  representation  of  the  currents  controlled  by  the 
measured  IMF  reproduces  the  majority  of  features  seen  in 
the  convection  observations.  Additional  accuracy  might  be 
achieved  with  a  more  accurate  conductivity  model  and  with 
additional  parameters  to  adjust  the  longitudinal  positions  of 
the  cleft  currents.  We  believe,  however,  that  the  simulation 
demonstrates  that  many  of  the  complex  features  seen  in  the 
high-latitude  convection  are  a  reflection  of  where  magneto- 
spheric  merging  boundaries  connect  to  the  ionosphere.  This 
appears  to  be  true  for  both  positive  and  negative  IMF  Bx . 

In  this  simulation,  we  have  chosen  to  treat  the  region 
1  and  region  2  currents  in  a  simple  manner,  particularly 


on  the  night  hemisphere.  They  are  continuous  around  the 
oval,  having  a  distribution  similar  to  the  lijtma  and  Potemra 
[1976]  distribution.  The  location  and  strength  of  the  region 
1  and  2  currents  are  constant  for  negative  B,  For  Bx  =  0 
the  region  1  and  region  2  current  strength  is  reduced  and 
for  Bz  >  0  the  currents  move  poleward.  The  location  and 
strength  are  constant  for  Bx  >  0  and  there  is  a  step  change 
at  Bt  =  0.  The  continuous  distribution  of  region  1  and  2 
currents  is  a  simplification  based  upon  the  early  statistical 
analysis  of  satellite  observations.  Since  our  primary  interest 
here  concerns  the  dayside  convection  and  its  apparent  direct 
control  by  the  IMF,  we  feel  justified  in  this  simplification. 
Recent  studies,  however,  suggest  that  the  dayside  and  night- 
side  region  1  and  region  2  currents  may  be  separated  and 
may  be  produced  by  separate  sources  [Phis- Christensen  et 
al.,  1985;  Spence  et  al.,  1987].  Therefore,  while  we  feel  that 
our  representation  of  the  dayside  high-latitude  auroral  zone 
and  polar  cap  convection  is  explained  by  the  model  used  in 
this  paper,  we  make  no  such  claim  for  the  convection  ob¬ 
served  on  the  nightside  hemisphere. 

Finally,  there  has  been  some  debate  concerning  the  con¬ 
figuration  of  the  large-scale  convection  pattern  during  in¬ 
tervals  of  northward  IMF.  Positive  Bx  convection  models 
which  contain  three  or  four  convection  cells  [e.g.,  Potemra 
et  al.,  1984;  Re  iff  and  Burch,  1985;  Banks  et  al.,  1984]  as 
wetl  as  a  model  which  achieves  sunward  polar  cap  convection 
through  a  distorted  two-cell  pattern  [Heppner  and  Maynard, 
1987]  have  been  discussed. 

To  address  this  discussion,  we  have  found  that  several  fac¬ 
tors  influence  the  dayside  high-latitude  convection  pattern: 
first,  the  relative  strengths  of  the  IMF  By  and  Bx  compo¬ 
nents,  and  second,  the  duration  of  the  positive  IMF  Bx. 
As  seen  here,  it  required  over  20  min  for  the  high-latitude 
reversed  convection  to  be  established  during  an  interval  of 
very  strong  positive  Bx.  (This  is  in  contrast  to  the  more 
rapid  reconfiguration  of  flow  at  lower  latitudes,  which  un¬ 
dergo  180°  reversals  within  about  12  min  of  the  IMF  po¬ 
larity  change  encountering  the  magnetopause,  in  response 
to  changes  in  the  IMF  Bv  component  [Clouer  and  Banks, 
1986]).  Also,  we  find  that  the  region  of  reversed  convection 
associated  with  northward  IMF  is  confined  to  very  high  lat¬ 
itudes.  In  the  case  examined  here,  the  convection  reversal 
boundary,  which  one  usually  considers  to  be  at  or  near  the 
polar  cap  boundary,  was  located  between  71°  and  75°  prior 
to  the  northward  turning  of  the  IMF.  The  reversal  boundary 
appears  to  remain  at  those  latitudes  for  some  time  follow¬ 
ing  the  northward  turning  of  the  IMF.  Reversed  convection, 
however,  appears  to  be  confined  to  invariant  latitudes  above 
80°. 

Some  additional  insight  can  perhaps  be  obtained  from  an 
examination  of  Figure  5.  The  top  row  of  this  figure  shows 
model  potential  patterns  for  strong  northward  IMF.  These 
potential  patterns  are  congruent  with  the  streamlines  for 
the  convection  of  ionospheric  plasma.  The  center  pattern, 
for  IMF  By  =  0  shows  a  four-cell  pattern.  However,  the 
combination  of  strong  IMF  positive  Bx  and  strong  positive 
or  negative  IMF  By  produces  a  pattern  that  more  resem¬ 
bles  the  distorted  2-cell  model.  Thus,  the  appearance  of  a 
distorted  two-cell  system  or  a  four-cell  system  during  north¬ 
ward  IMF  may  depend  upon  the  relative  strengths  of  the 
IMF  components  and  the  associated  DPY  and  DPZ  current 
systems.  The  effect  of  the  DPY  system  is  to  distort  or  skew 
the  conventional  two-cell  pattern  on  the  dayside.  The  ef- 
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feet  of  the  DPZ  system  is  to  add  two  additional  convection 
cells  which  either  enhance  the  conventional  systems  on  the 
dayside,  or  are  opposite  to  the  conventional  systems  and 
produce  either  distortion  or  appear  separately  as  reversed 
convection  cells  depending  upon  the  relative  strengths  of 
the  separate  driving  currents. 

Unfortunately,  our  observations  alone  are  not  sufficient  to 
resolve  the  controversy.  To  do  so  would  require  simultaneous 
observations  at  high-latitudes  on  the  nightside  to  observe  the 
closure  of  the  convection  cells.  Observations  of  the  northern 
very  high-latitude  convection  pattern  observed  in  the  0000 
to  0600  LT  sector  during  northward  IMF  and  negative  By 
conditions  appear  to  be  the  crucial  item  to  test  the  validity 
of  the  distorted  two-cell  pattern  for  strong  northward  IMF 
conditions. 
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E  -  region  ion  convection  in  the  vicinity  of  the  dayside  polar  cleft  is  strongly  controlled  by  the 
interplanetary  magnetic  field,  suggesting  a  direct  current-driven  coupling  between  the  solar  wind 
and  dayside.  high-  latitude  ionosphere.  This  coupling  has  been  explored  through  observations  of 
high-latitude  ion  convection  measured  by  the  Sondre  Strom  fjord  Greenland  incoherent  scatter 
radar  and  subsequent  computer  simulations.  The  computer  simulation  calculates  the  ionospheric 
electric  potential  distribution  for  a  given  configuration  of  field- aligned  currents  and  conductivity 
distribution.  Using  a  simple  model  of  the  field-aligned  currents  linking  the  solar  wind  with  the 
dayside.  high-latitude  ionosphere,  we  can  explore  the  consequences  of  variations  of  the  currents 
on  ionospheric  plasma  convection.  The  direction  and  strength  of  the  currents  are  set  according 
to  measurements  in  the  upstream  solar  wind  and  the  time  varying  ionospheric  electric  field  and 
resulting  plasma  convection  are  simulated.  Comparisons  of  the  simulated  plasma  convection  with 
the  ion  velocity  measurements  at  Sondrestrom  show  very  good  agreement.  Thus,  a  very  simple 
model  which  provides  a  direct  current-driven  connection  between  the  solar  wind  and  the  high- 
latitude  ionosphere  via  field-aligned  currents  in  the  polar  cleft  is  able  to  explain  a  large  portion 
of  the  ionospheric  electric  field  variations  observed  in  the  vicinity  of  the  polar  cleft. 


1.  Introduction 

The  regions  in  and  around  the  dayside  polar  cleft  offer 
an  ideal  location  in  which  to  investigate  the  electrodynamic 
coupling  between  the  solar  wind,  magnetosphere,  and  iono¬ 
sphere.  Observations  of  current  systems,  electric  fields,  and 
plasmas  within  other  regions  of  the  auroral  zone  seem  to 
involve  many  intervening  interactions  which  obscure  iden¬ 
tification  of  the  original  processes  at  the  magnetospheric 
boundary.  The  polar  cleft,  however,  is  known  to  contain  so¬ 
lar  wind  plasma  and  the  phenomena  within  it  are  expected 
to  reveal  a  direct  and  detailed  account  of  the  interaction 
of  the  solar  wind  with  the  magnetosphere.  Investigation  of 
this  region  of  the  magnetosphere  bears  on  the  fundamental 
questions  of  whether  the  magnetosphere  is  open  or  closed, 
the  electrical  connectivity  between  the  solar  wind,  magne¬ 
tosphere  and  ionosphere,  and  how  solar  wind  plasma  enters 
the  magnetosphere. 

It  is  largely  for  this  reason  that  the  Chatanika  incoherent 
scatter  radar  was  moved  to  Sondre  Stromfjord,  Greenland 
[Banks  and  Evans.  1979].  Now  located  at  an  invariant  lat¬ 
itude  of  74°,  the  Sondre  Stromfjord  radar  is  able  to  make 
measurements  of  electrodynamic  and  plasma  properties  at 
the  ionospheric  end  of  the  northern  polar  cleft  on  a  routine 
basis.  This  report  presents  results  from  a  series  of  measure¬ 
ments  obtained  from  the  Sondre  Stromfjord  radar  together 
with  correlative  observations  of  the  interplanetary  magnetic 
field  (IMF)  and  solar  wind  plasma  obtained  by  the  IMP  - 
8  satellite.  These  experiments  were  designed  primarily  to 
examine  the  electrical  connection  between  the  solar  wind, 
magnetosphere,  and  ionosphere. 

It  is  well  known  that  the  high-latitude  ionospheric  con- 
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vection  pattern,  and  therefore  also  the  distribution  of  high- 
latitude  electric  fields  and  field-aligned  currents,  depends 
strongly  on  the  orientation  of  the  IMF  [Heehs.  1984:  Zanetti 
et  a/..  1984:  Burch  et  a/..  1985:  Reiff  and  Burch.  1985:  Frits- 
Chnstensen  et  al..  1985].  Ground  based  magnetic  observa¬ 
tions  at  high-  latitudes  provided  the  first  indications  that 
the  interplanetary  magnetic  field  was  an  important  order¬ 
ing  parameter  of  high-latitude  current  systems  [ Svalgaard . 
1973:  Mansurov.  1969:  Fnis-Chnstensen  et  al..  1972].  High- 
latitude  magnetic  observations  near  local  noon  were  found 
to  be  related  to  the  Bv  component  of  the  IMF  and  were 
designated  as  ''DPY'  perturbations  [Fnis- Christensen  and 
Wilhjelm.  1975].  Satellite  observations  have  shown  the  ex¬ 
istence  of  field-aligned  currents  which  could  feed  the  DPY 
current  system  in  the  vicinity  of  the  magnetospheric  cleft 
} Iijima  et  al..  1978:  McDiamud  et  al.,  1978].  Wilhjelm  et  al. 
[1978]  suggested  that  the  DPY  magnetic  variation  observed 
by  ground  observatories  is  the  result  of  an  ionospheric  Hall 
current  flowing  between  a  pair  of  oppositely  directed  field- 
aligned  current  sheets. 

The  complexity  of  the  total  field-aligned  current  system, 
however,  has  made  it  difficult  to  determine  the  exact  mor¬ 
phology  of  the  currents  and  current  sources  in  the  cleft.  This 
complexity,  as  seen  from  individual  satellite  passes,  led  Mc- 
Diarmid  et  al.  j  1978]  to  describe  the  cleft  as  a  “zone  of  con¬ 
fusion."  Some  hold  that  the  DPY  observations  are  merely 
the  result  of  a  redistribution  of  dayside  auroral  electrojets 
j Rostoker.  1980:  McDiarmid  et  al..  1979].  while  others  be¬ 
lieve  the  DPY  to  be  produced  by  a  separate  current  which 
can  vary  independently  with  respect  to  the  auroral  electro¬ 
jets  [Fnis- Christensen.  1981]. 

Satellite  observations  have  established  that  distinct  field- 
aligned  currents  exist  poleward  of  the  region  1  Birkeland 
currents  in  the  cleft  region  between  0930  and  1430  MLT 
[Iijima  and  Potemra.  197Go.  b:  Iijima  et  al..  1978 .Doyle  et 
nl..  1981).  It  is  known  that  the  currents  move  spatially  and 
they  appear  to  vary  in  strength  and  polarity  depending  upon 
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the  By  and  Bt  components  of  the  IMF.  The  exact  config¬ 
uration  of  the  currents  in  the  cleft  region,  however,  is  not 
yet  clearly  determined.  For  example.  Bythrow  et  al.  (1982] 
suggest  that  the  region  1  currents  observed  near  local  noon 
may  have  a  source  different  from  that  of  region  1  currents 
at  other  local  times.  Thus,  region  1  currents  near  noon  may 
be  part  of  the  cleft  current  system. 

Rich  and  K amide  [1983]  conducted  computer  modeling 
experiments  to  investigate  the  consequences  of  various  pat¬ 
terns  of  field  aligned  current  in  the  vicinity  of  the  dayside 
cleft.  They  found  that  patterns  of  overlapping  currents  sim¬ 
ilar  to  those  suggested  by  McDtarmid  et  al.  [1978]  and 
Doyle  et  ai  [1981]  better  reproduced  the  known  DPT  iono¬ 
spheric  current  behavior  than  patterns  without  overlapping 
current  at  noon  such  as  the  pattern  suggested  by  Iijima  and 
Potemra  [1976b].  The  overlapping  current  patterns  were 
achieved  through  a  reconfiguration  of  the  region  1  and  re¬ 
gion  2  currents,  whereas  the  nonoverlapping  current  config¬ 
uration  contained  a  seperate  cleft  system  in  which  the  equa- 
torward  portion  was  continuous  with  the  region  1  currents. 

Vennerstrom  et  al.  [1984]  compared  the  latitudinal  loca¬ 
tion  of  the  DPY  current  system  determined  from  the  West 
Greenland  meridian  chain  of  magnetometers  to  the  location 
of  the  convection  reversal  boundary  and  the  maximum  F  - 
region  electron  temperature  measured  by  the  Sondrestrom 
radar.  It  was  found  that  while  the  maximum  of  the  F  - 
region  electron  '"mpcatu'-p  roughly  coincides  «ith  the  con¬ 
vection  velocity  reversal  boundary,  the  DPY  current  is  al¬ 
ways  located  more  poleward.  The  auroral  electrojets  are  ob¬ 
served  to  be  equatorward  of  the  convection  reversal  bound¬ 
ary.  Thus.  Vennerstrom  et  al..  conclude  that  the  DPY  is 
indeed  a  current  system  separate  from  the  auroral  electro¬ 
jets  and  that  it  is  directly  associated  with  the  primary  solar 
wind  driven  convection. 

Models  to  explain  the  observed  DPY  variation  have  been 
proposed  by  several  researchers  [ D'Angelo .  1980;  Pnmdahl 
and  Spangslev,  1981.  1983;  Barbosa.  1979.  1983,  1984;  Banks 
et  al.,  1984],  In  these  models,  the  dayside  cleft  and  polar  cap 
contain  currents  which  are  directly  related  to  the  By  com¬ 
ponent  of  the  IMF.  Some  of  the  models  also  contain  currents 
which  are  related  to  the  Bt  components  of  the  IMF  and  pro¬ 
duce  DPZ  variations  which  may  also  be  related  to  the  occa¬ 
sional  occurrence  of  a  narrow  high-speed  convection  throat 
directed  in  the  antisolar  direction  [D’Angelo,  1980:  Banks  et 
al..  1984],  While  the  current  configurations  are  similar  in 
these  models,  they  vary  in  terms  of  the  generation  mecha¬ 
nism.  Figure  1,  adapted  from  D'Angelo  [1980]  and  Banks 
et  al.  [1984]  schematically  shows  the  configuration  of  IMF 
Bt  (top)  and  By  (bottom)  associated  field-aligned  currents 
in  the  vicinity  of  the  polar  cleft.  In  the  models  proposed  by 
Barbosa  [1979],  D'Angelo  [1980]  and  by  Banks  et  al.  [1984], 
the  currents  in  the  polar  cleft  are  driven  by  the  interplane¬ 
tary  electric  field  apparent  at  the  magnetopause  as  the  mag¬ 
netized  solar  wind  plasma  flows  past  the  earth.  While  the 
exact  nature  of  the  electrical  connection  between  the  solar 
wind  and  magnetosphere  is  not  addressed  in  detail  in  these 
models,  it  is  shown  that  the  models  can  reproduce  most  of 
the  observed  features  in  the  vicinity  of  the  cleft  and  that 
the  magnitudes  of  the  fields  and  currents  are  well  within  the 
realm  of  possibility.  Pnmdahl  and  Spangslev  [1981;  1983], 
and  Barbosa  [1983]  consider  a  somewhat  different  approach 
for  producing  cleft  currents.  In  this  mechanism,  rather  than 


a  direct  electrical  connection  with  the  solar  wind  Bx.  the 
cleft  currents  (bottom  of  Figure  1)  are  part  of  the  magne¬ 
topause  shielding  currents  related  to  By  variations  in  the 
IMF. 

Due  to  the  complexity  of  the  observed  field-aligned  cur¬ 
rents  in  the  vicinity  of  the  dayside  polar  cleft,  it  has  not  been 
possible  to  verify  the  configuration  of  currents  sketched  in 
Figure  1.  Through  computer  simulation,  however,  one  can 
explore  the  effects  on  ionospheric  electric  fields  of  different 
current  distributions  in  the  cleft  [Rich  and  K amide,  1983]. 
These  can  then,  in  turn,  be  compared  with  observations  of 
convection  obtained  by  high-latitude  radars.  In  this  paper, 
the  model  discussed  by  Banks  et  al.  [1984]  has  been  adapted 
to  adjust  the  strength  and  direction  of  the  DPY  and  DPZ 
field-aligned  currents  in  the  polar  cleft  based  on  IMF  mea¬ 
surements  obtained  by  the  IMP  -  8  satellite  upstream  in 
the  solar  wind.  Then,  through  a  series  of  steady  state  cal¬ 
culations,  it  is  possible  to  simulate  the  temporally  varying 
ionospheric  convection  expected  during  a  particular  inter¬ 
val  based  on  this  model.  The  simulated  Sondrestrom  radar 
convection  observations  can  then  be  directly  compared  with 
the  actual  radar  observations. 

We  note  that  while  we  use  values  of  the  IMF.  our  model 
presumes  that  the  interplanetary  electric  field  is  the  impor¬ 
tant  driving  parameter.  Apart  from  the  solar  wind  velocity, 
one  should  realize  that  IMF  correlations  are  equivalent  to 
correlations  with  the  interplanetary  electric  field. 

2  Observations 

Measurement  Procedure  and  Data  Analysis 

The  radar  data  which  have  been  analyzed  were  collected 
using  pairs  of  stationary  antenna  positions  on  either  side  of 
the  Sondre  Stromfjord  magnetic  meridian.  The  analysis  uses 
line-of-sight  ion  velocities  measured  in  the  F  -  region.  At  F 
-  region  altitudes  collisions  are  unimportant  and  the  ion  ve¬ 
locity  is  assumed  to  be  purely  the  result  of  E*  B  drift.  For 
each  pair  of  positions  the  measured  line-of-sight  velocities 
are  reduced  to  obtain  the  N-S  and  E-W  components  of  the 
ion  drift  velocity  as  a  function  of  latitude. 

Reduction  of  the  line-of-sight  velocity  measurements  from 
the  pairs  of  antenna  positions  to  obtain  a  velocity  vector 
orthogonal  to  the  local  magnetic  field  depends  upon  the  as¬ 
sumption  that  the  velocity  is  uniform  over  the  spatial  range 
of  the  measurements.  If  this  assumption  is  not  valid,  that 
is  if  there  is  a  gradient  or  a  rotation  in  the  ion  flow  from 
one  measurement  position  to  the  second,  the  analysis  proce¬ 
dure  can  over  estimate  the  E-W  component  of  velocity  and 
underestimate  the  N-S  component  of  velocity. 

Assuming  a  vertical  magnetic  field  for  simplicity  of  illus¬ 
tration.  the  radar-deduced  ion  velocity  is  derived  from  the 
pair  of  measurements  in  the  following  manner:  Let  li  and 
>2  be  unit  vectors  in  the  antenna  azimuth  pointing  direction 
for  the  two  radar  positions.  Let  V  be  the  velocity  which  is 
uniform  over  the  two  measurement  points  Pi  and  Pj.  The 
magnitudes  of  the  measured  line-of-sight  velocities  O)  and 
«2  at  the  respective  points  are  given  by, 

»i=l|V 

and 

a2  =  12  V 


6G 


CLAUER  AND  BANKS:  CLEFT  CU  RRENTS.  OBSERVATIONS  AND  SIMULATION 


6%  I 
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B 

Fig.  I.  A  general  view  of  the  way  cleft  currents  driven  by  the 
interplanetary  electric  field  connect  across  the  magnetosphere  In 
Figure  1  a,  line  currents  arising  from  £„  connect  from  the  Banks 
of  the  magnetosphere  to  two  discrete  zones  on  either  side  of  local 
magnetic  noon  The  direction  of  the  current  is  determined  by  the 
sign  of  B:  and  the  case  for  Bz  >  0  is  shown.  In  figure  16  cur¬ 
rents  arising  from  Ez  connect  to  each  cleft  through  sheet  currents 
centered  at  local  magnetic  noon.  Completion  of  the  current  path 
through  the  opposite  hemisphere  is  accomplished  by  means  of  a 
sheet  current  along  closed  field  lines  lying  at  or  near  the  mag¬ 
netopause  The  direction  of  currents  shown  here  corresponds  to 
By  <  0.  Oppositely  directed  current  would  result  for  By  >  0. 
After  Banks  ft  at.  [198-11 


l\E 

h  n 

Ve 

_  ra* 

h  E 

h  n 

Yn 

l«2  . 

where  the  subscripts  E  and  N  refer  to  eastward  and  north¬ 
ward  components  respectively 

The  solutions  to  this  equation  are  given  by 


'fT  = 


-  hNCt2 

llEl2N  -  h\h E 


VN  = 


^1EQ2  ~  ^2Ea\ 
llEhN  ~  ll\h E 


We  have  examined  the  velocity  computations  using  a  qual¬ 
itative  method  to  estimate  the  validity  of  the  assumption 
that  the  velocity  be  uniform  over  the  spatial  dimensions 
of  the  measurement  points.  This  was  done  by  conducting 
an  experiment  using  four  beam  positions  at  30°  elevation 
in  the  northward  direction.  The  azimuths  of  the  positions 
were  -26°.  -13°.  +13°.  and  +26°  measured  from  magnetic 
north.  Velocities  were  computed  separately  using  the  inner 
pair  and  the  outer  pair  of  positions.  This  experiment  was 
run  on  April  25  and  April  27.  1983  during  period  of  moder¬ 
ate  geomagnetic  activity  (kp  =  3-  to  4). 

Before  presenting  the  results  of  these  experiments,  how¬ 
ever.  it  is  useful  to  consider  the  effects  of  a  nonuniformity 
in  the  flow.  Suppose  that  the  velocity  at  the  westward  mea¬ 
surement  point  is  400  m/s  to  the  east  and  the  velocity  at 
the  eastward  measurement  point  is  400  m/s  to  the  north. 
Such  a  situation  would  result  from  a  90°  rotation  in  the 
flow  between  the  measurement  points.  For  the  inner  pair 
of  positions.  f,v  =  0.9744  and  Ie  =  0.2250  which  lead  to 
ai  =  389.76  m/s  and  02  =  -90.0  m/s.  Then. 


389.76  +  90  0 
2.(0.2250) 


1066.1  in/s 


389.76  -  90.0 
2.(0.9744) 


153.8  m/s 


Similarly,  for  the  outer  positions  one  obtains  Ve  =  609 
m/s  and  V'v  =  102  m/s.  Thus,  the  eastward  velocity  is 
greatly  overestimated  and  both  the  inner  and  outer  pair  of 
points  give  essentially  an  eastward  velocity  but  of  different 
magnitudes. 

Suppose,  for  a  second  example,  that  the  velocity  at  the 
westward  measurement  point  is  300  m/sec  to  the  north  and 
at  the  eastward  measurement  point  is  600  m/s  northward. 
Thus,  there  is  a  gradient  in  the  flow.  Using  the  inner  points 
yields  Ve  =  649.5  m/s  and  V/v  =  501.3  m/s  and  using  the 
outer  points  yields  Ve  =  307.  m/s  and  V)y  =  450  m/s. 
Again  a  substantial  eastward  velocity  is  introduced  into  the 
result. 

Plates  la  and  lb  show  the  results  of  the  experiments  which 
were  performed  on  April  25  and  April  27.  The  F  -  region 
ion  velocities  orthogonal  to  B  are  plotted  as  a  function  of  in¬ 
variant  latitude  and  universal  time  (shown  around  the  inner 
circumference  of  the  plot).  Local  time  is  shown  around  the 
outer  circumference  of  the  plot  with  magnetic  local  noon  at 
the  top  of  the  plot.  The  black  vectors  were  obtained  using 
the  interior  antenna  positions  while  the  red  vectors  were  ob¬ 
tained  using  the  outer  antenna  positions.  The  sets  of  vectors 
obtained  from  the  inner  and  outer  pairs  of  measurement  po¬ 
sitions  have  been  superposed  for  comparison.  For  the  most 
part,  the  agreement  is  very  good.  However,  there  are  oc¬ 
casional  times  when  there  is  considerable  disagreement  be¬ 
tween  the  velocities  obtained  using  the  two  pairs  of  points, 
thus  indicating  that  occasionally  the  assumption  of  uniform 
velocity  breaks  down  It  is  not  clear  whether  this  is  the  re¬ 
sult  of  a  spatial  difference  in  the  velocities  or  the  result  of  a 
temporal  change. 
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IMP-8  MAGNETIC  FIELD 
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Fig  2  Interplanetary  magnetic  field  observations  shown  in  GSM  coordinates  for  (o)  July  22.  1983  and  (6)  June 
28.  1983. 
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Nevertheless,  the  good  general  agreement  between  the 
black  and  red  velocity  vectors  indicates  that  the  multipo¬ 
sition  method  of  determining  the  vector  ion  velocity  is  ac¬ 
curate.  There  are  occasional  small-  scale  features  in  the 
real  ion  velocity  which  cause  errors  in  the  velocity  deter¬ 
mination.  However,  the  points  which  we  shall  discuss  in 
this  paper  should  be  unaltered  by  the  uncertainties  in  the 
computed  velocities. 

During  the  course  of  our  experiments  we  used  several  vari¬ 
ations  of  the  positions  and  dwell  times.  These  changes  af¬ 
fected  the  cycle  time  of  the  observations,  the  latitudinal 
range  covered  and  the  uncertainty  in  the  line-of-sight  ve¬ 
locity  measurements.  Additional  details  of  the  radar  capa¬ 
bilities  and  analysis  techniques  can  be  found  in  W  tcku  aT  et 
al.  [I984j. 

Correlative  observations  of  the  solar  wind  speed,  density 
and  magnetic  field  are  provided  by  the  IMP  -  8  spacecraft 


These  data  together  with  the  observations  of  high-latitude 
ion  convection  obtained  using  the  Sondrestrom  radar  are 
studied  to  determine  the  solar  wind  control  of  high-  latitude 
dayside  currents  and  fields.  The  analysis  presented  here 
extends  the  preliminary  results  reported  by  Clauer  et  al. 

;  1 984]  and  Jorgensen  et  al.  [  1984]  . 

Data  Presentation 

It  is  well  known  that  variations  in  the  Y  component  of 
the  IMF  are  associated  with  changes  in  the  high-latitude 
ionosphrric  convection  pattern  [ Heelis .  1984).  This  is  par¬ 
ticularly  true  in  the  dayside  polar  cleft  region.  Plates  2a  and 
26  show  observations  of  F  -  region  ion  velocities  obtained  us¬ 
ing  tile  .Sondrestrom  radai.  Correlative  measurements  of  the 
IMF  obtained  from  the  IMP  -  8  satellite  are  shown  in  figures 
'In  and  26  On  July  22.  1983.  shown  in  Plate  2a.  the  IMF  By 
<  omponent  was  negative  for  most  of  the  observing  interval. 
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Fig.  3.  Histogram  showing  time  delay  from  when  a  change  in 
sign  of  IMF  By  encounters  the  magnetopause  to  the  time  when  a 
change  in  polarity  of  the  east  -  west  F-  region  ion  drift  is  observed 
in  the  vicinity  of  the  cleft. 

while  on  June  28,  1983,  shown  in  Plate  2b.  the  By  compo¬ 
nent  was  primarily  positive.  The  difference  in  the  sign  of  By 
appears  to  control  the  direction  of  the  dayside  high-latitude 
convection  with  Bv  >  0  associated  with  westward  convec¬ 


tion  and  By  <  0  associated  with  eastward  convection.  The 
magnitudes  of  the  By  and  Bz  components  are  also  different 
on  these  days  and  this  is  associated  with  the  differences  in 
magnitude  of  the  ion  velocities  and  the  latitudinal  location 
of  the  convection  reversal  boundary  on  these  two  days. 

The  basic  behavior  of  the  high-latitude  currents  has  pre¬ 
viously  been  described  based  on  magnetometer  observations 
[Frtis-Christensen.  1981;  Rostoker.  1980],  The  electric  field 
observations  (orthogonal  to  the  ion  velocity  measurements) 
confirm  the  hypothesis  that  the  DPY  current  system  is  a 
Hall  current  ( Frits- Christensen.  1986]. 

ClaueT  et  a l.  [1984]  showed  that  the  dayside  high-latitude 
ionospheric  electric  field  responds  very  directly  to  changes  in 
the  IMF  By  component.  Polarity  changes  in  the  east-west 
ion  convection  were  measured  to  occur  within  about  10  or 
15  minutes  of  the  time  that  an  IMF  By  polarity  change  was 
estimated  to  encounter  the  magnetopause.  Figure  3  shows 
a  histogram  which  displays  the  delay  times  obtained  for  all 
of  the  clear  correlated  polarity  changes  found  in  our  data 
base.  This  data  base  contains  264  hours  of  radar  observa¬ 
tions  in  the  vicinity  of  the  polar  cleft  during  1983  and  1984. 
For  events  to  be  included,  we  have  required  that  the  IMF 


FIELD  ALIGNED  CURRENT 


Fig.  4.  Polar  plot  of  the  distribution  of  field-aligned  current  used  in  the  computer  model.  Three  regions  of 
current  are  shown:  (1)  the  region  1  and  region  2  currents  around  the  auroral  oval  with  inward  current  denoted  by 
shading  (2)  the  DPY  sheet  currents  centered  at  local  magnetic  noon  and  lying  poleward  of  the  region  1  currents, 
and  (3)  the  DPZ  line  currents  at  the  longitudinal  ends  of  the  DPY  currents.  The  polarities  of  the  currents  in  ibis 
illustration  correspond  to  the  case  of  IMF  By  >0  and  Bz  >  0. 
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change  in  By  measured  by  IMP  -  8  in  the  upstream  solar 
wind  have  a  constant  polarity  for  20  minutes  on  both  sides 
of  the  sign  change.  We  have  then  estimated  the  delay  time 
between  the  satellite  and  the  magnetopause.  This  is  done 
by  first  computing  the  subsolar  magnetopause  standoff  dis¬ 
tance  based  on  a  balance  between  the  solar  wind  dynamic 
pressure  and  magnetospheric  magnetic  pressure.  The  tran¬ 
sit  time  between  IMP  -  8  and  the  magnetopause  is  computed 
using  the  measured  solar  wind  velocity.  An  additional  delay 
of  four  minutes  is  added  to  account  for  the  magnetosheath 
transit  time  based  on  the  gasdynamic  model  calculations  of 
Spretter  and  Stahara  [1980].  Two  calculations  are  made:  (1) 
assuming  that  the  plane  of  the  IMF  By  polarity  change  lies 
orthogonal  to  the  earth-sun  line,  and  (2)  assuming  that  it 
lies  along  the  Parker  spiral  angle.  This  gives  two  estimates  of 
the  delay  time  between  the  satellite  and  the  magnetopause 
which  can  be  quite  similar  or  quite  different  depending  on 
the  position  of  the  satellite.  For  the  purpose  of  the  plot,  we 
have  taken  the  average  of  the  computed  times.  The  time 
delay  for  the  high-latitude  ionospheric  convection  to  change 
polarity  is  then  measured  relative  to  the  estimated  time  that 
the  By  sign  change  encounters  the  magnetopause. 

While  we  have  264  hours  of  observation  time  at  Sondre 
Stromfjord,  the  relatively  few  observations  shown  in  the  his¬ 
togram  reflects  the  inconsistent  tracking  of  the  IMP  -  8  satel¬ 
lite  together  with  our  selection  criterion. 

The  peak  of  the  distribution  shown  in  Figure  3  is  at  a 
delay  time  of  about  14  minutes.  If  the  change  in  current 
to  the  ionosphere  from  the  magnetopause  is  carried  by  a 
tiaiii .efae  Alfven  wave,  this  suggests  that  several  bounces 
of  the  wave  probably  occur  since  the  Alfven  transit  time  is 
on  the  order  of  2  minutes  or  less.  This  may  be  reasonable 
since  the  dayside  ionosphere  would  be  a  good  reflector  of  the 
wave. 

Model  Simulation 

As  mentioned  in  the  introduction,  the  model  presented  by 
Banks  et  a l.  [1984]  has  been  adapted  to  perform  a  series  of 
computations  of  the  ionospheric  electric  field  based  on  mea¬ 
surements  of  the  IMF  obtained  by  the  IMP  -  8  spacecraft 
Using  these  computations,  we  may  simulate  the  ion  convec¬ 
tion  velocity  which  would  be  observed  by  the  Sondrestrom 
radar  under  the  conditions  described  by  the  model. 

In  Figure  4  we  show  a  polar  projection  of  the  distribution 
of  currents  which  is  input  to  the  model.  The  current  dis¬ 
tribution  includes  the  region  I  and  region  2  currents  using 
the  approximate  statistical  distribution  and  spatial  arrange¬ 
ment  from  hjtma  and  Potemra  [1976a]  and  the  suggested 
DPY and  DPZ  field-aligned  currents  in  the  polar  cleft  from 
Banks  et  al.  [1984].  This  figure  shows  the  current  flow  for 
the  case  when  Bz  >  0  and  By  >  0  When  Bz  —  By  =  0 
no  cleft  currents  would  flow  and  for  Bz  <  0  and  By  <  0  the 
polarities  of  the  cleft  currents  would  reverse.  In  all  of  the 
simulations  presented  here,  the  total  region  1  inward  and 
outward  currents  are  each  1.7  x  106  A  and  the  total  region 
2  inward  and  outward  currents  are  each  15*  lit6  A.  The 
simulations  use  the  conductivity  model  from  Kamide  and 
Matsushita  [1979]  with  a  conductivity  enhancement  in  the 
auroral  oval.  The  peak  potential  difference  across  the  polar 
cap  resulting  from  the  region  1  and  region  2  currents  alone  is 
128  kV.  The  region  1  and  2  currents  are  constant  in  location 
and  strength  in  all  of  our  simulations.  The  only  variations 
are  in  the  direction  and  strength  of  the  DPY  and  DPZ  cur¬ 
rents  in  the  cleft.  The  direction  of  the  DP  V  currents  depends 
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upon  the  sign  of  the  IMF  By  component.  The  magnitude 
of  the  current  varies  linearly  between  1.00*  10®  *  A  and 
0  A  for  ] By |  =  8  nT  and  0  nT  respectively.  The  currents 
saturate  at  1.00  ■  10®  A  for  larger  absolute  values  of  By. 
Similarly,  the  direction  of  the  DPZ  currents  depends  upon 
tl  '  sign  of  the  IMF  Bz  component.  The  magnitude  of  the 
DPZ  current  varies  linearly  between  3.26  *  10s  A  and  0  A 
for  | B 2 1  —  8  nT  and  0  nT  respectively.  The  currents  satu¬ 
rate  at  3.26  *  10s  A  for  larger  absolute  values  ol  Bz.  The 
location  of  the  cleft  currents  remains  constant  in  the  model 
simulations. 

For  each  measurement  time  at  Sondrestrom,  we  take  the 
measured  IMF  at  IMP  -  8  using  the  appropriate  delay  lime 
determined  from  the  actual  solar  wind  conditions  during  the 
event.  Based  on  these  IMF  values,  we  adjust  the  strength 
and  direction  of  the  DP)'  and  DPZ  currents  and  compute 
the  hemispheric  electric  field  distribution.  Using  the  com¬ 
puted  model  electric  field,  the  E  *  B  drift  at  the  radar  mea¬ 
surement  point  is  determined  and  the  simulated  vector  ve¬ 
locity  is  plotted. 

Plate  3a  and  36  shows  the  measured  and  simulated  convec- 
tion  velocities  for  April  25.  1983.  respectively,  while  Figure 
5  shows  the  IMF  observations  for  this  same  interval.  The 
IMF  measurements  show  a  gap  between  1113  and  1207  UT. 
However,  from  the  Greenland  meridian  chain  of  magnetome¬ 
ters.  one  can  infer  that  the  By  was  negative  for  most  of  that 
interval.  For  the  model  simulations,  therefore,  we  have  used 
the  values  of  the  IMF  measured  at  1207  UT  throughout  the 
data  gap 

In  comparing  Plate  3a  with  3b.  the  overall  agreement  must 
be  said  to  be  very  good.  Some  of  the  vectors  at  lower  lati¬ 
tudes  on  the  evening-side  do  not  match  well  and  the  location 
of  the  cleft  currents  in  the  model  appears  to  be  at  too  low 
a  latitude,  indicated  by  the  poleward  velocities  of  the  most 
northern  vectors  in  the  simulation. 

The  most  important  feature  in  Plate  3a  is  the  interval  be¬ 
tween  1230  UT  and  1400  UT  where  the  high-latitude  flow 
direction  changed  by  about  180  degrees  compared  with  the 
flows  on  either  side  of  this  interval.  Clauer  et  al.  [1984] 
suggested  that  these  two  changes  in  ion  flow  direction  were 
directly  related  to  changes  in  the  polarity  of  the  IMF  Bv 
component  at  1215  UT  and  1343  UT  respectively.  This  sug¬ 
gestion  appears  to  be.  for  the  most  part,  verified  by  the  re¬ 
sults  of  the  simulation  shown  in  Plate  36.  The  reversal  in  the 
east  -  west  ion  velocities  is  observed  between  1230  UT  and 
1400  UT.  There  is.  however,  a  set  of  poleward  velocities  at 
1300  UT  which  appear  to  be  associated  with  the  southward 
fluctuation  of  the  IMF  at  1249  to  1300  UT.  Such  poleward 
velocities  were  not  observed  in  the  actual  measurements. 
Perhaps  the  sampling  frequency  was  too  slow,  thus  caus¬ 
ing  us  to  miss  sampling  the  poleward  flow,  or  perhaps  the 
10-minute  southward  fluctuation  was  of  insufficient  duration 
to  establish  the  large-scale  cleft  current  system  which  links 
the  solar  wind  with  the  ionosphere.  The  magnetosphere 
has  heen  modeled  as  a  low-pass  filter  which  only  couples 
variations  of  the  solar  wind  electric  field  with  frequencies 
less  than  It)-'*  Hz  to  auroral  current  systems  [ Clauer  et  al.. 
1081.  1983j.  Perhaps  the  cleft  currents  also  respond  only  to 
variations  in  the  IMF  below  some  critical  frequency. 

Plate  4a  and  46  show  another  example  of  observations  and 
model  simulation.  March  16.  1984  is  a  particularly  interest¬ 
ing  day  because  of  the  extended  region  of  poleward  flow 
from  1.300  UT  through  1445  UT  This  region  of  strong  pole- 
ward  How  is  associated  with  a  large  southward  IMF  shown 


CLAUER  AND  BANKS.  CLEFT  CURRENTS.  OBSERVATIONS  AND  SIMULATION 


6965 


IMP-8  MAGNETIC  FIELD 

25- APR-83 

11:00:00  TO  10:00:00 


11  12  13  14  IS  W  17  18 

Universal  Time  (hours) 


Fig.  5.  Interplanetary  magnetic  field  on  April  25.  1983,  in  GSM  coordinated  measured  by  the  IMP  -  8  satellite. 
The  traces,  from  the  top,  show  the  Bz,  By,  B,  components,  and  the  total  field. 


in  Plate  6.  In  this  and  the  following  simulation,  simulated  One  region  of  disagreement  appears  at  about  A  <  75°  near 
convection  velocities  are  not  computed  during  times  when  local  noon.  This  may  suggest  some  different  configuration  of 
IMP  -  8  measurements  are  not  available.  We  find  that  the  the  dayside  auroral  electrojets  and  the  region  1  and  region 
simulation  reproduces  the  observed  poleward  flow  very  well,  2  currents  at  this  time. 

thus  suggesting  that  DPZ  currents  similar  to  those  of  the  A  final  example  is  shown  in  Plates  5a  and  56  with  the 
model  flow  during  southward  Bz  conditions.  Further  com-  IMF  observations  shown  in  Figure  7.  July  24,  1983  showed 
parison  of  Plate  4a  with  46  also  reveals  very  good  agreement,  more  complex  temporal  variations  than  the  previous  exam- 
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Fig  7  Interplanetary  magnetir  field  on  July  24.  1983. 


pies  and  is  a  very  interesting  case  to  attempt  to  simulate 
with  our  simple  model.  The  results  are  very  interesting. 
From  1230  UT  to  1330  UT  a  region  of  strong  poleward  flow 
is  observed  in  the  ion  convection.  This  is  associated  with  a 
strong  southward  IMF  as  were  the  observations  in  the  pre¬ 
vious  example.  However,  when  one  examines  the  timing  of 
the  southward  variation  and  the  onset  of  poleward  ion  flow, 
a  problem  appears  to  arise.  Bz  turns  southward  at  11:15 
UT  and  the  expected  delay  time  between  the  satellite  and 
magnetopause  is  between  8  and  12  minutes.  Thus,  the  time 
delay  from  when  the  southward  fluctuation  encounters  the 
magnetopause  and  the  onset  of  poleward  flow  is  over  50  min¬ 
utes  and  is  inconsistent  with  the  histogram  shown  in  Figure 
5.  Further  examination  of  the  IMF  variations  indicates  that 
was  largo  and  positive  at  the  time  B z  turned  southward. 
By  became  nearly  0  at  about  1300.  This  suggests  the  pos¬ 
sibility  that  the  relative  magnitudes  of  By  and  Bz  may  be 
important  in  determining  the  ionospheric  effects  of  the  DPY 
and  DPZ  currents.  Alternatively,  perhaps  the  effects  of  the 
DPZ  current  were  not  observed  around  1140  UT  because  of 
a  shift  in  the  spatial  distribution  of  the  DPZ  currents. 

It  is  very  interesting  to  see  that  the  simulation  of  ion  ve¬ 
locity  measurements  agrees  very  well  with  the  observations. 
The  interval  of  poleward  flow  between  1230  and  1330  UT 
is  reproduced.  We  note  that  the  simulation  only  varies  the 
strength  of  the  cleft  currents  keeping  the  location  constant 
Therefore,  the  results  of  the  simulation  indicate  that  adjust¬ 
ing  the  relative  strengths  of  the  DPY  and  DPZ  currents  is 
sufficient  to  explain  the  observations.  Following  1330  UT. 
the  simulation  and  measured  velocities  do  not  agree.  This 
region  of  disagreement  appears  to  occur  during  an  interval  of 
northward  IMF.  This  suggests  that  it  is  not  sufficient  to  sim¬ 
ply  reverse  the  polarity  of  the  DPZ  rurrent  during  intervals 
of  strongly  northward  IMF.  Other  regions  of  disagreement 
occur  at  lower  latitudes,  particularly  on  the  afternoonside 
and  at  local  noon.  Otherwise,  we  find  generally  good  agree¬ 


ment  between  the  results  of  the  simulation  and  the  observed 
ionospheric  convection  on  July  24.  1984. 

3.  Discussion  and  Conclusions 

This  report  has  described  experimental  observations  of 
high-latitudinal  ion  convection  using  the  Sondre  Stromfjord. 
Greenland,  incoherent  scatter  radar.  The  technique  for  mea¬ 
suring  F  -  region  ion  velocities  at  high  time  resolution  over  a 
large  range  of  latitudes  has  been  described  and  experiments 
have  indicated  that  the  technique  provides  reliable  measure¬ 
ments  of  the  ion  velocity  orthogonal  to  the  local  magnetic 
field  under  most  circumstances.  Small-scale  spatial  varia¬ 
tions.  however,  will  tend  to  cause  errors  in  the  computed 
velocity  vertor.  and  such  variations  probably  occur  at  some 
times.  Nevertheless,  our  results  indicate  that  the  pattern  of 
convection  obtained  from  our  measurements  is.  in  general, 
accurate. 

Correlative  observations  of  solar  wind  plasma  parameters 
and  IMF  from  the  IMP  -  8  satellite  have  also  been  stud¬ 
ied  These  data  indicate  that  the  dayside  high-latitude  iono¬ 
spheric  electric  field,  which  drives  the  observed  F  -  region 
plasma  convection,  responds  within  about  14  minutes  to 
IMF  variations  at  the  magnetopause.  There  is  also  indi¬ 
cation  that  IMF  fluctuations  of  10  minutes  or  less  do  not 
produce  changes  in  the  ionospheric  electric  field  configura¬ 
tion 

Finally,  we  have  found  that  the  use  of  a  very  simple  model 
of  the  polar  cleft  currents  driven  by  solar  wind  electric  fields 
is  able  to  explain  much  of  the  temporal  variations  in  the 
ionospheric  plasma  convection  pattern  observed  using  the 
Sondrestrom  radar.  In  addition  we  find  several  areas  of 
disagreement  which  will  require  refinement  of  the  simple 
model.  F'or  example,  observations  during  intervals  of  north¬ 
ward  IMF  are  not  consistent  with  the  simulations.  Also,  oc¬ 
casionally  the  convection  equatorward  of  the  cleft  currents 
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in  the  model  is  inconsistent  with  observations,  particularly 
near  local  noon.  This  suggests  that  the  region  1  and  region 
2  currents  near  noon  may  not  be  modeled  properly. 

When  beginning  these  simulations,  it  was  expected  that 
the  model  would  be  an  oversimplification  of  the  actual  cur¬ 
rent  systems  in  the  cleft.  For  example,  we  do  not  move  the 
currents  in  local  time  in  response  to  IMF  By  variations  as 
suggested  by  the  observations  of  Doyle  el  al.  [1981]  and  the 
antiparallel  merging  model  [ Crooker .  1979:  Ltihmann  el  al.. 
1984j.  In  addition,  the  currents  in  our  simulation  do  not 
move  in  latitude  as  one  would  expect  if  the  polar  cap  grows 
or  shrinks  in  response  to  IMF  Bz  variations.  We  have  made 
the  DPY  current  system  separate  and  above  the  region  1 
and  region  2  currents.  However,  as  suggested  by  Bythrow  et 
al..  [1982].  the  equatorward  portion  should  perhaps  overlap 
the  region  1  currents.  Finally,  perhaps  a  different  configu¬ 
ration  of  currents  during  northward  IMF  such  as  the  NBZ 
currents  suggested  by  Iijima  et  al.  [1984 j  is  required. 

Some  of  these  modifications  may  correct  some  of  the  in¬ 
consistencies  noted  above.  Nevertheless,  it  is  quite  remark¬ 
able  that  the  simulations  reproduce  the  observations  as  well 
as  they  do.  This  seems  to  indicate  that  the  structure  of  the 
currents  in  the  polar  cleft  is  indeed  rather  similar  to  that 
of  the  model.  It  also  argues  strongly  that  the  currents  are 
indeed  directly  driven  by  the  solar  wind. 
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The  Technique  of  Linear  Prediction  Filters  Applied  to  Studies  of 
Solar  Wind-Magnetosphere  Coupling 

C.  Robert  Clauer 

Space.  Telecommunications  and  Radioscience  Laboratory.  Stanford  Electronics  Laboratories 
Stanford  University.  Stanford,  California  94305,  U  S  A. 

Linear  prediction  filtering  is  a  powerful  empirical  technique  suitable  for  the 
study  of  stimulus-response  behavior.  The  technique  enables  one  to  determine  the 
most  general  linear  relationship  between  multiple  time-varying  quantities,  assuming 
that  the  physical  systems  relating  the  quantities  are  linear  and  time  invariant.  Several 
researchers  have  applied  linear  prediction  analysis  to  investigate  solar  wind- 
magnetosphere  interactions.  This  short  review  describes  the  method  of  linear 
prediction  analysis,  its  application  to  solar  wind-magnetosphere  coupling  studies 
both  in  terms  of  the  assumptions  and  in  terms  of  physical  processes,  and  the  results  of 
investigations  which  have  used  this  technique. 


1.  Introduction 

A  great  deal  of  scientific  effort  has  been  invested  over  the  past  several  decades  to 
understand  the  causes  of  geomagnetic  activity.  It  is  now  known  that  the  sun  through 
its  expanding  outer  atmosphere,  the  solar  wind,  is  the  primary  source  of  the  energy 
which  produces  geomagnetic  disturbances.  These  disturbances  recorded  at  our 
geomagnetic  observatories  result  from  electrical  currents  which  flow  around  the  earth 
in  space,  through  the  ionosphere  and  within  the  conducting  earth.  The  magnetic 
effects  produced  by  these  currents  are  characterized  by  various  simple  measures  called 
geomagnetic  indices  (Rostoker,  1972;  Mayaud,  1980;  Kamide  and  Akasofu, 
1983).  Careful  understanding  of  the  methods  used  to  create  indices  and  the 
consequent  (imitations  of  the  indices  is  required  when  making  a  physical  interpreta¬ 
tion  based  on  their  use.  The  benefit  of  the  indices,  however,  is  to  synthesize  a  variety  of 
complex  geomagnetic  disturbance  measurements  into  a  simple  time  series.  Because  of 
this  simplification,  geomagnetic  indices  have  been  used  extensively  in  studies  of 
magnetospheric  behavior.  Our  understanding  of  the  solar  wind’s  control  over 
geomagnetic  activity  has  been  largely  acquired  through  studies  based  on  correlations 
between  solar  wind  quantities  and  geomagnetic  indices.  These  correlations  have 
formed  the  framework  within  which  theories  of  solar  wind-magnetosphere  coupling 
have  been  developed. 

Most  correlative  studies  have  used  regression,  superposed  epoch,  or  cross 
correlation  techniques  to  establish  relationships  between  the  solar  wind  and 
geomagnetic  activity  and  a  discussion  of  these  techniques  is  given  by  Baker  (1986)  in 
this  volume.  Recently  another  very  powerful  technique,  that  of  linear  prediction 
filtering,  has  been  introduced  to  the  field  of  magnetospheric  physics  (IYEMORI  et  a! ., 
1979;  IYEMORI  and  Maeda,  1980).  This  technique  allows  one  to  empirically 
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determine  the  most  general  linear  relationship  between  a  solar  wind  input  function 
and  a  geomagnetic  disturbance  output  function  taking  into  account  time  delays  and 
frequency  response.  Such  a  functional  relationship  provides  a  more  detailed  set  of 
constraints  for  models  which  must  account  for  the  physical  mechanisms  associated 
with  the  coupling  of  solar  wind  energy  to  magnet ospheric  current  systems. 

The  technique  of  linear  prediction  filtering  is  based  on  least  squares  criteria  and 
leads  to  filters  (functions  which  transform  the  input  into  the  output)  which  are  linear. 
The  method  is  particularly  noteworthy  because  of  the  quality  of  the  results  obtained 
and  the  simplicity  of  the  concepts  involved.  The  theory  on  which  linear  prediction 
analysis  is  based  was  originally  developed  by  the  mathematician  Norbert  Wiener  and 
is  put  forth  in  his  book  on  linear  least  squares  prediction  and  signal  enhancement 
(Wiener,  1942).  The  technique  is  now  commonly  refered  to  as  Wiener  Filtering. 
While  the  original  work  by  Wiener  was  formulated  in  terms  of  continuous  time  series, 
the  important  adaptation  to  discrete  time  series  was  developed  by  Levinson  (1942) 
and  appears  as  an  appendix  to  Wiener’s  book.  Perhaps  one  of  the  most  useful  books 
to  deal  with  the  application  of  Wiener  filtering  to  geophysical  problems  is  by 
Robinson  (1967).  This  book  provides  a  clear  description  of  the  technique  plus 
Fortran  computer  code  for  the  subroutines  required  to  perform  the  analysis.  The 
more  recent  book  by  Robinson  and  Treitel  (1980)  also  provides  a  good  reference. 
In  the  following  sections  of  this  review,  I  will  provide  a  mathematical  introduction  to 
the  technique  followed  by  a  review  of  the  applications  which  have  used  Wiener 
filtering  in  studies  of  solar  wind-magnetosphere  coupling. 

2.  Principles  of  Digital  Wiener  Filtering 

The  basic  model  for  the  filtering  process  considered  here  consists  of  a  system 
having  an  input  signal,  a  desired  or  measured  output  signal  and  a  computed  output 
signal.  This  is  shown  schematically  in  Fig.  1  where  /(/)  is  the  input  time  series,  0(0  is 
the  desired  output  time  series  and  0(0  is  the  computed  output.  The  input  and  the 
computed  output  are  related  by  a  filter / which  must  be  determined  from  the  input  and 
the  desired  output  time  series.  A  basic  assumption  of  the  analysis  is  that  the  physical 
systems  characterized  by  /(f)  and  0(0  can  be  related  by  a  linear  time  invariant  filter /. 
This  can  be  described  mathematically  by  the  equation 
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0(0  =  j  ^fr)/(r  -  r)dr.  (1) 

For  the  case  of  discrete  time  series  data.  Equation  (1)  may  be  written  as 

00  -J  00 

0,=  I  =  I  +  I  fxlt-t.  (2) 

r=-oo  r=-oo  r=0 

Here  /,  is  the  prediction  filter  or  impulse  response  function,  t  denotes  time  and  r 
denotes  time  lag.  Note  that  Equation  (2)  can  be  broken  into  two  summations,  one 
with  r<0  and  the  second  with  r>0.  Negative  values  of  r  indicate  times  in  the  future. 
Thus,  applying  limits  to  the  summation  can  force  assumptions  about  causality.  If 
r>0,  the  filter  will  be  purely  causal  and  the  output  will  strictly  follow  the  input. 

Given  finite  and  noisy  measurements  of  /,  and  O,  our  task  is  to  solve  for  the 
predictive  filter /.  To  do  this,  we  must  assume  that  both  the  input  and  output  time 
series  are  stationary.  That  is,  their  statistical  properties  do  not  change  in  time.  We  also 
must  assume  that  errors  in  the  two  time  series  I,  and  O,  are  random  and  uncorrelated. 
And  finally  we  must  assume  that  they  are  related  by  a  linear  operator  as  in  Eq.  (2).  If 
these  assumptions  are  not  met  this  procedure  will  still  produce  the  best  linear  filter 
which  transforms  the  input  into  the  output.  This  filter,  however,  will  be  based  on  the 
finite  set  of  data  utilized  in  the  analysis  and  may  not  be  generally  valid  for  other 
subsets  of  the  infinite  data  domain  from  which  the  data  was  sampled. 

The  determination  of  the  filter  coefficients  (To,  fufi,  ...,  fm)  is  based  on  the 
least-squares  criterion:  one  minimizes  the  power  existing  in  the  difference  time  series 
obtained  by  subtracting  the  actual  or  measured  output  O,  from  the  computed  output 

6..  This  will  be  illustrated  through  the  example  of  a  single  channel  filter. 

2. 1.  Single  channel  filler 

The  single  channel  filter  has  one  input  and  one  output.  In  this  example  let 
X,= input, 

Ei=desired  output, 
f,= computed  output. 

Assuming  a  causal  relationship,  the  convolution  Equation  (2)  will  become 

m 

f,  =  foX,  +/,  fmX,  m  =  Zf,X,-i.  (3) 

t=0 

The  error  at  point  t  is  given  by  the  difference  between  the  measured  output  and 
the  computed  output.  This  is  expressed  by  e,  in  the  following  equation 

m 

e,~  Y,  —  Y,  —  I/A',  ,  (4) 

<-o 

The  mean-square-error  or  error  power  is  then  given  by 
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/  =  EWA 
=  E{[Y, 


(foX,  +fiX,  ,  +  ...  +  fmX,-m)]2} 


(5) 


where  £  is  the  Ensemble  Average  operator.  This  quantity  is  minimized  by  setting  its 
partial  derivatives  with  respect  to  the  filter  coefficients  f  equal  to  zero.  That  is 
dljdf=0.  Doing  this  produces  the  set  of  equations  called  the  normal  equations 


—  =  £{[T,  -  (foX,  +  +  ...  +  fmX,-m)\XA  =  0 

o/o 

~r  =  £{[ Y,  -  (foXi  +  +  ...  +  /mA',-m)]A',-i|  -  0 

Oj\ 


dl 

—  =  £{[T,  -  {foX,  +  y, A',- ,  +  ...  +  fmX,-m)]X,-„)  =  0. 

Ojrt t 


(6) 


where  both  sides  of  the  equations  have  been  divided  through  by  -2. 

These  equations  can  be  written  in  terms  of  the  autocorrelation  coefficients 
R,,=E{X,-,X,-j }  and  the  cross  correlation  coefficients  G,=£{  YtX,-,}  becoming 


foRdO  +/|£l0  +  fl  R20  +  ...  +  fmRmO  ~  Go 
foRol  +  f\Rl\  +  flRlI  +  ...  +  fmRm\  —  Gl 

foRom  +  f\R\m  +  flRlm  +  ...  +  fmRmm  =  O'*.  (7) 


Given  the  covariances  R,,  and  g,  we  can  solve  these  m+  1  equations  for  the  m+  I 
unknown  filter  coefficients  f. 

Since  R,,=  R,,  we  can  simplify  the  notation  by  letting  R,  equal  the  autocorrelation 
coefficient  at  lag  r.  The  normal  equations  form  a  set  of  m+  I  linear  simultaneous 
equations  which  can  be  written  in  matrix  form  as 


Ro  £,  Ri  ...  Rm  ' 

R\  Ro  R\  ...  Rm-\ 

I  fo 

/l 

h 

_ 

f  Go  \ 
G, 

Gi 

^  Rm  Rm- 1  Rm- 2  >..  Ro  ^ 

■  < 

\  1 

These  equations  could  be  solved  by  normal  techniques,  however,  because  of  the 
symmetries  in  the  matrices  a  much  more  efficient  method  is  possible.  The  matrix  R 
has  the  special  form  of  a  Tdplitz  matrix  (i.e.,  a  matrix  with  equal  submatrices  on  any 
diagonal).  Thus  the  matrix  involves  only  m+ 1  distinct  elements  namely  £0,  £i,  £2, 
Rm.  An  efficient  computational  solution  to  these  equations  is  the  Tdplitz  recursive 
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method.  This  method  is  essentially  the  same  as  the  Levinson  algorithm  (Levinson, 
1949).  A  description  of  the  recursion  will  not  be  given  here  but  can  be  found  in 
Robinson,  (1967)  and  Robinson  and  Treitel,  (1980).  Using  the  recursion 
technique,  one  usually  computes  filters  of  successively  increasing  length.  The  value  of 
the  minimum  mean-square-error  /m,„  decreases  with  each  recursion.  The  recursion 
from  n  to  n+1  may  be  stopped  whenever  a  certain  preset  criterion  is  met.  For 
example: 

1.  The  filter  reaches  a  maximum  preset  length. 

2.  The  normalized  mean-square-error  reaches  some  minimum  preset  value. 

3.  The  normalized  mean-square-error  sequence  levels  out  and  shows  little 
likelihood  of  further  appreciable  decrease. 

When  the  analysis  is  completed,  we  have  a  computed  filter  which  transforms  the 
input  into  the  output  time  series.  In  addition,  we  have  a  measure  of  the  filter’s 
efficiency  at  converting  the  input  to  the  output.  This  efficiency  measure  is  based  on  the 
power  in  the  residual  or  error  time  series  and  is  given  by  the  value  of  /=/m, „  from  Eq. 
(5). 

2.2.  Multichannel  filters 

Many  systems  have  multiple  inputs  and  outputs.  Such  systems  are  called 
multichannel  systems.  A  multichannel  system  involves  many  time  series  which  are 
related  The  theory  of  multichannel  filters  is  conceptually  the  same  as  the  theory  of  the 
single  channel  filter  and  mathematically  it  represents  only  an  extension  from  ordinary 
scalar  algebra  to  matrix  algebra.  That  is,  multichannel  theory  can  be  obtained  from 
single  channel  theory  by  the  appropriate  substitution  of  matrices  for  scalars.  If  the 
system  has  the  same  number  of  input  and  output  channels,  the  equation  will  contain 
square  matrices  and  vectors.  If  the  number  of  input  channels  and  output  channels  is 
different,  then  the  matrices  will  be  rectangular. 

For  example,  consider  a  system  having  two  input  channels  and  two  output 
channels  where 


and 


=  input. 


Y,  = 


—  output. 


The  basic  convolution  formula  is. 


Y,  ~  fnX,  +  fX,  ,  +  ...+  (9) 

where  t  is  again  an  integer  denoting  discrete  time.  Y  and  X  are  2X  1  column  vectors 
and  the  coefficients  f  are  given  by  2x2  matrices.  Expanding  Equation  (9)  yields. 
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In  general,  given  JV  input  channels  and  M  output  channels,  one  has: 


(MO) 

MO)  \  I 
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X,  =  XU)  ~ 


Y,  =  Y(t)  = 


Yt  =  K(r)  = 
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XiU) 
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YtU)  \ 
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Input, 


=  Desired  Output, 


=  Computed  Output, 


(10) 


and  a  set  of  A/X7V  matrix  valued  filter  coefficients 


f  =  As)  = 


Ms)  Ms) 

fu(s)  fn(s) 


Ms) 

fl  rAs) 


fm{s)  /m7(s)  •••  /mnO)  j 


The  error  e,  between  the  desired  output  and  the  actual  output  is  then  a  M X 1 
vector  valued  time  series 


(11) 


Ur)  -  >',(/)  \ 

e,U)\ 

-  Y,  = 

U0  ~  YAi) 

= 

e2(t) 

Mr)  -  >«(/)  I 

.  euU) 

The  mean-square-error  matrix  is  defined  as 
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E{e\(f)\  £{<?,(/)€2(/)|  ...  E{e,(t)e^t)\  \ 

E{e,eT,}=  ...  (12) 

E{eM(t)e\(t)\  E{eu(t)e2(i)  ...  £{«(f)}  /• 

The  mean-square-error  or  error  power  is  given  by  the  trace  of  the  mean-square-error 
matrix  which  is  defined  to  be  the  sum  of  the  diagonal  terms, 

m 

/  =  uE\e,e\\  =  I£(eJ,(f))|.  (13) 

Again  setting  the  partial  derivative  of  /  with  respect  to  each  filter  coefficient  equal  to 
zero  produces  a  set  of  normal  equations  given  by: 

AO)  1)  +  ..  +fm<t>iA~m)  =  <M  0) 

M'A  I)  +  /.<M  0)  +  ...  +fm<t>,A\-m)  =  <M0 

M<Am)  T/i^m-l)  +  ...  +/m<M0)  =  <M m)  (14) 

where  d>t;r(r)  are  the  autocorrelation  of  the  input  X  at  lag  r  and  tj>yA t)  are  the  cross 
correlation  of  the  desired  output  Y  with  the  input  X.  This  set  of  linear  simultaneous 
equations  can  also  be  solved  by  the  Toplitz  recursion  method  discussed  above.  The 
reader  is  directed  to  texts  by  Robinson  for  further  information  about  the  calculation 
of  the  multichannel  autocorrelation  and  cross  correlation  coefficients  and  the  Tftplitz 
recursive  solution. 

3.  Application  to  Studies  of  Solar  Wind-Magnetosphere  Coupling 

Wiener  filtering  has  been  applied  to  the  problem  of  solar  wind-magnetosphere 
coupling  using  various  parameters  measured  in  the  solar  wind  and  indices  of 
geomagnetic  activity.  In  these  studies,  solar  wind  parameters  represent  the  input  and 
geomagnetic  activity  the  output.  The  computed  filters  or  impulse  response  functions 
characterize  the  processes  which  link  solar  wind  energy  to  the  currents  monitored  by 
the  various  geomagnetic  indices  used  in  the  analyses. 

The  technique  of  Wiener  filtering  was  introduced  to  the  solar  wind-magneto¬ 
sphere  coupling  problem  by  iYEMORl  et  al.  (1979).  In  this  study  the  authors  used  300 
days  of  hourly  average  d ata  from  1 967  through  1 969.  Filters  relating  the  solar  wind 
quantities  VBS  and  A  to  auroral  magnetic  indices  AL,AU,AE  and  to  the  low-latitude 
ring  current  magnetic  index  £)„  were  computed.  The  D»  filter  develops  in  1  hour  and 
decays  with  a  time  constant  of  about  8  hours  which  is  consistent  with  the  previous 
empirical  model  developed  by  Burton  et  al  (1975).  The  AE  and  AL  filters  are  a 
narrow  pulse  peaking  at  1  hour  lag.  The  A  (/filter  is  a  broader  peak  at  I  hour  lag.  Since 
the  AU  filter  is  significantly  different  than  the  AL  filter  in  shape,  IYEMORl  et  al, 
suggest  that  the  generation  for  the  currents  monitored  by  A  U  may  be  different  than 
those  for  AL  or  else  the  generation  may  be  complex,  consisting  of  multiple 


mechanisms.  The  prediction  efficiency  was  fairly  high  for  all  of  the  predictions. 

The  Wiener  Filtering  technique  was  extended  to  studies  of  high  time  resolution 
data  by  a  number  of  subsequent  researchers  Cl  AUER  et  al.,  1981, 1983;  McPHERRON 
etal.,  1984;  BarGATZE  era/.,  1985;  and  Akasofu  et  al.,  1985).  Figure  2  from  CLAUER 
et  al.  ( 1 98 1 )  shows  high  time  resolution  Wiener  filters  which  relate  3  forms  of  the  solar 
wind  input  function  to  the  A  L  index.  The  three  input  functions  have  a  KZKsin4/#  2), 
F&,  an  dependence  respectively  where 

tan"1  -^77,  for  B:  >  0; 

l&l 

180°  —  tan'1  forB.- <0. 

Iftl 

and 

(0,  for  B:  >  0 
‘  ~  for  B:<  0. 

Two  time  intervals  were  analyzed.  The  filters  at  the  top  of  Fig.  2  were  created 
from  an  interval  of  moderate  level  goemagentic  activity  while  the  bottom  panel  filters 
were  created  from  an  interval  of  strong  activity.  The  horizontal  axis  of  the  plots  is  time 
lag.  A  smooth  curve  is  drawn  through  the  filters  and  is  shown  only  as  an  aid  to  the 
reader.  The  filters  are  very  spiky.  This  is  due  to  uncertainties  which  result  from 
correlations  based  on  a  small  number  of  active  periods  in  the  10  and  12  day  intervals 
used  in  the  analysis.  Less  noisy  filters  result  from  analysis  on  long  intervals  which 
include  many  repeated  correlated  variations  between  the  input  and  output  time  series. 
For  example,  compare  the  K&— /ILfilters  in  Fig.  2  with  the  VB,—  AL  filter  shown  in 
Fig.  5  which  was  created  using  56  days  of  data. 

The  dominant  feature  of  the  filters  in  Fig.  2  is  a  delayed  negative  pulse  indicative 
of  a  low  pass  filter  with  a  time  delay.  It  is  interesting  to  see  that  the  filters  for  the  strong 
activity  and  moderate  activity  intervals  show  different  time  delays  suggesting  some 
form  of  nonlinearity  in  the  KB,— AL  relationship. 

The  variation  in  the  Kft—  AL  impulse  response  as  a  function  of  geomagnetic 
activity  level  has  been  examined  in  greater  detail  by  Bargatze  et  al.  (1985).  The  data 
set  examined  contains  73  days  of  2.5  minute  resolution  data.  An  activity  level  based 
on  the  AL  index  was  assigned  to  each  of  the  34  intervals  comprising  the  data  set.  These 
intervals  were  sorted  by  activity  level  to  produce  a  time  series  of  increasing  level  of 
geomagnetic  activity.  A  set  of  running  filters  was  determined  over  the  entire  data  set 
by  first  using  5  consecutive  intervals  to  compute  a  Wiener  filter  and  then  stepping 
ahead  one  interval  at  a  time.  Figure  3  shows  the  set  of  30  VBS— AL  filters  thus 
obtained  with  the  weak  activity  filters  at  the  bottom  ranging  to  the  strong  activity 
filters  at  the  top. 

Two  peaks  can  be  identified  in  these  filters,  the  first  at  20  minutes  lag  and  the 
second  at  60  minutes  lag.  The  relative  importance  of  the  peaks  varies  as  a  function  of 
activity  level,  with  the  second  peak  dominant  at  middle  activity  levels,  the  first  peak 
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TIME  LAG  (hours) 

Fig.  3.  Slack  plot  of  Wiener  Filters  for  varying  levels  of  geomagnetic  activity  ranging  from  quiet  at  the 
bottom  to  severe  at  the  top  (Bargatze  ei  at..  1985) 


dominant  at  strong  activity  levels  and  both  peaks  about  equal  at  weak  activity  levels. 

It  is  suggested  that  the  two  peaks  correspond  to  separate  response  modes  of  the 
magnetosphere  and  that  the  relative  dominance  of  these  modes  varies  with  the  level  of 
activity.  The  physical  interpretation  given  to  the  two  peaks  regards  the  first  peak  as 
corresponding  to  auroral  activity  which  is  directly  driven  by  the  solar  wind- 
magnetospheric  dynamo  while  the  second  peak  corresponds  to  auroral  activity  which 
results  from  the  release  of  stored  energy  from  the  magnetotail 

While  auroral  currents  have  obtained  the  greatest  attention  in  linear  prediction 
investigations,  four  studies  have  presented  analysis  which  relate  solar  wind  para¬ 
meters  to  the  ring  current  index  D».  lYEMORi  et  al  (1979);  lYEMORI  and  Maeda 
(1980);  McPherron  el  al.  (1984);  and  indirectly  AKASOFU  et  al.  (1985).  The  results 
of  these  studies  indicate  that  D„  is  the  best  determined  geomagnetic  index  using  solar 
wind  information  alone.  Usii.g  5-minute  data  and  a  multi-channel  filter  with  solar 
wind  dynamic  pressure  and  electric  field  inputs,  McPHERRON  et  al.,  (1984) 
reproduced  70%  of  the  variance  in  the  DS1  index. 

Most  of  the  studies  using  linear  prediction  filtering  have  concluded  that  only  a 
portion  of  the  output  index  is  predictable  from  solar  wind  information  alone  lYEMORi 
et  al.,  1979,  1980;  Clauer  et  al.,  1981,  1983;  McPherron  et  al.  1984;  and 
Bargatze  et  al.,  1985).  This  is  illustrated  for  the  A  L  index  in  Fig.  4  which  shows  a 
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COMPARISON  OF  PREDICTED  WITH  OBSERVED  AL  INDEX 


Fig.  4  Comparison  of  predicted  and  observed  41  index  for  periods  of  moderate  (top)  and  strong 
(bottom)  geomagnetic  activity  Each  frame  'hows  a  plot  of  the  solar  wind  input  function  f 'B,  at  the  top 
with  the  predicted  (dark  line)  and  ooscrved  (light  line)  A  L  index  underneath,  (figure  courtesy  of  R  L 
McPherron) 


comparison  between  the  measured  AL  index  and  the  computed  AL  index  The  top 
panel  shows  an  interval  of  mode’-'te  level  activity  and  the  bottom  panel  shows  an 
interval  of  strong  activity.  In  each  case  the  computed  AL  looks  like  a  smoothed  or 
low-pass  filtered  version  of  the  observed  index,  it  is  suggested  mat  unpredicted 
portions  of  the  index  can  be  attributed  to  noise  in  the  measurements  tT  the  proper 
input  parameters  and  in  the  indices  which  characterize  the  magnetosphei ic  current 
systems,  and  also  to  real  variations  in  the  currents  which  are  not  consistently  related 
to  the  solai  wind  driving  function  That  is,  currents  which  a  e  controlled  by  internal 
magnetospheric  processes  There  are  several  examples  of  large,  sudden,  negative 
excursions  of  the  AL  index  which  arc  net  predeted  and  which  ar  ■  thought  to  result 
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from  intensifications  of  the  westward  electrojet  caused  by  substorm  expansions 
initiated  within  the  magnetotail  by  internal  processes. 

The  analysis  presented  by  Claler  ei  at.  ( 1983)  provides  a  fairly  strong  case  for 
this  point  of  view.  Fifty-six  days  of  10-minute  resolution  data  were  used  in  the 
analysis.  Figure  5  shows  the  filters  computed  in  that  study.  Note  that  both  future  and 
past  time  lags  were  used  in  this  analysis.  The  solar  wind  input  function  used  in  the 
analysis  was  VB ,  and  the  magnetic  indices  examined  were  the  auroral  A  L  index  and 
the  low-latitude  ASYM  index  defined  to  be  the  difference  between  the  dawn  and  dusk 
perturbations  in  the  X  component  of  the  geomagnetic  field.  ASYM\s  thus  an  index  of 
the  partial  ring  current  magnetic  perturbation. 

The  FSs— ASYM  and  VB^—AL  filters  in  Fig.  5  are  essentially  zero  prior  to  0 
time  lag,  thus  they  are  purely  causal.  Both  are  delayed  pulses  and  are  essentially  low 


Linear  Prediction  Filters  Relating 
Solar  Wind  Electric  Field  To 
AL  and  ASV  Indicies  and  to  Each  Other 

(56  days  in  1967  -  1968) 


Time  Lag  (hours) 


Fig  5.  Three  Wiener  Filters.  AL  —  ASYM,  I  fl, -A.' V  YM  and  IB,—  AL.  plotted  together  for  compar¬ 
ison  (Cl  xi  er.  el  al.  1983) 
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pass  filters  with  a  time  delay.  The  A  L—ASYM  filter  is  also  show  n  in  Fig.  5.  This  filter 
was  created  to  test  the  notion  that  the  partial  ring  current  is  the  result  of  individual 
substorm  expansions  injecting  bursts  of  energetic  ions  into  the  ring  current  region. 
These  ions  would  drift  westward  to  produce  a  partial  ring  current  (BogoiT  and 
Mozer,  1973).  We  find,  however,  that  the  filter  is  nonzero  at  negative  lags  indicating 
that  development  of  ASYM  precedes  AL.  This  feature  of  the  Wiener  filter  is 
consistent  with  the  analysis  of  several  individual  events  which  showed  the  develop¬ 
ment  of  the  partial  ring  current  preceding  the  substorm  expansion  onset  determined 
from  both  the  AL  index  and  the  midlatitude  positive  bay  signature  (Claier  and 
McPherron,  1980). 

The  measures  of  the  prediction  efficiency  provide  some  interesting  information. 
About  459c  of  the  variance  of  the  A  S  YM  index  is  predictable  from  VB„  while  56%  is 
predictable  from  AL.  Also,  41%  of  the  variance  of  A L  is  predictable  using  VB,.  This 
suggests  that  the  A  L  index  contains  information  that  is  correlated  with  ASYM  but  is 
unrelated  to  the  solar  wind.  Cl  At  ER  et  al.  (1983)  hypothesize  that  both  AL  and 
ASYM  are  produced  by  two  uncorrelated  processes,  one  which  directly  involves  the 
solar  wind  and  another  which  is  due  to  internal  magnetospheric  processes. 

Additional  analysis  was  performed  on  the  error  (or  residual)  time  series.  It  was 
reasoned  that  the  solar  wind  related  portion  of  the  signal  is  absent  in  the  residual  time 
series.  There  should,  however,  still  remain  the  signal  due  to  internal  magnetospheric 
processes  which  are  not  directly  determined  by  the  solar  wind  input  alone.  To  see  if  the 
error  time  series  for  AL  and  ASYM  both  contained  this  internal  magnetospheric 
signal,  a  filter  relating  the  prediction  residual  of  AL  to  the  prediction  residual  of 
ASYM  was  computed.  It  was  found  that  24%  of  the  variance  of  the  ASYM  residual 
was  related  to  the  AL  residual. 

Figure  6  shows  a  schematic  model  of  the  hypothesized  system.  The  solar  wind 
input  VB,  is  transformed  to  the  output  A  L,  by  filter  g] .  A  second  component  A  L,  is 
assumed  to  be  produced  by  the  internal  magnetospheric  process  through  filterg,  and 
is  added  to  A  L,  to  produce  A  L.  The  A  L  index  in  turn  is  transformed  by  filter#;  into 
the  ASYM  index.  Therefore,  in  this  model  the  solar  wind  input  VB,  is  related  to 
A  S  YM  by  a  cascade  of  two  filters  gi  and#;.  Provided  that  both  systems  are  linear,  the 
cascade  can  be  represented  by  a  single  filter  g„  which  is  the  convolution  of  g i  and  g;. 
Thus, 


g'hm  S.  g\.ng2,lm  m ■ 
n  ') 

Figure  5  shows  the  empirically  determined  filters  g„  and  g;  The  output  of  filter  g,. 
AL,  is  not  directly  observed.  By  hypothesis,  however,  AL ,  is  uncorrelated  with  the 
solar  wind  so  that  gi  is  the  filter  determined  from  I  B,  and  AL.  As  a  check  of  this 
hypothesis  the  I B,  —  AL  filler  was  convolved  with  the  AL  —  ASYM  filter  for 
comparison  with  the  VB,—  ASYM  filter.  The  comparison  is  shown  in  fig.  7.  A  dashed 
line  shows  the  empirical  l  B,  —  AS  YM  Idler  and  a  solid  line  show  s  the  filter  produced 
by  the  convolution.  The  agreement  between  the  two  curves  is  quite  good,  thus 
supporting  the  hypothesis  that  Al.  is  composed  of  both  a  solar  wind  driven 
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A  possible  Unear  System  Relating  vbs,  AL.  ASYM 


I 


•  ASYM 


Fig.  6.  Schematic  diagram  of  a  model  in  which  the  solar  wind  input  function  FB,  is  transformed  into  the 
ASYM  index  by  filter  gu.  The  AL  index  is  transformed  into  ASYM  by  filter  g;.  The  AL  index, 
however,  is  composed  of  two  components  One  component  AL.  is  related  to  the  solar  wind  FB, 
through  filter  gi  while  the  other  component  AL,  is  related  to  another  process  X  by  filter  g,.  AL  is  the 
sum  of  these  two  signals  (Cl  At  t R  ei  a/.,  1983). 


I 
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COMPAR I S I  ON  OF  VBS  TO  ASYM  FILTERS 


TIME  (HOURS) 

Fig.  7.  Comparison  of  the  FB,  .45  YM  filter  (go  in  Fig  5)  with  the  convolution  of  the  FB,  -AL  and 
AL— ASYM  filters  (gi  and  g:  in  Fig.  5)  (Cl  At  fr  el  a! .  1983). 
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component  and  another  component  uncorrelated  with  the  solar  wind. 

A  somewhat  different  conclusion,  however,  was  reached  by  Akasofu  et  al. 
(1985)  in  an  analysis  which  computed  Wiener  filters  relating  the  c  parameter  with  the 
so-called  total  magnetospheric  energy  dissipation  rate  parameter,  U r.  This  investiga¬ 
tion  used  5-minute  resolution  data  during  magnetic  storm  intervals.  The  solar  wind 
input  parameter  £<*  Vlfsin*  (8i2)  was  measured  by  the  I  SEE-3  satellite  located  240  R< 
upstream  from  the  subsolar  magnetopause.  The  output  parameter  Ut  is  given  by: 

L'i(erg  sec)  =  4- 10'"  (  +  —  I  -i-  3- 10 'iAE 

\  di  Tr  / 

where  ZX,  is  a  solar  wind  pressure  adjusted  D„  index  and  tr  denotes  a  ring  current 
decay  time  scale. 

Figure  8  show  plots  from  3  of  the  6  storms.  Each  column  of  plots  shows  data  from 
one  storm.  The  top  panel  shows  the  input  c  parameter,  the  next  panel  shows  the 


Fig  8.  Fxamplo  at  linear  predation  analysis  for  3  geomagnetic  storms  on  September  2?  1978  F.omihc 
top  are-,hown  the  input  parameter* ,  the  output  I  unction  I  /.  the  prediction  of  IT.  and  the  prediction 
filter  I  Ak  tsnu  el  ul .  19X5) 
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output  Uj  parameter,  the  next  panel  shows  the  predicted  0 t  parameter  and  the 
bottom  panel  shows  the  computed  Wiener  filter.  Figuer  9  shows  a  superposition  of  all 
six  of  the  computed  Wiener  filters. 

The  authors  note  that  for  the  3  largest  storms  studied  the  predicted  0 r  is  a  good 
fit  to  the  measured  U t  while  the  smaller  3  storms  were  not  well  fit.  The  filters  shown  in 
Figs.  8  and  9  are  essentially  S  functions  at  /  =  0.  These  results,  the  zero  time  lag  and  the 
good  prediction,  are  interpreted  to  indicate  that  for  periods  of  strong  activity 
magnetospheric  energy  dissipation  is  directly  driven  with  little  or  no  time  lag. 

The  zero  time  lag  in  the  Akasofu  el  al.  filters,  however,  is  troubling  and  seems 
physically  implausible,  particularly  when  the  uncertainty  of  the  transit  time  between 
ISEE-3  and  the  magnetopause  is  greater  than  the  accuracy  of  the  5-minute  data 
resolution  (RUSSELL  et  al,  1980).  Further,  spikes  at  the  ends  of  filters  are  often 
characteristic  of  errors  in  removing  the  mean  or  trends  from  time  series.  Alternatively, 
it  is  likely  that  the  methods  used  to  compute  U r  are  cause  of  this  artifact.  The 
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Fig.  9.  Superposition  (top)  and  average  (bottom)  of  6  impulse  response  functions  computed  for  6 
geomagnetic  storms  (Akasofi'  el  al..  1985). 
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derivation  of  U^  is  based  on  several  solar  wind  parameters.  The  value  of  At  is  adjusted 
for  the  solar  wind  pressure  and  the  ring  current  decay  “constant”  rR  is  a  function  of  the 
input  function  c.  The  output  parameter  therefore  contains  terms  from  the  measured 
input  thereby  confusing  the  correlation  and  producing  a  zero-lag  artifact.  One  must 
also  note  that  the  application  of  Wiener  filtering  to  single  event  studies  may  be 
inappropriate.  The  spikey  character  of  the  filters  obtained  is  an  indication  of  this.  The 
ability  to  compute  physically  meaningful  filters  relies  on  the  least  squares  aspects  of 
the  analysis  and  requires  many  correlated  variations  in  the  input  and  output  time 
series. 

The  pressure  effect  could  more  properly  be  taken  care  of  by  using  a  multichannel 
analysis  which  inputs  both  solar  wind  dynamic  pressure  and  c  to  an  output  composed 
only  of  geomagnetic  measurements.  To  do  this,  rR  must  be  obtained  independently 
from  the  input  time  series  c.  For  example,  FELDSTEIN  etal.  (1984)  estimate  values  for 
rR  based  on  the  magnitude  of  A,.  Such  a  multichannel  technique  would  also  provide  a 
better  method  of  determining  the  transit  time  from  1SEE-3  to  the  magnetopause  since 
thepF2  —  A,  filter  was  shown  to  be  a  narrow  pulse  at  a  delay  time  consistent  with  the 
expected  propagation  time  between  the  satellite  and  magnetopause  (McPHERRON  et 
al.,  1984). 

4.  Summary  and  Conclusions 

The  technique  of  linear  prediction  filtering  has  provided  another  powerful 
statistical  tool  for  the  data  analyst  to  use  in  an  attempt  to  resolve  some  of  the  complex 
relationships  which  exist  between  various  geomagnetic  current  systems  and  solar 
wind  input  energy.  It  appears  clear  from  all  of  the  Wiener  filtering  studies  cited  above 
that  a  substantial  portion  of  the  geomagnetic  disturbances  observed  on  the  earth  can 
be  directly  related  to  solar  wind  measures  by  linear  impulse  response  functions.  This 
suggests  that  the  responsible  current  systems  are  very  directly  coupled  to  (or  driven 
by)  the  solar  wind-magnetosphere  MHD  dynamo.  In  addition,  however,  there 
appears  to  be  a  substantial  portion  of  the  A  L  index  which  results  from  currents  not 
strictly  determined  by  the  solar  wind  input. 

Investigations  have  used  Wiener  analysis  to  compute  the  linear  prediction  filters 
which  relate  solar  wind  input  parameters,  primarily  c  and  K&,  with  the  current 
systems  manifested  in  A  L,  A  U,  A,  and  A  S  YM.  In  general,  these  filters  have  low  pass 
characteristics  with  linear  phase  and  time  delays.  The  different  filters,  however, 
quantitatively  show  that  different  indices  respond  differently  to  the  solar  wind  electric 
field  and  hence,  that  different  physical  process  produce  the  electric  currents  which 
influence  each  index.  It  also  seems  clear  that  the  assumptions  of  linear  and  time 
invariant  relationships  between  the  input  and  output  are  only  approximately  valid. 

The  A,  index  appears  to  be  the  most  predictable.  Solar  wind  inputs  can 
reproduce  over  70%  of  the  variance  in  the  A,  index  (lYEMORl  et  al.,  1979; 
McPHERRON  et  al.,  1984).  The  A  l.  and  ASYM  indices  can  be  predicted  relatively  well 
with  about  40%  of  the  variance  of  each  predicted  by  the  solar  wind  electric  field  or  r.. 
This  seems  impressive  in  view  of  the  errors  present  in  the  measurements  of  both  the 
input  and  output  time  series.  The  A  V  index  is  only  poorly  predicted  and  this  is  a  major 
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concern  to  this  author.  This  may  indicate  either  a  major  misunderstanding  of  the 
sources  controlling  the  eastward  electrojet  or  serious  problems  with  the  interpretation 
that  A  U  in  fact  characterizes  the  eastward  electrojet. 

The  low  pass  characteristics  of  the  VB.^AL  filters  indicate  that  the  magneto¬ 
sphere  attenuates  solar  wind  electric  field  frequencies  above  10  4  hertz  in  the 
generation  of  the  current  systems  linked  with  the  westward  auroral  electrojet.  Also, 
careful  examination  of  the  residual  or  error  time  series  suggests  that  there  are 
variations  in  the  output  indices  which  are  produced  by  processes  not  directly  related 
to  the  solar  wind  input  function.  It  has  been  suggested,  for  example,  that  these 
separate  sources  of  current  in  the  westward  electrojet  are  controlled  by  processes  in 
the  magnetotail  which  are  associated  with  the  release  of  stored  energy. 

The  results  of  these  analyses  provide  an  empirical  description  of  the  solar 
wind-magnetospheric  system.  This,  in  turn,  provides  insight  to  the  physical  processes 
which  operate  in  the  system.  Most  of  the  studies  reported  above  have  utilized  Wiener 
filtering  in  this  way.  The  technique,  however,  may  eventually  be  employed  in  the  more 
pragmatic  way  of  providing  the  ability  to  forecast  short  term  geomagnetic  activity  if 
the  proper  inputs  were  provided  sufficiently  far  upstream  in  the  solar  wind.  A  satellite 
at  the  Lagrangian  point  240  R,  upstream,  for  example,  might  provide  such  short  term 
input.  If  proper  filters  could  be  produced  to  relate  observations  at  this  location  with 
geomagnetic  activity  indices,  predictions  could  easily  be  computed.  Based  on  present 
results,  it  seems  reasonable  to  expect  that  such  filters  could  be  obtained  and 
reasonably  accurate  predictions  could  be  produced.  Such  predictions  could  have 
great  benefit  to  a  variety  of  commercial  and  scientific  activities  which  are  influenced 
by  the  level  of  geomagnetic  disturbance. 
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IONOSPHERIC  TRAVELING  CONVECTION  VORTICES  OBSERVED  NEAR  THE  POLAR 
CLEFT:  A  TRIGGERED  RESPONSE  TO  Sl'DDEN  CHANGES  IN  THE  SOLAR  WIND 

E.  Friis-Christenscn1 ,  M.  A  McHenry2.  C.  R.  Clauer2.  S.  Vennerstrom1 


Abstract.  Analysis  of  20-second  resolution  magnetometer  data 
from  an  array  of  temporary  stations  operated  around  Sondre 
Stramfjord,  Greenland  during  the  summer  of  1986  shows  the 
signatures  of  localized  ionospheric  traveling  convection  vor¬ 
tices.  An  example  of  an  isolated  event  of  this  kind  observed 
near  08  local  time  is  presented  in  detail.  This  event  consists  of 
a  twin  vortex  pattern  of  convection  consistent  with  the  pres¬ 
ence  of  two  field-aligned  current  filaments  separated  by  about 
600  km  in  the  east-west  direction.  This  system  of  currents  is 
observed  to  move  westward  (tailward)  past  the  array  of  sta¬ 
tions  at  about  4  km /sec.  The  event  is  associated  with  relative 
quiet  time  ionospheric  convection  and  occurs  during  an  interval 
of  northward  IMF.  It  is,  however,  associated  with  a  large  fluc¬ 
tuation  in  both  the  Z  and  Y  components  of  the  IMF  and  with 
a  large  sudden  decrease  in  the  solar  wind  number  density.  The 
propagation  of  the  system  is  inconsistent  with  existing  models 
of  FTE  current  systems,  but  nevertheless  appears  to  be  related 
to  a  readjustment  of  the  magnetopause  boundary  to  a  sudden 
change  in  the  solar  wind  dynamic  pressure  and/or  to  a  change 
in  reconnection  brought  about  bv  a  sudden  reorientation  of  the 
IMF. 

Introduction 

The  immediate  response  of  the  large  scale  polar  magneto- 
spheric  plasma  convection  to  changes  in  the  solar  wind  electric 
field  apparent  at  the  dayside  magnetopause  is  well  established 
through  a  number  of  studies  [Clauer  et  al..  1984;  Rishbeth 
et  al.,  1985;  Friis-Christensen.  1986;  Clauer  and  Banks.  1986; 
Clauer  and  FYiis-Christensen,  1988].  The  most  direct  response 
to  changes  in  the  solar  wind  is  seen  at  high  latitudes  around 
magnetic  local  noon  in  the  vicinity  of  the  polar  cleft.  This  re¬ 
gion  is  generally  thought  to  be  near  the  ionospheric  projection 
of  the  boundary  between  open  and  closed  magnetic  field  lines 

This  is  also  a  region  where  irregular  magnetic  pulsations 
often  occur.  The  intensity  of  these  pulsations,  called  ipcl-tvpe 
pulsations  (irregular  pulsations,  continuous,  long)  increases 
with  increasing  negative  values  of  IMF  fi,  [Bolshakova  et  al. . 
1975],  and  the  latitude  of  maximum  ipcl-pulsations  is  close  to 
the  equatorward  boundary  of  the  cusp  statistically  defined  by 
means  of  particle  data  [Kleymenova  et  al.  (1982). 

Kleymenova  et  al.  [1985]  divided  the  ipcl-pulsations  into 
two  classes;  “noise  pulsations”  (np)  within  the  3-8  minute 
period  range  and  “long  period”  pulsations  (lp)  with  a  quasi- 
period  from  15-20  minutes.  It  was  proposed  that  the  Ip  -  pul¬ 
sations  occur  on  open  field  lines  because  of  their  intimate  rela¬ 
tionship  to  IMF  variations,  while  the  np-pulsations  were  pro¬ 
posed  to  be  generated  in  the  boundary  layer  on  closed  field- 
lines.  It  was  found  that  np-pulsations  occur  for  both  IMF 
B,  <  0  and  B,  >  0  and  that  the  maximum  occurrence  is 
in  the  prenoon  sector  (0600  to  1100  MLT).  When  B  <  0 
the  magnitude  of  the  pulsations  is  larger  and  the  presence  of 
np-pulsations  is  often  accompanied  by  the  occurrence  of  Pc5 
pulsations  at  lower  latitudes. 
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The  latitude  of  occurrence,  relationship  to  IMF  B;.  and 
temporal  duration  of  ipcl  pulsations  seem  consistent  with  the 
properties  of  pulsations  which  are  expected  to  be  produced 
by  sporadic  reconnection  near  the  subsolar  point  of  the  mag¬ 
netopause.  Thus.  Bolshakova  and  Troitskaya.  [1982;  suggest 
that  ipcl  pulsations  may  be  related  to  “pulsed  reconnection" 
or  flux  transfer  events  (FTE’s)  described  by  Russell  and  Elphic 
[1979]. 

Two  different  models  of  FTE  field-aligned  current  systems 
have  been  proposed.  The  first  by  Saunders  et  al  [1984;  and 
Lee  [1986]  consists  of  a  central  core  of  field-aligned  current 
which  spreads  radially  in  the  ionosphere  and  closes  through 
a  local  or  distant  sheath  of  field-aligned  current.  The  second 
system  which  consists  of  two  oppositely  directed  field-aligned 
currents  on  the  flanks  of  the  FTE  flux  tube  producing  flux  tube 
motion  perpendicular  to  a  line  connecting  the  vortex  centers, 
was  proposed  by  Southwood  [1985:  1986] 

Considerable  effort  has  been  made  to  uniquely  relate  iono¬ 
spheric  signatures  to  the  existence  of  FTEs  Goertz  et  al 
[1985]  associated  STARE  ionospheric  radar  observations  of  lo¬ 
calized  regions  of  poleward  ionospheric  flow  with  the  existence 
of  an  ionospheric  signature  of  FTE's.  Lanzerotti  et  al.  [1986] 
used  magnetometer  observations  from  a  station  near  the  south¬ 
ern  cusp  (South  Pole)  to  support  the  idea  of  the  existence  of  a 
localized  Hall  current  system  consistent  with  the  FTE  signa¬ 
ture  of  a  Saunders  et  al.  -  Lee  system  [Saunders  et  al  .  1984: 
Lee,  1986]  Results  from  a  theoretical  study  of  FTE'  current 
systems  indicate  that  a  closely  spaced  array  of  high  time  res¬ 
olution  magnetometer  observations  is  needed  in  order  to  ad¬ 
equately  resolve  the  question  of  the  existence  of  ionospheric 
FTE  signatures  [McHenry  and  Clauer,  1987!. 

The  purpose  of  this  paper  is  to  report  on  the  first  results 
of  an  experiment  conducted  to  investigate,  in  detail,  the  irreg¬ 
ular  pulsations  measured  by  ground-based  magnetometers  in 
the  vicinity  of  the  dayside  cleft.  In  particular,  we  report  on  the 
observation  of  traveling  small  scale  two-cell  ionospheric  current 
systems.  These  sy  ’terns,  which  seem  to  explain  a  number  of 
the  observed  irregular  pulsations  and/or  spike-like  variations  in 
the  ground  based  magnetic  records,  tire  the  ionospheiic  mani¬ 
festations  of  a  pair  of  oppositely  directed  field-aligned  currents 
moving  rapidly  near  the  convection  reversal  boundary.  These 
events,  however,  do  not  appear  to  be  consistent  with  existing 
models  of  the  ionospheric  signatures  of  FTEs 

Observations 

To  obtain  high  tempoial  and  spatial  resolution  observa¬ 
tions  of  the  polar  cleft  ionosphere  the  Greenland  Magnetome¬ 
ter  Array  was  supplemented  with  additional  temporary  sta¬ 
tions  near  Sondre  Stromfjord.  During  the  summer  of  1986  a 
total  of  8  stations  within  the  invariant  latitude  range  from  08' 
to  76'  were  operated  with  a  time  resolution  of  20-seconds. 

A  map  showing  the  magnetometer  stations  in  Greenland 
used  in  this  study  is  shown  in  Figure  1.  The  stations  to  the 
south  of  and  including  GDH  recorded  20-second  magnetic  field 
measurements.  While  primarily  a  latitudinal  array,  some  east- 
west  separation  exists,  in  particular  between  the  SKT  and 
DYB  stations. 

Figure  2  shows  a  stacked  plot  of  H-component  magne- 
lograms  for  0700  to  1300  UT  on  June  28.  1966.  Around  1000 
to  1020  UT  (4  hours  prior  to  magnetic  local  noonj  a  spike-like 
disturbance  is  seen  at  nearly  all  the  stations,  with  a  maximum 
peak  to  peak  disturbance  of  approximately  200  nT  near  the 
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Figure  1.  Map  of  Greenland  showing  the  magnetometer  sta¬ 
tions  used  in  this  study.  Godhavn  (GDH)  and  all  stations  on 
the  west  coast  to  the  south  of  GDH  were  equipped  to  record 
20-second  resolution  magnetic  field  measurements.  The  broken 
lines  indicate  invariant  latitude  circles. 


station  ATU  at  75  3°  invariant  latitude.  In  Figure  3  the  same 
event  is  shown  in  all  three  components  with  higher  time  reso¬ 
lution  using  the  20-second  data.  From  this  figure  it  is  evident 
that  the  current  system  responsible  for  the  variations  is  small 
scale  or  localized  in  character  (large  gradients  in  the  H-  and 
Z-components)  and  that  there  seems  to  be  a  phase  shift  be- 
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Figure  2.  Stacked  plot  of  H-component  magnetograms  from 
the  Greenland  chain  of  magnetometers  for  0700  to  1300  UT 
on  June  28,  1986.  The  location  of  the  stations  are  shown  in 
Figure  1.  Magnetic  local  noon  occurs  at  1400  UT 
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Figure  3.  Stacked  plot  of  20  second  resolution  magnetograms 
from  the  Greenland  chain  for  0930  to  1030  UT  on  June  28, 
1986.  The  data  have  been  filtered  using  a  high-pass  filter 
(corner  frequency  corresponding  to  a  period  of  30  minutes). 
The  H-component  (northward)  is  shown  on  the  top  panel,  D 
(eastward)  in  the  middle-,  and  the  Z- component  (downward) 
is  given  in  the  bottom  panel. 
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tween  SKT  and  DYB.  This  latter  point  indicates  a  westward 
(anti-sunward)  movement  of  about  3  to  5  km/sec  of  the  current 
system  which  is  producing  the  pulsations. 

Assuming  that  the  perturbations  are  caused  by  a  moving 
current  structure  which  changes  very  little  in  time,  we  are  able 
to  use  the  meridian  chain  of  magnetometers  to  “scan"  the  pat¬ 
tern  as  it  moves  across  the  line  of  stations.  The  velocity  of  the 
current  structure  is  estimated  from  the  phase  shift  between 
SKT  and  DYB.  Using  an  estimated  velocity  of  4  km/sec  we 
have  plotted  in  Figure  4  the  total  horizontal  magnetic  pertur¬ 
bation  vectors  in  a  coordinate  system  moving  with  the  current 
pattern.  For  simplicity  the  perturbation  vectors  have  been 
rotated  by  90°  counterclockwise  to  be  consistent  with  the  con¬ 
vection  direction.  That  is,  we  assume  the  magnetic  signature 
to  result  from  ionospheric  Hall  currents.  Since  the  event  oc¬ 
curred  in  the  dayside  oval  at  midsummer  the  ionospheric  Hall 
conductivity  is  probably  sufficiently  large  and  uniform  to  sup¬ 
port  this  assumption. 

The  observed  convection  pattern  consists  of  a  twin-vortex 
moving  westward  (tailward)  roughly  along  a  line  connecting 
the  centers  of  the  two  vortices.  The  consonance  of  the  pattern 
derived  from  the  distribution  of  measurements  in  space  and 
time  confirms  the  velocity  determination  and  the  stability  of 
the  pattern  during  its  passage  across  the  magnetometer  array. 

McHenry  and  Clauet  [1987]  calculated  the  distribution  of 
currents  and  fields  corresponding  to  a  single  filament  of  current 
impressed  on  the  ionosphere.  Two  field  aligned  currents  of 
nearly  equal  magnitude  and  opposite  sense  would  generate  a 
distribution  of  magnetic  fields  below  the  ionosphere  close  to  the 
observed  distribution.  The  observed  pattern,  however,  seems 
to  be  distorted  so  that  the  derived  convection  directions  are 
not  exactly  perpendicular  to  the  line  connecting  the  vortex- 
centers  as  would  be  expected.  McHenry  and  Clauer  ignore 
both  ground-induction  and  conductivity  gradients.  It  is  not 
yet  known  if  these  effects  are  the  cause  of  the  distortion  of  the 
observed  perturbation  pattern  or  if  the  distortion  is  due  to  a 
more  complex  distribution  of  the  field-aligned  currents. 

In  Figure  5  solar  wind  data  measured  by  IMP-8  are  shown. 
B;  is  positive  during  the  event  and  the  solar  wind  velocity  is 
around  450  km/sec.  A  sudden  drop  in  the  plasma  density  from 
11  per  cc  to  6  per  cc  occurs  at  0955.  Considering  the  propa¬ 
gation  time  between  IMP-8  and  the  magnetopause  (about  10 
minutes)  fluctuations  in  the  sign  of  a  large  sudden  change 
in  the  value  of  the  positive  B,  and  the  sudden  drop  in  plasma 
density  are  proximate  to  the  observed  vortex  event. 

While  we  show  here  detailed  observations  for  one  event, 
we  have  examined  15  cases  Analysis  has  concentrated  on  rela- 
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Figure  4.  Total  horizontal  magnetic  perturbation  vectors  cor¬ 
responding  to  the  variations  in  Figure  3  have  been  rotated  by 
90°  counterclockwise  and  plotted  every  20  seconds  during  the 
interval  from  10:06:00  to  10:18:40  UT  For  each  time  the  posi¬ 
tion  of  the  vectors  have  been  off  set  to  the  right  by  a  distance 
corresponding  to  80  km  to  account  for  an  assumed  4  km/sec 
westward  motion  (to  the  left)  of  the  pattern 


IMP-8  PLASMA  AND  FIELD  MEASUREMENTS 


Figure  5.  Solar  wind  data  measured  by  IMP-8.  FYom  top  to 
bottom  is  shown  the  interplanetary  magnetic  field  components 
in  GSM  coordinates  Bz,  B„,  and  Bx<  the  solar  wind  velocity, 
and  the  number  density.  The  bottom  trace  shows  the  derived 
solar  wind  dynamic  pressure  plotted  in  units  of  nanopascals. 


tively  quiet  intervals  having  isolated  spike-like  pulsations,  since 
during  disturbed  or  even  moderately  disturbed  conditions  the 
pulsations  are  difficult  to  extract  from  the  background  field 
Events  were  chosen  according  to  this  criterion  as  well  as  ac¬ 
cording  to  the  magnitude  of  the  perturbation.  A  high-pass 
filter  (30  minutes  cut-off)  was  applied  to  the  data  to  extract 
the  events  from  the  background  field. 

The  pulsations  are  not  always  associated  with  a  twin  vor¬ 
tex  pattern.  Our  study  has  shown  that  the  events  can  be 
sorted  into  different  categories  corresponding  to  the  follow¬ 
ing  types.  1)  isolated  twin  vortices,  2)  single  vortices,  and  3) 
multiple  vortices.  Some  of  the  pulsations  could  not  be  iden¬ 
tified  in  a  consistent  way  as  any  of  the  mentioned  types.  We 
are  presently  conducting  a  more  comprehensive  analysis  of  the 
different  types  of  events. 

Discussion 

Irregular  pulsations  are  a  common  phenomenon  in  ground- 
based  magnetometer  data  measured  near  the  polar  cleft.  The 
irregularity,  however,  makes  them  difficult  to  extract  from  the 
usually  noisy  background  field.  During  intervals  with  a  rela¬ 
tively  quiet  background  field  it  is  possible  to  extract  a  relatively 
unique  current  pattern  which  produces  the  observed  perturba¬ 
tions.  Because  of  the  scale  size  of  the  phenomena,  a  relatively 
dense  array  of  stations  near  the  polar  cleft  is  necessary  to  re¬ 
solve  the  current  structures  and  their  motion. 

To  put  this  event  into  a  more  global  context,  we  find  that 
it  is  associated  with  a  smail,  world  wide  SI  (sudden  impulse) 
related  to  the  solar  wind  dynamic  pressure  change.  Stations  on 
the  east  coast  of  Greenland  located  near  local  magnetic  noon 
exhibit  large  variations.  It  is.  however,  not  possible  to  identify 
a  signature  corresponding  to  the  twin  vor'ex  pattern  observed 
on  the  west  coast.  Stations  in  the  afternoon  sector  at  high  lati¬ 
tudes  which  we  examined  show  only  a  small  SI  variation.  This 
suggests  to  us  that  the  observed  traveling  vortex  current  sys¬ 
tem  is  a  result  of  a  generation  near  local  noon  with  subsequent 
westward  motion 


97 


256 


Friis-Christensen  et  &1.:  Ionospheric  Traveling  Convection  Vortices 


Thp  [  '.in  voit'-x  odservation  which  we  describe  is  similar 
to  “nrs  observed  in  the  auroral  zone  (Honisch  and  Gia&s- 
meier.  1986:  Glassn.eier  (private  communication)  1987].  Using 
the  Scandinavian  Magnetometer  Array,  their  events  had  am¬ 
plitudes  of  order  40  nT.  but  only  viewed  a  portion  of  the  2-cell 
svstem  well  below  the  vortex  centers.  They  also  found  west¬ 
ward  (tadward)  motion  at  about  4  km/sec. 

The  magnetic  signature  of  the  event  presented  here  is  sim¬ 
ilar  to  those  previously  reported  by  Lanzerotti  et  ah.  [1986. 
1987].  Our  modeled  current  system  to  explain  these  variations, 
however,  is  quite  different.  Having  only  one  station  Lanzerotti 
et  ah.  [1986]  modeled  the  current  structure  as  a  Saunders- 
Lee  type  of  distribution  moving  poleward.  Using  a  meridional 
chain  of  stations  it  is  possible,  even  without  detailed  model¬ 
ing.  to  test  the  assumption  of  a  poleward  movement,  and  with 
the  additional  east-west  distribution  of  stations  it  is  possible 
to  estimate  the  east-west  movement  of  the  currents.  The  event 
presented  tiere.  and  several  others  which  have  been  examine  ! 
show  a  predominant  east-west  movement  of  the  currents. 

The  two  pre-noon  events  (June  03.  1982  and  October  16. 
1985)  discussed  by  Lanzerotti  et  ah  [1986;  1987]  are  very  sim¬ 
ilar  to  the  evert  discussed  in  the  present  paper  We  have 
examined  the  available  Greenland  magnetic  data  and  found 
them  to  be  consistent  in  all  components  with  a  westward  mov¬ 
ing  twin- vortex.  The  data  are  not  consistent  with  a  poleward 
moving  single  vortex.  Although  Lanzerotti  et  al.  J 1 986;  1987] 
modeled  the  events  with  a  poleward  moving  Hall  current  pat¬ 
tern  they  did  in  tart,  for  the  October  16,  1985  event  find  a 
westward  5  km/sec  motion  in  the  pattern  when  comparing  the 
magnetic  signatures  at  Sondre  Stramfjord  and  Frobisher  Bay. 
Using  conjugate  data  from  the  northern  and  southern  hemi¬ 
sphere  Lanzerotti  et  al.  [1987]  suggested  that  the  October  16. 
1985  event  occurred  on  closetl  field  line#  and  might  be  associ¬ 
ated  with  the  convection  of  an  elongated  plasma  cloud  in  the 
boundary  layer. 

The  studv  of  the  June  28.  1986  event  and  a  reexamina¬ 
tion  of  the  previously  reported  events  suggested  to  be  the 
ionospheric  signatures  of  FTE's.  leads  to  the  conclusion  that 
none  of  the  examined  events  corresponds  to  either  of  the  two 
proposed  models  of  FTE's.  the  Saunders-Lee  type  of  a  sin- 
de  poleward  moving  Hall  current  vortex,  or  the  Southwood 
type  of  a  twin-vortex.  Although  the  event  presented  here  is 
consistent  with  a  twin-vortex  pattern  of  plasma  convection, 
its  direction  of  movement  is  perpendicular  to  the  motion  pre¬ 
dicted  by  the  model.  The  high  velocities  (3-6  km/sec)  and 
a  direction  of  motioL  parallel  to  the  electric  field  impressed 
on  the  ionosphere  by  the  two  oppositely  directed  field-aligned 
currents  means  that  the  event  does  not  correspond  to  a  net 
transport  of  plasma  from  closed  to  open  field-lines  since  this 
would  require  an  electric  field  perpendicular  to  the  direction  of 
motion  to  create  the  necessary  J  x  B  force.  The  event  seems 
more  probably  to  correspond  to  a  motion  of  a  structure  in  the 
magnetopause  boundary  possibly  created  by  either  a  sudden 
change  in  the  solar  wind  density  and/or  by  a  sudden  change  in 
the  reconnection  regions,  caused  by  an  abrupt  change  in  the 
IMF  orienrarion.  Such  a  readjustment  might  propagate  down¬ 
stream  associated  with  field-aligned  currents  mapping  down  to 
the  ionosphere  as  observed  with  the  chain  of  magnetometers. 
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Modeled  Ground  Magnetic  Signatures  of  Flux  Transfer  Events 

Mark  A.  McHenry  and  C.  Robert  Clauer 

STAR  Laboratory,  Stanford  University,  Stanford,  California 


The  magnetic  field  on  the  ground  due  to  a  small  (<  200  km  scale  size)  localized  field-aligned 
current  system  interacting  with  the  ionosphere  is  calculated  in  terms  of  an  integral  over  the 
ionospheric  distribution  of  field-aligned  current.  Two  different  candidate  current  systems  for  flux 
transfer  events  (PTEs)  are  examined.  The  first  system  has  current  flowing  down  the  center  of 
a  cylindrical  flux  tube  with  a  return  current  uniformly  distributed  along  the  outside  edge.  The 
second  system  has  upward  current  on  one  half  of  the  perimeter  of  a  cylindrical  flux  tube  with 
downward  current  on  the  opposite  half.  The  peak  magnetic  field  on  the  ground  is  found  to  differ 
by  a  factor  of  2  between  the  two  systems,  and  the  magnetic  perturbations  are  in  different  directions 
depending  on  the  observer’s  position.  Assuming  the  current  system  moves  in  a  constant  linear 
poleward  direction,  we  predict  the  ground  magnetic  field  versus  time  which  would  be  measured 
by  a  ground  observatory.  Using  FTE  detection  statistics  at  the  magnetopause,  we  estimate  the 
detection  rate  of  FTEs  at  a  ground  station  under  the  dayside  convection  reversal  to  be  between 
12  and  60  min  per  sighting  depending  on  which  FTE  current  system  is  considered. 


1 .  Introduction 

The  transfer  of  energy  from  solar  plasma  to  the  Earth’s 
magnetospheric  plasma  is  a  topic  which  has  received  many 
years  of  study  but  which  still  holds  many  mysteries.  While 
several  mechanisms  have  been  proposed  to  account  for  the 
transport  of  momentum  across  the  magnetopause  from  the 
solar  wind  to  the  magnetosphere,  magnetic  reconnection 
[Dungey,  1961]  has  been  considered  one  of  the  principal 
electrodynamic  mechanisms  to  produce  this  momentum  ex¬ 
change.  Observations  from  the  ISEE  spacecraft  have  re¬ 
cently  suggested  that  this  process  may  occur  in  both  a  quasi- 
steady  form  [Paschmann  et  al.,  1979]  and  also  in  a  more  im¬ 
pulsive  and  sporadic  form  termed  flux  transfer  events  (FTE) 
[Russell  and  Elphtc,  1978;  Sonnerup  et  al.,  1981].  At  times 
the  occurrence  rate  of  FTEs  appears  to  be  sufficient  to  ac¬ 
count  for  a  major  portion  of  the  momentum  transfer  across 
the  magnetopause  [Cowley,  1982;  Goertz  et  al.,  1985). 

The  magnetic  structures  associated  with  FTEs  observed 
at  the  magnetopause  are  approximately  1  Re  in  cross  sec¬ 
tion  and  also  contain  a  unique  twisting  magnetic  field  in¬ 
dicating  the  presence  of  significant  field-aligned  currents. 
Satellite  particle  detectors  also  have  found  that  inside  the 
FTE,  the  plasma  is  a  mixture  of  magnetospheric  and  mag- 
netoeheath  plasma  (Russell  and  Elphic,  1979].  This  plus  the 
observation  that  FTEs  occur  primarily  when  the  interplan¬ 
etary  magnetic  field  (IMF)  is  largely  opposite  to  the  Earth's 
magnetic  field  [Rijnbeek  et  al.,  1984;  Berchem  and  Russell, 
1984]  suggests  that  FTEs  are  the  result  of  localized  recon¬ 
nection  between  interplanetary  and  magnetospheric  mag¬ 
netic  fields  at  the  dayside  magnetopause. 

Because  of  the  potentially  important  role  of  FTEs  in  cou¬ 
pling  energy  to  the  magnetosphere  from  the  solar  wind,  the 
study  of  FTEs  has  received  recent  attention  both  experi¬ 
mentally  and  theoretically  (see  the  recent  review  by  Cowley 
[1982]).  Among  the  areas  to  receive  considerable  interest 
recently  is  a  search  for  the  signatures  of  FTEs  in  the  iono- 
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sphere  [Goertz  et  al.,  1985;  Todd  et  al.,  1986;  Lanzerotti 
et  al.,  1986].  Such  independent  observations,  if  consistent 
with  the  interpretation  of  FTEs  at  the  magnetopause,  would 
also  provide  considerable  strength  to  the  arguments  support¬ 
ing  the  FTE  interpretation,  and  would  also  provide  another 
abundent  source  of  observational  information  regarding  day- 
side  magnetopause  reconnection. 

If  the  interpretation  of  the  FTE  observations  at  the  mag¬ 
netopause  is  correct,  one  would  expect  that  the  stress  ap¬ 
plied  to  the  reconnected  flux  tube  as  it  is  pulled  tailward 
must  be  transmitted  to  the  ionosphere.  This  must  be  accom¬ 
plished  by  closing  field-aligned  currents,  associated  with  the 
moving  flux  tube,  through  the  ionosphere.  The  ionospheric 
signature  which  might  be  expected  from  a  FTE  will  depend 
upon  the  distribution  of  current  associated  with  the  moving 
flux  tube. 

There  are  two  proposed  FTE  current  distributions  which 
have  been  qualitatively  discussed.  The  first  consists  of  a 
central  core  of  field-aligned  current  which  spreads  radially  in 
the  ionosphere.  The  current  closes  with  field-aligned  return 
currents  in  a  diffuse  shell  around  the  central  core  current 
or  through  distant  return  currents  [Sounders  e<  al.,  1984; 
Lee,  1986].  One  should  note  the  symmetry  of  this  system, 
which  is  shown  schematically  in  Figure  la.  A  second  sys¬ 
tem  which  is  asymmetric  has  been  suggested  by  Southwood 
[1985,  1987],  This  system  consists  as  a  first  approximation 
of  two  oppositely  directed  field-aligned  currents  on  oppos¬ 
ing  edges  of  the  FTE  flux  tube.  A  refined  version  distributes 
the  currents  over  the  circumference  of  the  flux  tube.  This  is 
shown  schematically  in  Figure  1  b. 

In  this  paper  we  examine  the  magnetic  effects  of  these 
two  postulated  FTE  current  systems.  This  is  accomplished 
by  first  obtaining  the  Green's  function  solution  for  a  single 
field-aligned  current  filament.  The  solution  for  a  distributed 
FTE  current  system  is  then  obtained  by  numerically  sum¬ 
ming  the  single-filament  solutions.  We  compute  the  ground 
magnetic  perturbation  of  each  FTE  system  and  then,  as¬ 
suming  a  1  km/s  poleward  velocity,  determine  the  expected 
time-varying  field  which  would  be  measured  by  a  ground 
observatory  located  100  km  from  the  ground  track  of  the 
moving  current  system. 
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Saunders  et  ai.  -  Lee  Southwood 


FTE  Current  System  FTE  Current  System 

Fig.  1  (a)  Postulated  [Meanders  et  ai,  1984;  Lee,  1986]  symmet¬ 
ric  FTE  current  system  composed  of  a  central  core  with  return 
current  around  the  outside  edge.  (6)  Postulated  [Southtcood, 
1985,  1987]  asymmetric  FTE  current  system  composed  of  cur¬ 
rent  flowing  along  the  outside  edge.  The  direction  of  the  current 
reverses  along  the  perimeter. 


application  of  Ohm’s  law  and  current  continuity  [ Kamide 
and  Matsushita.  1979]  gives 

V(E  VG>)  =j||(r)sinx  (1) 

where  the  derivatives  are  in  the  X  -  Y  plane,  E  is  the  height- 
integrated  conductivity  tensor,  <I>  is  the  electric  potential, 
j||(r)  is  the  field-aligned  current  density  toward  the  Earth 
and  x  *s  the  angle  of  Earth’s  magnetic  field  line  in  the  X  -  Z 
plane  with  90°  being  vertical.  It  can  easily  be  shown  [Rat- 
cltffe,  1972]  that  E  has  antisymmetric  off-diagonal  elements 
even  when  the  field  lines  are  not  vertical.  This  with  the 
assumption  of  uniform  conductivities  gives 

V24»  =  j11(r)sinX/Ep  (2) 

where  Ep  is  the  height-integrated  Pedersen  conductivity.  To 
have  a  solution  which  has  no  currents  far  away  from  the 
source,  we  require 
a* 

— - •  0  as  r  — *  oc  (3) 

or 

Green’s  Function  Solution 


2.  Calculation  of  the  Ground  Magnetic  Field  of 
a  Small-Scale  Field-Aligned  Current  System 

Formulation  of  Problem 

Our  goal  in  this  section  is  to  find  the  magnetic  field  on 
the  ground  created  when  a  small  (20O-km  scale  size)  current 
system  is  impressed  on  the  ionosphere.  Because  of  the  small 
size  of  the  current  system  and  because  the  fields  drop  off 
quickly  with  distance,  we  can  assume  the  ionosphere  is  flat 
and  parallel  to  the  Earth’s  surface.  Similarly,  we  assume  the 
Earth’s  magnetic  field  lines  are  straight  but  uniformly  tilted 
in  a  local  north/south  plane.  We  also  use  the  assumptions 
commonly  used  to  find  the  effect  of  larger-scale  current  sys¬ 
tems  [Kamide  and  Matsushita,  1979];  the  ionosphere  is  rep¬ 
resented  by  a  two-dimensional  layer  with  height-integrated 
Hall  and  Pedersen  conductivities,  and  only  steady  state  so¬ 
lutions  are  considered. 

We  also  assume  the  height-integrated  conductivity  is  uni¬ 
form  in  a  region  around  the  small-scale  current  system. 
Radar  measurements  around  the  dayside  convection  reversal 
boundary  where  FTEs  are  likely  to  occur  [Robinson  et  ai, 
1984]  have  been  used  to  determine  the  height-integTated  con¬ 
ductivity  and  found  this  to  be  a  good  approximation.  When 
nonuniformity  does  exist,  one  can  expect  the  magnetic  field 
on  the  ground  to  be  higher  than  otherwise  calculated  be¬ 
cause  of  the  loss  of  symmetry  in  the  driven  ionospheric  cur¬ 
rents. 

The  final  approximation  is  to  ignore  any  ground  induction 
effects.  Typically,  they  are  modeled  by  placing  a  perfectly 
conducting  layer  approximately  200  km  below  the  Earth's 
surface  [Kisabetk  and  Rostoker,  1977],  Image  currents  are 
then  placed  below  the  perfectly  conducting  layer  at  the  same 
distance  as  the  true  currents  are  above  the  layer.  It  will  be 
shown  later  that  the  magnetic  field  from  our  small-scale  cur¬ 
rent  patterns  decreases  rapidly  with  distance,  so  observers 
on  the  Earth’s  surface  (on  the  order  of  500  km  away  from 
the  image  sources)  will  to  first  order  not  see  the  induction 
effects. 

We  start  with  a  coordinate  system  where  the  ionosphere 
is  in  the  X  -  Y  plane  and  Z  points  downward.  A  simple 


A  straight-forward  application  of  Green’s  functions 
[Morse  and  Feshback,  1953]  gives  the  potential  as  an  in¬ 
tegral  over  the  source  field-aligned  currents  in  a  manner 
similar  to  the  familiar  solution  of  Poisson’s  equation  for  a 
two-dimensional  distribution  of  charge. 

*(ro)  =  -| /  jXnrh^dA  (4) 

Here  ro  is  some  point  in  the  ionosphere,  r  is  the  distance 
from  ro  to  the  source  point  (where  the  field-aligned  current 
intersects  the  ionosphere),  and  A  is  the  area  containing  the 
source  points.  The  electric  field  is 

and  the  height-integTated  current  density  i  is 

i(ro)  =  £  ■  E(r0)  =  J  J  X  ■  £j„(r)dA  (6) 

If  we  temporarily  put  the  origin  at  a  source  point,  we  find 

£  ■  J  =  ZPf  -  £/,<*>  (7) 

where  is  the  height-integTated  Hall  conductivity  and  f 
and  <{>  are  unit  vectors  which  change  during  the  integration. 
Using  this  expression  in  the  above  current  expression  i(ro), 

i(ro)  =  £  '  E(ro)  =  fgr  j  J (T.pt  - 'Lhd>)1^-dA  (8) 

we  see  the  familar  Pedersen  and  Hall  currents.  To  get  the 
magnetic  field  of  currents  in  the  ionosphere,  we  use  the  Biot- 
Savart  law 

B(r')  =  S//i(ro)x  hiA'  (9) 

where  r'  is  the  location  of  the  observer,  6  is  a  unit  vector 
from  the  observer  on  the  ground  to  the  integration  point  at 
the  ionosphere,  R  is  the  distance  from  the  observer  to  the 
integration  point  at  the  ionosphere,  and  the  area  of  integra¬ 
tion  A'  is  the  entire  ionosphere.  The  field-aligned  current 
which  acts  as  the  point  source  in  the  above  equations  is  a 
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Fig.  2.  Unit  vectors  and  angles  used  in  calculations.  X  is  northward.  V  is  eastward  and  Z  is  down.  The 
background  magnetic  field  is  in  the  X  -  Z  plane.  The  axes  have  been  repeated  for  clarity  and  have  their  origins 
at  the  point  where  the  field-aligned  current  intersects  the  ionosphere 


semi-infinite  line  current,  and  using  Biot-Savart,  its  contri¬ 
bution  can  be  found.  After  integrating,  we  find 


B  =  ®FAC  +  BPed  Arsen  +  BHall 


(10) 


with 


B'«  =  ~TJJ  (I.) 

(I2> 


the  angles  o  and  o0  are  the  same  and  the  magnetic  fields 
of  the  field-aligned  current  cancels  that  of  the  Pedersen  cur¬ 
rent,  leaving  only  the  field  due  to  the  Hall  currents.  This 
demonstrate'  what  has  become  known  a'  Fukushima’s  the¬ 
orem  [Fukushima.  1969],  Consequently,  if  we  assume  the 
Hall  and  Pedersen  conductivities  are  comparable  and  the 
magnetic  field  lines  are  nearly  vertical,  we  can  neglect  con¬ 
tributions  to  the  magnetic  field  from  the  combined  effect  of 
field-aligned  currents  and  Pedersen  currents. 

Discussion  of  Assumptions 


®H«1I 


Ep  4*  J  J  sin  a  r 


(13) 


where  nq  is  the  angle  between  a  line  connecting  the  observer 
with  the  source  point  and  the  background  magnetic  field  line 
(greater  than  90°  if  the  observer  is  below  the  ionosphere). 
a  is  the  vertical  angle  of  the  line  connecting  the  source  and 
observer,  and  r'  is  the  distance  from  the  source  point  to  the 
observer.  The  azimuthal  unit  vector  0O  is  about  the  field- 
aligned  current,  while  0  is  about  the  vertical  line  passing 
through  the  point  where  the  field-aligned  current  intersects 
the  ionosphere.  The  unit  vector  ft  is  along  the  horizontal 
projection  of  r'  pointing  from  the  source  to  the  observer. 
The  magnetic  field  is  in  the  X  -  Z  plane  and  tilted  \  from 
the  X  axis.  In  Figure  2,  a  repeated  set  of  axes  shows  these 
coordinates  and  unit  vectors.  The  integration  is  to  be  car¬ 
ried  out  over  the  area  where  field-aligned  currents  intersect 
the  ionosphere.  Note  that  if  the  field  is  vertical  (x  =  *72). 


A  critical  assumption  which  would  probably  receive  lit¬ 
tle  debate  is  that  field-aligned  currents  on  closed  field  lines 
must  interact  with  the  lower  ionosphere  where  significant 
cross  B  conductivity  exists.  However,  in  recent  experiments. 
Stcnzel  and  Urrutia  [1986]  have  found  that  when  the  elec¬ 
tric  field  across  B  is  large,  instabilities  occur  which  lead  to 
unexpected  currents  across  B.  This  might  occur  in  the  F 
region  with  FTE  current  systems.  If  this  is  the  case,  the  cur¬ 
rents  will  not  come  as  close  to  the  surface  of  the  Earth,  and 
hence  the  ground  magnetometer  signature  would  be  greatly 
reduced.  In  this  paper,  we  will  neglect  this  possibility  and 
assume  field-aligned  currents  close  in  the  E  region. 

Another  assumption  was  the  use  of  a  quasi-static  solu¬ 
tion.  If  the  induction  electric  field  created  by  the  time- 
varying  magnetic  field  is  large  compared  to  the  impressed 
elect ru  field  m  the  ionosphere,  perhaps  a  different  pattern 
of  currents  will  exist  in  the  ionosphere.  To  estimate  these 
induction  electric  fields,  we  use  Faraday  s  law.  a  scale  size  of 
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Fig.  3.  Magnetic  field  strength  on  the  ground  due  to  a  single  field-aligned  current  filament  of  2  x  10s  A  interacting 
with  the  ionosphere.  The  X  axis  measures  ground  distance  in  the  north  direction,  and  the  Y  axis  measures  ground 
distance  in  the  east  direction.  Each  contour  represents  10  nT,  with  the  dashed  contours  for  negative  values.  The 
ground  magnetic  field  perturbations  are  shown  separately  for  (top  row)  the  field-aligned  current  portion,  (middle 
row)  the  ionospheric  Pedersen  current,  and  (bottom  row)  the  ionospheric  Hall  current  Each  row  shows  contours 
of  magnetic  perturbation  in  the  X  (magnetic  northward),  Y  (eastward),  and  Z  (downward)  components  of  the 
magnetic  held  as  a  function  of  ground  distance  from  the  field-aligned  current  filament.  The  background  magnetic 
field  is  in  the  X  -  Z  plane. 


100  km,  and  a  magnetic  field  change  of  100  nT  in  60  s.  This 
results  in  an  electric  field  strength  of  0.2  mV/m.  A  rough 
estimate  of  the  impressed  electric  field  can  be  made  using 
equation  (5)  assuming  a  current  system  of  two  very  small 
filaments  separated  by  100  km  with  current  in  the  filaments 
in  opposite  directions. 


E(  r0)  = 


sin  x 

2rEp 


1 


r2 


(14) 


Here  /  is  the  current  in  each  filament,  and  r\  and  r%  are 
the  distances  from  where  each  filament  intersects  the  iono¬ 
sphere  to  some  other  point  in  the  ionosphere.  If  x  — 
90°.  =  o  mhos  I  —  2  x  105  A,  and  the  net  dis¬ 

tance  term  is  50  km,  the  impressed  electric  field  would  be 
60  mV/m.  which  is  in  the  range  postulated  by  Lee  [1986]. 
Hence  induction  electric  fields  are  negligible  for  current 
changes  taking  place  in  the  order  of  minutes. 
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Results  for  a  Single  Field-Aligned  Current  Filament 


For  a  single  field-aligned  current  filament,  j(|  (r)  becomes  a 
delta  function,  and  the  magnetic  field  below  the  ionosphere 

due  to  this  perturbation  is 

BFAC--r(1Tco8ao)M, 

4*  r'  sin  ao 

(15) 

„  no  (1  +  cosa)  rl 

oP«d«fMn  —  a  » -  ^  *  $ 

4w  r'  sin  a 

(16) 

T>  _  MO  r/(l  +  cosa) 

B«*»-  Ep4xr'/(  sina  P  +  *> 

(17) 

where  I  is  the  total  current  in  the  filament.  Figure  3  shows 
contour  plots  of  the  components  of  the  magnetic  field  in 
the  X ,  Y  and  Z  magnetic  coordinate  directions  for  each 
constituent  part  of  the  current.  Here  /  is  2  x  10s  A  cor¬ 
responding  to  what  a  FTG  might  carry  [Saunders  et  al., 
1984].  The  current  is  directed  toward  the  Earth  and  does 
not  close  locally  but  through  the  ionosphere  at  r  —  oc.  Here 
Eh/Ep  =  2  and  the  conducting  layer  is  assumed  be  100  km 
above  the  ground.  The  choice  of  positive  Efc  is  for  the  north¬ 
ern  hemisphere.  In  the  southern  hemisphere  the  field  due  to 
the  Hall  currents  would  be  opposite  to  the  direction  shown. 
If  the  current  in  the  filament  was  directed  away  from  the 
Earth,  all  the  fields  would  be  reversed.  The  magnetic  field 
lines  are  tilted  10°  from  the  vertical  toward  the  south  or 
negative  X  direction  corresponding  to  the  edge  of  the  polar 
cap.  The  top  row  of  Figure  3  shows  plots  of  the  magnetic 
field  on  the  ground  due  to  the  field-aligned  current  (FAC) 
alone.  Contours  are  drawn  at  10-nT  intervals,  and  negative 
values  are  dashed.  Note  the  asymmetry  due  to  the  tilt  of  the 
field  line.  The  middle  row  shows  the  magnetic  field  due  to 
the  ionospheric  Pedersen  currents.  As  described  above,  the 
Pedersen  currents  produce  a  field  which  cancels  that  of  the 
field-aligned  current  when  the  magnetic  field  is  vertical,  and 
here  we  see  the  fields  are  almost  identical  and  opposite  in 
direction.  In  the  single-filament  case  shown  in  Figure  3,  the 
magnetic  field  due  to  the  Pedersen  currents  and  FAC  cur¬ 
rents  together  is  less  than  10-nT  in  most  areas,  with  a  peak 
value  around  20-nT.  Their  effect  on  a  FTE  current  system 
of  similar  strength  would  be  only  several  nanoteslas  because 
of  the  large  cancellation  due  to  local  closure  of  the  current. 
Hence  in  the  remaining  portion  of  this  paper,  we  will  assume 
the  contributions  of  the  field-aligned  current  and  Pedersen 
currents  cancel  exactly,  and  only  the  Hall  current's  contri¬ 
bution  to  the  ground  magnetic  field  will  be  considered. 

The  bottom  row  of  Figure  3  shows  the  field  strength  on 
the  ground  due  to  the  Hall  currents  generated  by  the  single 
field-aligned  current  filament.  The  Hall  currents  are  a  set  of 
circular  loops  in  a  horizontal  plane  about  the  point  where  the 
field-aligned  current  intersects  the  ionosphere.  They  act  like 
a  soleniod,  creating  a  strong  vertical  field.  Because  the  sole- 
niod  is  flat,  the  field  lines  are  not  confined  and  spread  out¬ 
ward  away  from  the  vertical  axis.  This  creates  the  outward 
radial  field  shown  in  the  first  two  panel?  of  the  bottom  row. 
Underneath  the  source  point  we  see  a  peak  field  strength 
of  400  nT  representing  the  maximum  field  strength  for  a 
given  amount  of  field-aligned  current.  Distributed  current 
systems  or  current  systems  with  local  field-aligned  closure 
(positive  and  negative  j||(r))  having  the  same  size  currents 


will  have  a  ground  signature  reduced  depending  on  the  exact 
distribution  of  current  density. 

From  the  1-Re  size  of  FTEs  at  the  magnetopause  and 
conservation  of  magnetic  flux  between  the  magnetopause 
and  the  ionosphere  it  is  believed  FTEs  should  have  current 
systems  on  the  order  of  200  km  [ Cowley ,  1984;  Goertz  et  al., 
1985)  across  at  the  ionosphere.  The  size  of  the  Green's  func¬ 
tion  or  field-aligned  current  delta  function  response  shown 
in  Figure  3  is  roughly  100  km,  indicating  that  the  current 
distribution  which  a  FTE  has  at  the  ionosphere  will  strongly 
influence  the  ground  magnetic  signature  and  is  as  important 
as  the  peak  current  strength. 

3.  Results  for  Postulated  Flux 
Transfer  Event  Current  Systems 

Currently,  there  are  two  FTE  current  systems  discussed 
in  the  literature.  One  possibility  [South wood,  1985,  1987) 
is  that  the  field-aligned  currents  form  a  cylinder  with  one 
side  directed  upward  and  the  other  side  directed  downward 
as  shown  in  Figure  16.  Southwood  suggests  the  asymme¬ 
try  is  required  to  have  Pedersen  currents  create  a  J  x  B 
force  which  pushes  the  flux  tube  through  the  lower  colli¬ 
sions!  ionosphere.  Hence  the  dividing  line  between  upward 
and  downward  currents  is  parallel  to  the  direction  of  travel  of 
the  disturbance.  As  a  first  approximation,  we  assume  this 
dividing  line  is  in  the  north/south  direction  and  the  FTE 
moves  in  a  poleward  direction.  With  these  assumptions, 
7||  (r)  at  the  ionosphere  becomes 

7n(r)  =  ^*(r  -  °)s>n(^)  (18) 

where  /  is  the  integral  of  the  upward  (or  downward)  current, 
a  is  the  radius  of  the  current  system  and  u>  is  the  azimuthal 
angle  around  the  inside  of  the  current  system.  (Motion  other 
than  poleward  would  introduce  a  sin(u/  +  u>o)  dependence.) 
The  poleward  edge  of  the  current  system  would  have  w  equal 
to  zero.  The  sin(u>)  dependence  and  circular  shape  were  cho¬ 
sen  for  simplicity.  The  top  row  of  Figure  4  shows  the  ground 
magnetic  ‘footprint’  of  this  current  system  using  equation 
(13)  when  a  =100  km,  /  =  2  x  105  A,  E/,/Ep  =  2,  and  the 
conducting  layer  is  100  km  above  the  ground.  The  horizon¬ 
tal  Y  axis  measures  ground  distance  in  the  east  direction, 
and  the  vertical  X  axis  measures  ground  distance  in  the 
north  direction. 

A  more  symmetric  current  system  [Saunders  et  al.,  1984; 
Lee,  1986]  is  composed  of  a  central  core  current  with  an  out¬ 
side  return  current  shown  in  Figure  la.  The  central  core  cur¬ 
rent  provides  an  azimuthal  magnetic  field  which  when  added 
to  the  background  field  creates  a  twisted  field  agreeing  with 
in  situ  satellite  observations.  The  exact  radial  distribution 
of  field-aligned  current  in  this  model  is  not  known,  but  for 
simplicity  we  consider  a  simple  delta  function  representation 

J||(r)  =  -6(r  -  a)  -  -/>(r)  (19) 

"a  t 

where  /  is  again  the  total  upward  (or  do- n ward)  current. 
The  direction  of  the  currents  could  be  reversed  depending 
on  the  Y  component  of  the  IMF  [Lee.  1986).  The  ‘footprint' 
of  this  case  is  shown  in  the  bottom  row  of  Figure  4  with  / 
and  a  the  same  as  the  Southwood  system.  Here  the  value 
of  the  vertical  field  on  the  ground  directly  below  the  center 
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Fig.  4.  (Top)  Magnetic  field  strength  on  the  ground  due  to  a  FTE  assuming  the  Southwood  current  configuration 
with  a  radius  of  100  km  at  the  ionosphere.  The  X  axis  measures  ground  distance  in  the  north  direction,  and  the 
Y  axis  measures  ground  distance  in  the  east  direction.  The  net  upward  (or  downward)  current  is  equal  to  2  x  10® 
A.  The  contours  are  dashed  for  negative  values  and  spaced  20  nT.  (Bottom)  Magnetic  field  strength  on  the  ground 
due  to  a  FTE  assuming  the  Saunders  et  al.-Lee  current  configuration  with  a  radius  of  100  km  at  the  ionosphere. 
The  net  upward  (or  downward)  current  is  equal  to  2  x  10s  A. 


of  the  system  can  be  found  by  direct  integration  of  equation 
(13)  using  the  current  density  of  equation  (19). 


R  _  _ j _ 1 

*  EP  V'vSi+TS*  h 

If  h  >>  a,  this  becomes 


(20) 


B, 


Ea  no  .  a2 
Ep  4ir  2 P 


(21) 


which  is  the  result  of  [Laruerottt  et  at.,  1986].  However,  the 
assumption  of  h  >>  a  is  not  valid  in  this  situation  since  h 
and  a  are  comparable.  Using  equation  (20),  we  find  a  peak 
value  of  B ,  of  120  nT  with  E^/Ep  =  2,  a  =  A  =100  km 
and  /  =  2  x  10®  A  while  equation  (21)  yields  200  nT  for  B, 
in  the  same  situation.  The  large  difference  between  these 
estimates  is  not  unexpected,  because  the  Hall  current  is  a 
distributed  current  system  which  is  quite  different  from  a 
current  loop  when  the  observer  on  the  ground  is  relatively 
cloee. 

To  find  the  magnetic  field  at  an  arbitrary  location  caused 
by  the  above  FTE  field-aligned  current  models,  equation 
(13)  must  be  numerically  integrated  using  a  known  field- 
aligned  current  system.  To  obtained  the  top  view  of  Fig¬ 


ure  4,  Southwood’s  field-aligned  current  system  was  used 
(equation  (18)),  while  for  the  bottom  row  of  Figure  4,  the 
Saunders  et  al.-Lee  system  (equation  (19))  was  used.  Each 
contour  line  in  these  two  figures  represents  a  change  of  20 
nT,  and  negative  values  are  shown  with  dashed  contours. 
Comparing  the  footprints  of  the  two  proposed  FTE  current 
systems,  we  see  they  differ  not  only  in  peak  value  but  also 
in  bask  configuration.  The  magnetic  perturbation  produced 
by  the  cylindrical^  symmetric  system  proposed  by  Saunders 
et  al.  and  Lee  has  an  area  less  than  half  that  produced  by 
the  system  suggested  by  Southwood.  It  also  has  a  signifi¬ 
cantly  smaller  peak  value  and  a  smaller  horizontal  magnetic 
field  at  the  ground  as  opposed  to  Southwood’s  strong  hor¬ 
izontal  component.  Note  that  symmetry  requires  the  hori¬ 
zontal  magnetic  field  point  toward  or  away  from  the  center 
of  the  Saunders  et  al.-Lee  current  system  unlike  the  South- 
wood  system,  where  the  horizontal  field  points  in  a  roughly 
uniform  direction  at  all  points  on  the  ground.  The  Z  com¬ 
ponent  of  the  magnetic  field  of  the  Southwood  system  also 
has  regions  of  strong  negative  and  strong  positive  values, 
while  the  Saunders  et  al.-Lee  system  has  no  sign  change  in 
the  Z  component.  If  the  FTE  footprint  moved  in  a  direc¬ 
tion  other  than  poleward,  the  Southwood  patterns  would  be 
rotated  corresponding  to  the  angle  away  from  poleward  the 
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Fig.  5.  Ionospheric  convection  c rented  by  the  Southwood  model 
of  FTE  field-aligned  currents.  The  net  downward  current  and 
upward  current  is  2  x  10*  A.  All  the  field-aligned  currents  intersect 
the  ionosphere  along  a  100- km  radius  circle.  The  scale  vector  has 
length  corresponding  to  a  velocity  of  3500  m/s. 


FTE  footprint  moved.  The  Saunders  et  al.-Lee  perturba¬ 
tions  would  change  sign  if  the  direction  of  the  field-aligned 
currents  were  reversed. 

Perhaps  a  more  realistic  version  of  the  symmetric  Saun¬ 
ders  et  al.-Lee  field-aligned  current  system  would  be  to  re¬ 
place  the  delta  function  core  with  a  uniform  current  distribu¬ 
tion  avoiding  singularities  in  the  electric  field  in  the  center. 
Although  not  shown,  the  ground  magnetic  field  of  this  sys¬ 
tem  is  similar  to  the  delta  function  core  version  shown  in 
the  bottom  row  of  Figure  4  with  the  horizontal  field  purely 
radial  and  with  the  peak  Z  perturbation  underneath  the 
center.  The  strength  of  the  field,  however,  has  a  slightly 
different  radial  dependence.  As  done  in  equation  (20),  the 
vertical  field  directly  under  the  center  of  the  current  system 
can  be  calculated. 


B, 


_  Eft  PO  r( 

Ep4)riVi 
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-2 


y/a1  +  h2-h> 


(22) 


Va’  +  A3 

Assuming  a  ~  A,  this  gives  a  vertical  field  strength  di¬ 
rectly  under  the  current  system  which  is  approximately  half 
that  of  the  system  with  a  delta  function  core.  The  large  sen¬ 
sitivity  of  the  ground  magnetic  field  due  to  the  field-aligned 
current  configuration  is  due  to  the  size  of  the  current  sys¬ 
tem  being  roughly  the  same  as  the  height  of  the  ionosphere 
and  the  large  amount  of  cancellation  when  calculating  the 
magnetic  field  when  the  field-aligned  currents  close  locally. 


Combination  of  Postulated  Flux  Transfer  Event  Current 
Systems 

The  field-aligned  current  system  associated  with  a  FTE 
must  be  consistent  with  the  flux  tube  translating  in  the 
ionosphere  and  with  twist  in  the  magnetic  field  inside  the 
tube.  Since  the  two  systems  studied  here  each  produce  only 
one  of  the  desired  effects,  it  is  likely  that  a  combination  of 


them  describe  the  true  FTE  situation.  First,  we  require  t  hat 
as  the  combined  field-aligned  current  system  moves  across 
the  ionosphere,  flux  must  remain  inside  the  flux  tube,  and 
plasma  outside  must  be  moved  out  of  the  way.  To  find  the 
ionospheric  convection  pattern  of  the  Southwood  system  of 
field-aligned  currents  of  equation  (18)  with  /  =  2  x  10*  A, 
we  assume  E  x  B  drift  and  use  equation  (5)  to  get  E.  The 
result  is  a  large  twin-vortex  pattern  as  shown  in  Figure  5. 
Inside  the  100-km  radius  ionospheric  footprint  of  the  FTE 
flux  tube,  the  flow  is  uniform  with  northward  velocity  ( vx )■ 


vz 


I  sin2  x 
4a  Ep  B 


(23) 


If  the  field-aligned  current  system  does  not  move,  plasma 
in  the  ionosphere  circulates  with  no  net  motion,  since  the 
field-aligned  currents  cannot  apply  a  net  force  in  a  uniform 
ionosphere.  If  we  now  assume  the  field-aligned  current  sys¬ 
tem  translates  and  the  plasma  inside  the  flux  tube  moves 
with  the  field-aligned  current  system,  then  equation  (23)  re¬ 
lates  the  amount  of  current  in  the  Southwood  system  to  the 
velocity  of  the  field-aligned  current  system.  Adding  field- 
aligned  current  in  a  cylindrically  symmetric  system  such  as 
the  Saunders  et  al.-Lee  system,  with  its  radial  convection 
pattern,  does  not  cause  flux  to  leave  the  inside  of  the  FTE. 
Second,  the  observed  twist  in  the  magnetic  field  inside  a 
FTE  near  the  magnetopause  adds  the  requirement  that  the 
Saunders  et  al.  Lee  portion  of  the  field-aligned  current  be 
much  larger  than  the  Southwood  portion.  We  now  can  use 
equation  (18)  along  with  the  observed  speed  of  FTEs  (of  or¬ 
der  1  km/s  [ Goertz  et  al.,  1985;  Lanzerotti  et  al.,  1986;  Todd 
et  al.,  1986])  at  the  ionosphere  to  estimate  the  amount  of 
field-aligned  current  of  a  FTE  assuming  a  Southwood  con¬ 
figuration  with  100-km  radius.  Using  a  vertical  magnetic 
field  of  0.6  G  with  *  =  80°  and  Ep  =  2.5  mhos,  the  net 
downward  field-aligned  current  (/)  would  be  0.6  x  10®  A. 
Assuming  this  current  system  was  added  to  a  Saunders  et 
al.-Lee  system  of  2  x  10*  A,  a  twisted  magnetic  field  would 
be  expected  inside  the  FTE,  although  it  would  be  distorted. 
The  contributions  of  the  two  current  systems  to  the  ground 
magnetic  field,  however,  would  be  comparable  even  though 
the  Southwood  current  strength  is  a  factor  of  3  smaller  be¬ 
cause  of  the  small  ground  perturbation  of  the  Saunders  et 
al.-Lee  system. 

For  the  remaining  analysis,  however,  the  Southwood  and 
Saunders  et  al.-Lee  current  systems  will  continue  to  be 
considered  separately  because  estimation  of  the  relative 
amounts  of  current  in  each  system  in  a  combined  model 
depends  on  the  total  current,  the  conductivities,  and  the 
ionospheric  speed  of  FTEs,  which  are  not  well  known  and 
possibly  variable. 


Ground  Magnetometer  Signatures  of  Postulated  Flux 
Transfer  Event  Current  Systems 

To  predict  the  time  record  of  a  ground  magnetometer 
when  a  FTE  passes  near  by.  the  current  strength  as  a  func¬ 
tion  of  time  must  be  known  as  well  as  how  fast  and  in  what 
direction  the  disturbance  propagates.  If  we  assume  a  FTE's 
field-aligned  current  system  is  stable  and  moves  in  a  straight 
line  across  the  ionosphere,  its  ground  magnetic  signature  at 
a  fixed  observatory  can  be  approximated  by  a  cross  section 
of  the  ground  magnetic  ‘footprint'  taken  along  the  overhead 
track.  Figure  6  shows  for  each  current  system  the  compo- 
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Fig.  6.  (Left)  Magnetic  field  versus  time  of  a  FTE  assuming  the  Saunders  et  al.  Lee  current  system  moving 
poleward  at  a  rate  of  1  km/s  and  passing  100  km  ground  distance  away  to  the  east.  The  strength  and  configuration 
of  the  field-aligned  currents  are  assumed  to  be  constant.  (Right)  Magnetic  field  versus  time  of  a  FTE  with  the 
Southwood  current  system  moving  poleward  at  a  rate  of  1  km/s  and  passing  100  km  ground  distance  away  to  the 
east.  The  strength  and  configuration  of  the  field-aligned  currents  are  assumed  to  be  constant. 


nents  of  the  magnetic  field  on  the  ground  which  would  be 
measured  for  a  FTE  moving  poleward  at  1  km/s  and  with  its 
center  passing  100  km  to  the  east  of  a  ground  magnetometer 
located  in  the  northern  hemisphere.  The  central  core  current 
in  the  Saunders  et  al.-Lee  system  was  directed  away  from 
the  Earth.  Note  the  large  differences  in  the  basic  shape  of 
each  component  and  how  the  asymmetric  Southwood  cur¬ 
rent  model  has  more  than  twice  the  peak  magnitude  of  the 
symmetric  Saunders  et  al.-Lee  model  for  equal  field  aligned 
currents. 

Unfortunately,  there  have  only  been  a  few  possible  ob¬ 
servations  of  FTEs  in  the  ionosphere  to  verify  the  above 
assumption  of  a  stable  FTE  current  system  moving  across 
the  ionosphere  with  constant  velocity  [Goertz  et  al.,  1985; 
Todd  et  al,  1986;  Lanzerotti  et  al,  1986].  These  experiments 
do  indicate  that  the  ionospheric  signature  of  a  FTE  moves 
with  a  speed  of  0.5  to  2.0  km/s  and  has  a  duration  on  the 
order  of  4  min.  While  this  4-min  duration  could  represent 
the  time  it  took  for  the  FTE  to  cross  the  field  of  view  of 
the  experiment,  it  is  possible  it  could  represent  the  lifetime 


of  the  ionospheric  signature.  Comparing  this  4-min  lifetime 
with  the  time  it  takes  a  FTE  to  pass  a  ground  station  as 
noted  in  Figure  6,  we  see  the  assumption  of  constant  field- 
aligned  current  is  only  marginally  met.  The  true  ground 
magnetometer  signature  would  be  similiar  to  that  in  Figure 
6  during  the  time  period  with  significant  field-aligned  cur¬ 
rents  and  greatly  reduced  during  the  time  intervals  when 
the  FTE  field-aligned  currents  cease. 

One  way  to  identify  a  FTE  using  its  ground  magnetic 
signature  without  assumptions  of  its  ionospheric  motion  or 
field-aligned  current  strength  is  to  use  simultaneous  mea¬ 
surements  using  an  array  of  closely  spaced  magnetometers. 
The  observed  horizontal  magnetic  perturbation  vector  plot¬ 
ted  at  the  location  of  the  station  can  then  be  easily  compared 
to  theoretical  models.  In  Figure  7,  the  horizontal  magnetic 
perturbation  is  shown  for  both  the  Southwood  and  Saunders 
et  al.-Lee  models.  The  length  of  the  scale  vector  shown  is  200 
nT  in  the  Southwood  plot  but  only  100  nT  in  the  Saunders 
et  al.-Lee  data.  Note  the  the  large  Y  perturbation  normal 
to  the  assumed  X  -directed  velocity  in  the  Southwood  model 
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Fig.  7.  (Top)  Horizontal  magnetic  perturbation  due  to  a  FTE 
assuming  the  Southwood  field-aligned  current  system.  The  length 
of  the  scale  vector  represents  a  200-nT  perturbation  The  field 
is  largely  perpendicular  to  the  assumed  X-directed  motion  of  the 
FTE.  (Bottom)  Horizontal  magnetic  perturbation  due  to  a  FTE 
assuming  the  Saunders  et  al.-Lee  field-aligned  current  system. 
The  length  of  the  scale  vector  represents  a  100- nT  perturbation. 
The  field  is  radial  about  the  center  of  the  field-aligned  current 
system  because  of  symmetry. 


while  the  magnetic  perturbation  is  radial  from  the  center  of 
the  Saunders  et  al.-Lee  system.  As  mentioned  above,  the 
true  FTE  signature  will  be  a  combination  of  the  Southwood 
and  Saunders  et  al.-Lee  models,  with  the  ratio  depending 
on  the  conductivities,  the  FTE  ionospheric  speed,  and  the 
total  amount  of  current. 


Estimate  of  Ground  Magnetometer  FTE  Detection  Rate 

An  important  factor  in  observing  the  ground  signature  of 
a  FTE  is  the  detection  rate.  A  statistical  survey  [ Rtjnbeek 
et  al.,  1984]  of  satellite  detections  near  the  dayside  magne¬ 
topause  has  shown  that  FTEs  exist  uniformly  over  a  large 
region.  At  a  given  location  in  this  region.  FTEs  recur  every 
7  to  8  min  on  average.  These  data  are  consistent  with  the 
idea  that  the  observed  FTEs  are  caused  by  the  antisunward 
motion  of  reconnected  flux  tubes  moving  along  the  magne¬ 
topause  carried  by  magnetosheath  flow. 

A  fundamental  unanswered  question  regarding  FTEs  is 
the  global  generation  rate.  A  variety  of  assumptions  can  be 
made  to  estimate  this  rate.  For  example,  let  us  consider 
only  the  northern  hemisphere  and  assume  that  each  FTE  is 
initially  formed  near  the  equatorial  magnetopause.  Further 
assume  that  the  reconnected  flux  tube  is  pulled  tailward 
across  the  magnetopause  either  toward  dawn  or  dusk  de¬ 
pending  on  the  orientation  of  the  IMF  and  the  reconnection 
location.  A  satellite  located  at  the  dayside  magnetopause 
should  then  detect  a  large  fraction  of  the  antisunward  mov¬ 
ing  FTE  flux  tubes.  If  we  assume  that  on  average  half  the  re¬ 
connected  tubes  move  toward  the  dawn  flank  and  half  move 
toward  the  dusk  flanks,  we  might  expect  a  satellite  to  detect 
at  most  half  of  the  globally  generated  FTEs.  If  FTEs  are 
generated  globally  at  least  every  4  min,  what  would  be  the 
expected  detection  rate  of  a  ground  magnetometer? 

To  find  the  ionospheric  density  of  FTEs,  their  starting 
location  and  ground  track  must  be  known  and  we  can  only 
speculate  on  these  at  present.  Estimating  that  FTEs  can  be 
formed  on  the  front  third  of  the  magnetosphere,  we  map  this 
region  to  roughly  the  dayside  third  of  the  auroral  oval.  This 
is  a  band  perhaps  3000  km  wide  in  longitude.  The  normal 
dimensions  of  FTEs  at  the  magnetopause  when  mapped  to 
the  auroral  oval  give  a  width  of  several  hundred  kilometers, 
producing  a  detection  band  which  is  3000  km  in  longitudi¬ 
nal  length  and  about  300  km  in  latitudinal  width  located 
around  the  convection  reversal  boundary.  This  band  forms 
the  most  probable  area  where  the  ionospheric  signature  of  a 
FTE  would  start. 

If  we  define  50  nT  to  be  the  minimum  ground  magnetic 
field  perturbation  required  to  detect  a  FTE,  we  find  from 
Figure  4  that  the  Southwood  configuration  produces  a  600- 
km  detection  range,  while  the  Saunders  et  al.-Lee  configu¬ 
ration  produces  a  100-km  detection  range.  Thus,  roughly 
three  observatories  suitably  distributed  over  the  detection 
band  would  be  sufficient  to  observe  the  global  generation  of 
FTEs  of  the  Southwood  model,  while  15  stations  would  be 
required  for  the  Saunders  et  al.-Lee  model.  If  the  ground 
track  of  FTEs  is  composed  of  a  significant  longitudinal  com¬ 
ponent,  fewer  stations  may  be  required.  Of  course,  more 
stations  distributed  in  latitude  will  be  required,  since  the 
detection  band  will  move  in  latitude  as  the  auroral  oval  ex¬ 
pands  and  contracts. 

It  follows  from  the  above  estimates  that  if  a  station  was 
known  to  be  under  the  dayside  convection  reversal  bound¬ 
ary  and  the  IMF  was  southward,  FTE  signatures  should  be 
observed  in  the  data  roughly  every  12  (60)  min  on  the  aver¬ 
age  for  the  Southwood  (Saunders  et  al.-Lee)  current  model. 
This  represents  an  upper  limit  on  the  time  between  detec¬ 
tions  because  of  the  conservative  global  FTE  generation  rate 
assumed  above.  It  is  also  a  long-term  average,  since  the  lo¬ 
cation  of  reconnection  could  remain  steady  for  some  period 
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of  time;  hence  the  ionospheric  signatures  of  the  accompa¬ 
nying  FTEs  could  follow  nearly  the  same  path  repeatedly, 
allowing  a  higher  detection  rate  for  a  given  ground  station. 

4.  Summary  and  Conclusions 

Making  the  assumptions  of  a  thin  conducting  region  in 
the  ionosphere,  near-uniform  conductivity,  and  small  induc¬ 
tion  electric  fields,  the  ground  magnetic  field  of  an  arbitrary 
small-scale  field-aligned  current  system  has  been  reduced 
to  a  simple  integral  over  the  field-aligned  current  filaments. 
This  approximation  includes  the  tensor  nature  of  the  con¬ 
ductivity  of  the  ionosphere.  These  results  have  been  applied 
separately  to  two  models  of  FTE  field-aligned  current  sys¬ 
tems.  The  resulting  ground  magnetic  fields  are  of  the  order 
of  100  to  150  nT  depending  on  the  size  and  strength  of  the 
current  systems  at  the  ionosphere.  The  two  models  have 
significantly  different  vertical  and  horizontal  magnetic  fields 
on  the  ground  as  indicated  in  Figure  4.  To  meet  the  re¬ 
quirements  of  twist  in  the  magnetic  field  inside  the  FTE 
near  the  magnetopause  and  movement  of  the  FTE  footprint 
in  the  ionosphere,  the  most  likely  field-aligned  current  ar¬ 
rangement  is  a  combination  of  the  Southwood  and  Saunders 
et  al.-Lee  models.  By  comparing  the  calculated  velocity  of 
the  ionospheric  convection  to  the  observed  FTE  speed,  we 
find  the  fraction  of  total  current  in  the  Southwood  model 
roughly  a  factor  of  3  less  than  in  the  Saunders  et  al.-Lee 
model,  but  because  the  Southwood  system  has  a  much  larger 
magnetic  signature,  the  ground  magnetic  effects  are  compa¬ 
rable. 

From  the  above  detection  rate  estimates,  it  appears  that 
using  ground  magnetometers  is  a  plausible  method  to  study 
FTEs  if  an  area  underneath  the  convection  reversal  could 
be  monitored  with  a  chain  of  magnetometers.  It  is  likely 
that  the  relatively  strong  magnetic  signature  produced  by 
the  Southwood  asymmetric  configuration  would  be  easily 
detected.  If,  however,  FTE  currents  are  predominantly  of 
the  Saunders  et  al.-Lee  symmetric  configuration,  the  smaller 
magnetic  effects  are  much  less  likely  to  be  detected,  and  a 
more  dense  magnetometer  array  will  be  required. 

If  FTEs  could  be  reliably  identified  from  the  ground,  not 
only  could  the  current  configuration  and  time  history  be 
studied  using  the  results  of  this  paper,  but  the  field  line 
where  reconnection  started  could  be  found  as  well  as  check¬ 
ing  the  above  estimates  in  the  total  worldwide  FTE  rate. 
Since  a  single  FTE  does  not  transfer  enough  flux  to  explain 
the  existing  magnetospheric  potential  [Rijnbeek  et  ai.  1984], 
their  rate  of  occurrence  is  critical  in  deciding  their  overall 
importance. 
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The  pattern  of  ionization  and  electric  Helds  measured  in  the  davside  auroral  zone 
by  the  Sondrestrom  incoherent  scatter  radar  has  been  compared  with  interplanetary 
magnetic  field  (IMF)data  from  the  IMP-8  spacecraft.  The  electric  field  configuration 
is  most  sensitive  to  the  B ,  component  of  the  IMF.  For  8,  <0  there  is  a  strong  reversal 
in  convective  flow  between  1300  and  1600  MLT  with  a  region  of  enhanced  F  region 
ionization  associated  with  low  energy  particle  precipitation  at  the  reversal.  For  8,  >0 
the  convection  reversal  is  most  apparent  on  the  morning  side.  Ionization  produced  by 
particle  precipitation  for  8,  >0  is  much  less  conspicuous  than  that  observed  on  the 
afternoon  side  for  8,  <0.  We  show  radar  measurements  of  ionization,  electric  fields 
and  electron  temperatures  associated  with  two  distinct  types  of  particle  precipitation 
in  the  pre-noon  local  time  sector  for  8.  >0.  The  patterns  of  precipitation  for  8,  <0and 
8.  >0  are  consistent  with  the  distribution  of  field-aligned  currents,  in  that 
precipitation  occurs  primarily  in  the  upward  current  regions.  With  this  interpretation, 
only  the  precipitation  associated  with  the  cusp  field-aligned  currents  on  the  dawn  side 
may  be  intense  enough  to  produce  ionospheric  signatures  observable  from  the 
ground. 

I.  Introduction 

The  large-scale  pattern  of  electric  fields  and  field-aligned  currents  in  the  auroral 
zone  and  its  response  to  changes  in  the  interplanetary  magnetic  field  have  been 
studied  extensively  in  past  years  (see,  for  example.  HeppnER,  1972;  Friis- 
CHRISTENSEN  et al.s  1972;  llJlMAand  Potemra,  1978;  Heelis.  1984).  Bl'RCH  el  al. 
(1985)  have  shown  that  the  electric  field  pattern  can  be  understood  in  terms  of  the 
superposition  of  several  cells  referred  to  as  either  merging,  viscous  or  lobe  cells.  Each 
cell  responds  differently  to  changes  in  the  interplanetary  magnetic  field.  Precipitation 
is  expected  in  various  parts  of  the  convection  pattern  depending  on  the  connection  of 
these  field  lines  to  various  plasma  regimes. 

Since  operations  began  in  Greenland  in  February.  1983.  numerous  experiments 
were  performed  in  which  the  Sondrestrom  radar  was  used  to  measure  the  distribution 
of  ionization  and  electric  fields  at  very  high  latitudes.  The  radar  data  provide  a 
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complimentary  view  to  that  obtained  with  polar  orbiting  satellites,  because  patterns 
of  ionization  and  electric  fields  can  be  constructed  from  measurements  made  at 
different  local  times.  However,  the  interpretation  of  these  observed  patterns  must  be 
done  carefully  because  electric  fields  may  change  substantially  during  the  time  over 
which  the  measurements  are  made.  Interplanetary  magnetic  field  data  are  presently 
available  for  many  of  the  radar  experiments  so  that  variations  of  electric  fields  in 
response  to  IMF  conditions  can  be  identified. 

In  this  paper  we  describe  the  results  of  a  study  in  which  we  examined  radar 
measurements  in  the  davside  auroral  oval  and  polar  cap  made  during  24  experiments. 
After  describing  the  data  used  in  the  study,  we  present  a  short,  qualitative  description 
of  the  observed  electric  field  patterns  and  their  dependence  on  the  IMF.  We  then 
discuss  the  ionization  that  occurs  in  association  with  these  patterns  and  look  in  detail 
at  two  examples  of  ionization  enhancements  produced  by  particle  precipitation  in  the 
morning  sector. 

2.  Data  Base 

The  radar  data  used  in  this  study  were  obtained  during  experiments  on  24  days 
chosen  on  the  basis  of  the  availability  of  IMF  data.  Some  of  these  experiment:,  were 
performed  as  part  of  the  World  Day  observations;  others  were  performed  as  part  of 
campaigns  to  study  electrodynamics  in  the  polar  cusp.  Preliminary  results  of  these 
experiments  have  been  described  in  a  number  of  papers  (see.  for  example, 
JORGENSEN  el  at. .  1984.  Ci.AI'ER  el  al. ,  1984;  Kelly  and  Vickrey,  1984).  For  the 
present  study,  we  examined  the  data  set  for  the  specific  purpose  of  identifying  IMF 
B,  -dependent  patterns  of  particle  precipitation  and  electric  fields  in  the  prenoon  and 
postnoon  local  time  sectors.  The  radar  operating  mode  for  most  of  the  24  days 
consisted  of  some  combination  of  fixed-position  measurements  and  scans.  The  fixed- 
position  measurements  at  a  sequence  of  elevation  angles  and  azimuths  provided 
electric  field  data  over  a  broad  latitudinal  range,  while  the  scans  yielded  high  time-and 
spatial-resolution  measurements  of  ionization  and  drifts  over  a  more  limited  region. 

Table  1  shows  the  number  of  times  different  orientations  of  S,  occurred  during 
radar  measurements  in  the  morning  and  afternoon  sectors.  All  IMF  components 
given  in  this  paper  refer  to  the  GSM  coordinate  system.  To  be  counted  as  an 
occurrence.  B,  had  to  maintain  the  same  sign  for  longer  than  an  hour,  which  is  greater 
than  the  15  minutes  to  one-half  hour  required  for  electric  fields  on  the  dayside  to 


T able  I  Occurrence  frequency  of  various  orientations  of  I M  F  B.  in  morning  and  afternoon  sectors  for 
the  data  set  used  in  this  study 
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respond  to  changes  in  the  IMF  at  the  location  of  IMP-8  (CLALER  et  at .,  1984).  No 
sorting  according  to  the  magnitude  of  B,  was  performed  and  the  data  include  cases  of 
both  positive  and  negative  values  of  B..  Although  Bt  was  highly  variable,  there  were 
many  days  when  this  component  of  the  IMF  remained  fairly  constant  over  long  time 
intervals.  Out  of  the  24  days  there  were  eight  in  which  B,  was  positive  throughout  the 
morning  and  afternoon  sectors,  while  there  were  five  in  which  B,  remained  negative 
during  this  local  time  interval. 

3.  Large-Scale  Electric  Fields  and  Their  Relation  to  IMF  B, 

In  general,  the  electric  fields  measured  by  the  radar  during  the  local  times 
between  0600  and  1800  are  quite  consistent  with  the  results  of  previous  studies  based 
on  satellite  data.  Of  particular  note  is  the  strong  dependence  of  the  electric  field  on  the 
y  component  of  the  IMF  (DE  l.A  BEAUJARDIERE  et  al..  this  volume).  For  By< 0  the 
radar  measures  predominantly  eastward  drifts  in  the  morning  sector,  while  in  the 
evening  sector  both  eastward  and  westward  convecting  plasma  is  observed.  The 
westward  (sunward)  drifting  plasma  is  confined  to  the  equatorward  part  of  the  oval, 
with  the  eastward  (antisunward)  drifting  plasma  poleward  of  it.  The  antisunward 
drifting  plasma  is  typically  quite  strong  and  confined  to  the  region  immediately 
poleward  of  the  convection  reversal.  As  pointed  out  by  JORGENSEN  et  al.  ( 1984)  the 
convection  reversal  is  typically  rotational;  however,  the  change  in  the  east-west 
component  of  the  velocity  is  particularly  large  and  is  thus  an  easily  identifiable  feature 
in  latitudinal  plots  of  electric  field  in  the  post-noon  sector  (Robinson  et  al..  1984). 

For  fl,> 0  the  convection  reversal  in  the  afternoon  sector  disappears  and 
westw  ard  drifts  dominate  at  all  latitudes.  A  corresponding  reversal  appears,  however, 
on  the  morningside  with,  again,  sunward  drifts  on  the  equatorward  side  and 
antisunward  drifts  on  the  poleward  side.  The  reversal  in  velocity  on  the  morning  side 
is  not  as  dramatic  as  on  the  afternoon  side,  but  is  still  quite  obvious  in  the  radar  data. 

The  convection  observed  near  noon  is  highly  variable,  but  tends  to  be  dominated 
by  westward  drifts  for  B,  >0  and  eastward  drifts  for  B <0.  Occasionally,  antisunward 
drifting  plasma  at  noon  can  be  observed,  especially  when  Bz  is  large  and  negative. 

These  general  statements  about  the  patterns  of  electric  fields  and  their 
dependence  on  B,  underscore  the  importance  of  IMF  data  in  the  interpretation  of 
radar  measurements.  Under  conditions  when  B,  is  changing  rapidly,  radar  measure¬ 
ments  through  the  prenoon  to  postnoon  local  time  sector  reveal  multiple  reversals  in 
the  zonal  convection  of  plasma.  Fortunately,  there  are  enough  periods  of  fairly  stable 
IMF  conditions  to  allow  unambiguous  mapping  of  the  local  time  and  latitudinal 
distribution  of  electric  fields  from  radar  measurements.  These  studies  are  currently  in 
progress. 

4.  Ionization  Produced  by  Particle  Precipitation  and  Its  Dependence  on  IMF  B, 

Because  much  of  the  electron  precipitation  that  occurs  in  the  dayside  oval  and 
polar  cap  is  very  soft,  ionization  is  predominantl>  produced  at  altitudes  above  200 
km.  At  these  altitudes,  ionization  has  lifetimes  from  several  tens  of  minutes  to  hours 
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and  can  be  transported  great  distances  from  the  location  at  which  it  is  produced. 
Thus,  the  appearance  of  ionization  enhancements  does  not  necessarily  indicate  the 
presence  of  local  production  by  particle  precipitation.  Although  there  are  other 
measureable  indicators  of  local  production  that  can  be  used,  we  show  here  only 
ionization  that  can  be  definitely  attributed  to  particle  precipitation,  either  because  the 
precipitation  was  hard  enough  to  ionize  at  altitudes  below  200  km,  or  because 
simultaneous  satellite  data  were  available  to  confirm  the  association.  Because  we 
limited  ourselves  to  these  obvious  indicators  of  particle  precipitation,  we  may  have 
overlooked  regions  of  precipitation  in  which  the  energy  of  the  particles  was  too  low  to 
produce  measureable  enhancements  in  ionization.  Kofman  and  WickwaR  (1984) 
and  WICKWAR  (this  volume)  describe  a  method  based  on  solving  the  electron  heat 
balance  equation  by  which  regions  of  local  production  can  be  identified  regardless  of 
the  accompanying  ionization.  This  technique  may  prove  useful  for  future  studies. 

The  most  obvious  region  of  particle  precipitation  observed  with  the  radar  at 
Sondrestrom  is  that  associated  with  the  convection  reversal  in  the  afternoon  sector 
when  £, <0  (Fostfr.  et  al .,  1983;  Claier  ct  at.,  1983).  Rorisson  ei  ai  (1984) 
presented  radar  and  satellite  measurements  in  this  precipitation  region  to  show  that  it  is 
associated  with  the  region  1  Birkeland  current  and  precipitating  electrons  accelerated 
through  parallel  potential  drops  of  about  1  kV  from  a  source  of  density  3.5  cm  '  and 
temperature  350  eV.  Evans  (1984)  studied  precipitation  in  this  local  time  sector  using 
data  from  the  NOAA  satellites  and  confirmed  the  persistence  of  this  feature  even  in 
periods  of  very  low  magnetic  acitivity.  The  ionospheric  signature  and  induced  effects 
of  this  precipitation  have  been  examined  by  Wit  kwar  (this  volume).  For  B,< 0  on 
the  morning  side  no  localized  region  of  enhanced  ionization  was  identified  in  the 
radar  data  that  could  be  attributed  to  particle  precipitation. 

For  B,> 0,  a  reversal  in  the  convective  flow  of  plasma  appears  on  the  morning 
side.  In  examining  the  fourteen  periods  of  measurements  during  which  B ,  was  positive 
in  the  morning  hours  we  were  able  to  identify  two  examples  of  ionization  produced  by 
particle  precipitation.  Figure  I  shows  data  obtained  on  13  July  1983  while  the  radar 
was  performing  an  elevation  scan  in  the  plane  of  the  magnetic  meridian.  At  the  time  of 
these  measurements,  fl,  was  + 10  nT  and  B-.  was  -  10  nT.  The  top  panel  in  Fig.  I  is  a 
contour  plot  of  electron  density  in  the  meridian  plane  constructed  from  the  radar 
measurements.  The  center  of  the  plot  corresponds  to  the  location  of  the  radar  which  is 
at  74  degrees  invariant  latitude.  Directly  above  the  radar  a  region  of  ionization  was 
observed  with  a  peak  density  of  about  3x  (O'  cm  ’  at  an  altitude  of  200  km.  Because 
this  region  of  ionization  extends  downward  to  altitudes  well  below  200  km.  we  can  be 
reasonably  sure  that  local  production  by  precipitating  particles  was  occurring  at  this 
location.  The  second  panel  in  Fig.  I  shows  the  latitudinal  profile  of  electric  field 
measured  during  the  scan.  The  eastward  electric  field  is  small  to  the  south  and 
increases  gradually  to  the  north.  The  meridional  electric  field  is  southward  with  a 
peak  value  of  about  75  mV  mjust  poleward  of  the  ionization  enhancement.  Assuming 
no  significant  gradients  in  the  east-west  electric  field  in  the  zonal  direction  and 
uniform  ionospheric  conductivity,  the  north-south  gradients  in  the  meridional 
electric  field  are  consistent  with  the  presence  of  upward  field-aligned  currents  within 
the  ionization  enhancement  and  downward  field-aligned  currents  poleward  of  the 
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DISTANCE  (km) 

Fig.  I  Radar  data  obtained  during  a  5-minute  scan  in  (he  magnetic  meridian  on  13  July.  1983  Distances 
are  in  km  magnetically  north  and  south  from  the  radar 


enhancement.  The  sense  of  these  currents  is  consistent  with  that  of  the  region  1  and 
region  2  Birkeland  currents  (llJlM  A  and  PoTEMRA,  1976).  The  ionization  feature  at  74 
degrees  latitude  was  apparent  at  the  same  location  in  the  data  from  the  subsequent 
radar  scan  (not  shown)  which  was  made  about  a  half  hour  later.  Thus,  the 
precipitation  producing  this  ionization  was  probably  extended  at  constant  invariant 
latitude  along  the  auroral  oval  over  about  one-half  hour  in  local  time  and  persisted  for 
at  least  this  long  in  Universal  Time.  The  latitudinal  structure  of  electric  field  observed 
during  the  second  scan  was  somewhat  different.  The  precipitation  was  associated  with 
a  local  perturbation  in  a  wide  region  of  fairly  uniform  southward  fields.  However,  the 
gradient  in  electric  field  associated  with  the  perturbation  was,  as  before,  consistent 
with  the  presence  of  upward  field-aligned  currents.  In  summary,  the  observations  on 
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13  July  are  consistent  with  arc-like  precipitation  embedded  in  sunward  drifting 
plasma  and  associated  with  upward  field  aligned  currents. 

The  second  type  of  precipitation  identified  in  the  morning  sector  for  positive 
I M  F  B>  is  shown  in  Plate  1  The  top  panel  shows  a  color  plot  of  electron  density  in  the 
magnetic  meridian  plane  observed  on  13  April  1983  at  about  1000  magnetic  local 
time.  At  this  time  B  was  about  +8  nT  and  B.  was  about  -5  nT.  The  radar  data  show 
several  regions  of  latitudinally  confined  ionization  enhancements  with  electron 
densities  of  about  3  x  10  cm  ’.  Because  these  enhancements  were  situated  at  altitudes 
of  about  400  km.  they  cannot  definitively  be  associated  with  particle  precipitation. 
Electron  temperature  measurements  from  the  same  scan  are  shown  in  the  second 
panel.  These  data  show  that  the  electron  temperature  between  73  and  74  degrees 
invariant  latitude  was  higher  than  that  in  the  surrounding  plasma,  reaching 
tempertures  of  over  4000  K.  This  suggests  local  production  of  heat  and  ionization  by 
precipitating  particles.  This  conclusion  was  confirmed  by  data  from  the  NOAA-6 
satellite  during  a  north  polar  pass  between  1 137  and  1 154  UT.  Radar  and  satellite  data 
at  this  time  are  shown  in  Fig.  2.  The  top  panel  is  a  plot  of  the  total  energy  flux  carried 
into  the  atmosphere  by  precipitating  ions  and  electrons.  These  data  show  a  region  of 
enhanced  precipitation  within  the  same  latitude  range  as  the  ionization  enhancement 
detected  by  the  radar.  The  peak  energy  flux  exceeded  10  ergs  cm2-s  and  was  carried 
almost  exclusively  by  precipitating  electrons.  Spectral  information  about  the 
precipitating  electrons  indicated  average  energies  less  than  about  300  eV.  Electric 
fields  measured  by  the  radar  at  about  the  time  of  the  NOAA-6  pass  are  shown  in  the 
bottom  panel  of  Fig.  2.  Latitudinal  profiles  of  the  northward  electric  field  at  two 
different  times  are  shown.  The  similarity  between  the  two  sets  of  electric  field 
measurements  confirms  the  constancy  of  the  electric  field  pattern  over  the  duration  of 
the  observations.  The  electric  field  over  most  of  the  region  scanned  by  the  radar  was 
strong  and  northward  indicating  antisunward  drifting  plasma.  The  particle 
precipitation  is  located  just  poleward  of  where  the  electric  field  maximizes.  The 
gradient  in  the  electric  field  assuming  uniform  ionospheric  conductivity  is  consistent 
with  the  presence  of  upward  field-aligned  currents.  The  location  of  these  currents  is 
well  within  the  polar  cap  as  determined  by  the  direction  of  the  plasma  drift.  Thus,  the 
ionization  shown  in  Plate  1  was  probably  the  ionospheric  signature  of  a  polar  cap  arc. 
Note  that,  as  pointed  out  by  Wickwar  and  KOFMAN  (1984),  the  high  electron 
temperatures  along  with  the  enhanced  ionization  can  lead  to  substantial  thermal 
excitation  of  the  O' D  state  of  atomic  oxygen,  so  that  630.0  nm  emissions  probably 
exceeded  a  few  kR  in  intensity  from  this  source  alone. 

5.  Discussion 

By  examining  Sondrestrom  radar  measurements  made  during  the  morning  and 
afternoon  hours,  we  have  attempted  to  relate  ionization  produced  by  particle 
precipitation  to  boundaries  of  the  convective  flow  and  to  they  component  of  the  IMF. 
The  most  obvious  signature  of  electron  precipitation  is  the  F  region  ionization 
enhancement  associated  with  the  afternoon  sector  convection  reversal  during  times 
when  IMF  B>  is  negative.  Evans  (1984)  showed  that  this  precipitation  probably 
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INVARIANT  LATITUDE 

Fig  1.  NOAA-6  measurements  of  total  energy  flux  and  Sondrestrom  radar  measurements  of  the 
northward  electric  field  obtained  on  13  April  1983  at  about  1000  magnetic  local  time. 


sector  that  are  obviously  produced  by  discrete  particle  precipitation  are  rare.  In 
examining  the  radar  data  for  the  22  mornings  during  which  iMF  data  were  available, 
we  found  only  two  examples  of  ionization  produced  by  particle  precipitation.  The 
absence  of  such  ionospheric  enhancements  at  this  local  time  may  be  due  to  several 
effects.  First,  as  noted  above,  the  precipitating  electrons  may  be  very  low  in  energy, 
thus  depositing  most  of  their  energy  at  high  altitudes  where  ionization  production  is 
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slow.  Second,  the  precipitation  may  be  weak  and  diffuse  so  that  it  is  indistinguishable 
from  ionization  produced  by  solar  illumination  which  is  a  significant  source  at  these 
local  times.  Finally,  the  discrete  precipitation  at  these  local  times  may  occur  outside  of 
the  radar  field  of  view. 

In  this  paper  we  have  shown  two  cases  of  ionization  in  the  morning  sector 
traceable  to  discrete  particle  precipitation.  Both  occurred  between  1000  and  1200 
M  LT  during  times  when  S,  was  positive  and  B:  was  negative.  One  was  observed  in  the 
region  of  sunward  convecting  plasma  while  the  other  was  observed  in  the  region  of 
antisunward  convecting  plasma. 

While  the  event  in  the  sunward  convecting  plasma  can  probably  be  identified  as  a 
morning  sector  auroral  arc,  the  precipitation  in  the  region  of  antisunward  moving 
plasma  must  have  been  situated  within  the  polar  cap.  The  presence  of  discrete  aurora 
at  very  high  latitudes  in  the  morning  sector  has  been  noted  by  Lassen  and 
Damelsen  ( 1978)  and  Gussenhoven  (1982).  That  these  features  are  not  observed 
more  often  from  Sondrestrom  may  be  a  result  of  the  size  of  the  auroral  oval  relative  to 
the  latitude  range  within  which  the  Sondrestrom  radar  can  make  measurements.  On 
April  13  the  A.',,  index  for  the  three-hour  period  from  0900  to  1200  UT  was  6  so  that  the 
oval  was  probably  greatly  expanded  at  this  time.  For  most  of  the  other  mornings 
included  in  our  data  set,  Kr  was  between  2  and  4.  When  a  convection  reversal  was 
present  on  these  days  (B,  >0)  it  was  usually  located  between  74°  and  76°  invariant 
latitude.  Thus,  the  radar  field  of  view  only  extended  a  few  hundred  km  into  the  polar 
cap.  On  13  April,  the  convection  reversal  was  located  at  71°  invariant  latitude  and  a 
greater  portion  of  the  polar  cap  was  observed. 

In  trying  to  account  for  the  precipitation  features  that  have  been  observed  from 
Sondrestrom,  it  is  significant  to  note  that  all  are  associated  with  gradients  in  the 
electric  field  that  are  consistent  with  upward  field-aligned  currents.  In  the  case  of  the 
afternoon  precipitation,  this  corresponded  to  the  region  1  Birkeland  current.  For  the 
precipitation  occurring  in  the  sunward  drifting  plasma  on  the  morning  side,  we 
associated  the  upward  current  with  region  2  Birkeland  current.  The  precipitation 
observed  on  antisunward  drifting  field  lines  probably  coincided  with  the  polar  cusp 
currents.  These  currents  have  been  studied  by  Iijima  etal.  ( 1 978)  and  McDiARMlDez 
al.  (1978)  and  were  found  to  depend  on  the  B,  component  of  the  IMF. 

With  these  considerations  we  can  predict  the  locations  at  which  intense  particle 
precipitation  is  likely  to  occur  based  on  the  convection  patterns  that  have  been 
observed.  Figure  3  shows  the  general  pattern  of  convection  for  B,< 0  and  B,> 0. 
Upward  field-aligned  currents  occur  where  the  divergence  of  the  electric  field  is 
negative  (Lyons,  1978).  These  regions  on  the  dayside  are  shown  by  shading  in  Fig.  3. 
Each  of  the  regions  of  particle  precipitation  thus  far  identified  from  the  radar 
measurements  can  be  associated  with  one  of  these  shaded  areas.  In  the  nomenclature 
of  BURCH  et  al.  ( 1 985),  the  precipitation  at  the  convection  reversal  in  the  afternoon  for 
B,  < 0  and  equatorward  of  the  convection  reversal  in  the  morning  for  B,  >0  are  both 
associated  with  the  viscous  cell.  The  poleward  region  of  precipitation  for  B,  >0  on  the 
morning  side  is  associated  with  the  merging  cell 

It  should  be  noted  that  although  we  did  not  sort  our  data  by  B,,  this  component 
of  the  IMF  was  never  large  and  positive  for  a  sustained  period  for  any  of  the  data  sets 
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Fig.  3.  Relationship  between  observed  precipitation  and  the  large-scale  convective  flow  of  plasma.  The 
shaded  regions  show  the  locations  on  the  dayside  where  ionospheric  effects  of  particle  precipitation 
were  observed 


and  there  were  no  examples  in  which  the  convection  at  very  high  latitude  was 
sunward.  Thus,  the  type  of  polar  cap  arcs  that  appear  under  sustained  B:  north 
conditions  were  not  observed  on  these  days.  On  the  other  hand,  the  data  shown  for 
April  1 3  clearly  indicates  precipitation  in  the  polar  cap.  The  distinction  between  these 
two  types  of  polar  cap  arcs  may  be  that  the  former  is  associated  with  the  lobe 
convection  cell,  while  the  latter  is  associated  with  the  merging  cell  (Burch  el  al., 
1985). 

6.  Conclusion 

Many  claims  have  been  made  about  identification  of  polar  cusp  precipitation 
from  ground  observations.  An  identification  of  a  precipitation  feature  as  originating 
from  the  polar  cusp  cannot  be  made  solely  on  the  proximity  of  the  feature  to  local 
noon.  One  should  instead  consider  the  possible  sources  of  the  precipitation.  The 
precipitation  associated  with  the  post-noon  convection  reversal  boundary  appears  to 
originate  from  the  low  latitude  boundary  layer,  and  thus  is  not  polar  cusp 
precipitation.  However,  the  precipitation  observed  on  13  April  in  a  region  of 
antisunward  convecting  plasma  is  probably  associated  with  the  merging  cell  (Burch 
etal.,  1985).  For  6,  <0,  a  corresponding  region  of  antisunward  drifting  plasmaexists 
on  the  afternoon  side  of  the  polar  cap.  However,  as  the  polar  cusp  currents  are 
downward  in  this  region,  we  should  not  expect  the  precipitation  to  be  as  intense  as  on 
the  dawnside  for  B,  >0.  Because  the  polar  cusp  field-aligned  currents  change  direction 
on  opposite  sides  of  the  polar  cap.  it  follows  that  the  fluxes  of  precipitating  electrons 
associated  with  the  polar  cusp  must  exhibit  a  corresponding  asymmetry.  This  may 
account  for  the  results  of  studies  of  the  distribution  of  polar  cap  arcs  which  show  that 
these  precipitation  features  appear  more  often  on  the  dawn  side  in  the  northern 
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hemisphere  (see,  for  example,  Lassen  and  Danielsen,  1978). 
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